comparison normalize.r @ 1:710627b47962 draft

Uploaded
author ynewton
date Wed, 26 Sep 2012 19:09:23 -0400
parents 31cfcab40d8f
children 600872152be6
comparison
equal deleted inserted replaced
0:31cfcab40d8f 1:710627b47962
124 normal_mean <- mean(normal_values) 124 normal_mean <- mean(normal_values)
125 normal_sd <- sd(normal_values) 125 normal_sd <- sd(normal_values)
126 normalized_normals <- unlist(lapply(normal_values, function(x){return((x - normal_mean)/normal_sd);})) 126 normalized_normals <- unlist(lapply(normal_values, function(x){return((x - normal_mean)/normal_sd);}))
127 127
128 tumor_values <- data_row[tumors_list] 128 tumor_values <- data_row[tumors_list]
129 tumor_mean <- mean(tumor_values) 129 normalized_tumors <- unlist(lapply(tumor_values, function(x){return((x - normal_mean)/normal_sd);}))
130 tumor_sd <- sd(tumor_values)
131 normalized_tumors <- unlist(lapply(tumor_values, function(x){return((x - tumor_mean)/tumor_sd);}))
132 130
133 return(append(normalized_normals, normalized_tumors)) 131 return(append(normalized_normals, normalized_tumors))
134 } 132 }
135 } 133 }
136 134
145 return(apply(ranked_data_matrix, c(1,2), qexp)) 143 return(apply(ranked_data_matrix, c(1,2), qexp))
146 } 144 }
147 else if(dist == 'NORMAL'){ 145 else if(dist == 'NORMAL'){
148 ranked_data_matrix <- apply(data_matrix,2,rankNA) 146 ranked_data_matrix <- apply(data_matrix,2,rankNA)
149 #return(apply(ranked_data_matrix, c(1,2), function(x){return(qnorm(mean=mean(x), sd=sd(x)));})) 147 #return(apply(ranked_data_matrix, c(1,2), function(x){return(qnorm(mean=mean(x), sd=sd(x)));}))
150 return(apply(ranked_data_matrix, c(1,2), qnorm, mean=0, sd=2)) 148 return(apply(ranked_data_matrix, c(1,2), qnorm, mean=0, sd=1))
151 } 149 }
152 else if(dist == 'WEIBULL_0.5'){ 150 else if(dist == 'WEIBULL_0.5'){
153 ranked_data_matrix <- apply(data_matrix,2,rankNA) 151 ranked_data_matrix <- apply(data_matrix,2,rankNA)
154 return(apply(ranked_data_matrix, c(1,2), qweibull, scale=1, shape=0.5)) 152 return(apply(ranked_data_matrix, c(1,2), qweibull, scale=1, shape=0.5))
155 } 153 }