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 } |
