Mercurial > repos > ynewton > matrix_normalization
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 } |