R/getCorchild.R

In jaspershen/MSannotator: xMSannotator: R package for network based annotation of metabolomics data.

getCorchild <-
function(cur_mzdata,data_mt,cor.method){
	
	
	
	pearson_res<-lapply(1:dim(data_mt)[2],function(j){
		 return(WGCNA::cor(as.numeric(cur_mzdata),data_mt[,j],method=cor.method,use="pairwise.complete.obs"))
		
	
	})
	
	pvalues_list<-{}
	pearson_resmat<-{}
	pearson_list<-{}
	
	for(i in 1:length(pearson_res))
	{
		pearson_list<-c(pearson_list,pearson_res[[i]][1])
		nval<-length(which(is.na(data_mt[,i])==FALSE))
		
		pvalues_list<-c(pvalues_list,docortest(nval,pearson_res[[i]][[1]]))
	}
	
		
	return(list(cormat=pearson_list,complete_pearsonpvalue_mat=pvalues_list)) #,complete_pearsonqvalue_mat=qvalues_list))
}