function corr_pval=vassarp(vector_r,num_cells) %corr_pval=vassarp(vector_r,num_cells) % 2017 A.Leo tdist=nan(length(vector_r),1); corr_pval=nan(length(vector_r),1); tdist=vector_r./(sqrt(((1-(vector_r.^2))./(num_cells-2)))); corr_pval=1-tcdf(tdist,num_cells-1); end