Mercurial > repos > bgruening > music_inspect_eset

annotate scripts/estimateprops.R.orig @ 5:2b58bd6becc7 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
author bgruening
date Mon, 28 Oct 2024 17:32:42 +0000
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
5
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
1 suppressWarnings(suppressPackageStartupMessages(library(xbioc)))
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
2 suppressWarnings(suppressPackageStartupMessages(library(MuSiC)))
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
3 suppressWarnings(suppressPackageStartupMessages(library(reshape2)))
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
4 suppressWarnings(suppressPackageStartupMessages(library(cowplot)))
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
5 ## We use this script to estimate the effectiveness of proportion methods
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
6
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
7 ## Load Conf
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
8 args <- commandArgs(trailingOnly = TRUE)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
9 source(args[1])
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
10
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
11 ## Estimate cell type proportions
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
12 est_prop <- music_prop(
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
13 bulk.eset = bulk_eset, sc.eset = scrna_eset,
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
14 clusters = celltypes_label,
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
15 samples = samples_label, select.ct = celltypes, verbose = T)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
16
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
17
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
18 estimated_music_props <- est_prop$Est.prop.weighted
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
19 estimated_nnls_props <- est_prop$Est.prop.allgene
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
20 ##
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
21 estimated_music_props_flat <- melt(estimated_music_props)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
22 estimated_nnls_props_flat <- melt(estimated_nnls_props)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
23
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
24 scale_yaxes <- function(gplot, value) {
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
25 if (is.na(value)) {
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
26 gplot
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
27 } else {
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
28 gplot + scale_y_continuous(lim = c(0, value))
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
29 }
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
30 }
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
31
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
32 sieve_data <- function(func, music_data, nnls_data) {
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
33 if (func == "list") {
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
34 res <- list(if ("MuSiC" %in% methods) music_data else NULL,
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
35 if ("NNLS" %in% methods) nnls_data else NULL)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
36 res[lengths(res) > 0] ## filter out NULL elements
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
37 } else if (func == "rbind") {
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
38 rbind(if ("MuSiC" %in% methods) music_data else NULL,
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
39 if ("NNLS" %in% methods) nnls_data else NULL)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
40 } else if (func == "c") {
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
41 c(if ("MuSiC" %in% methods) music_data else NULL,
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
42 if ("NNLS" %in% methods) nnls_data else NULL)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
43 }
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
44 }
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
45
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
46
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
47 ## Show different in estimation methods
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
48 ## Jitter plot of estimated cell type proportions
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
49 jitter_fig <- scale_yaxes(Jitter_Est(
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
50 sieve_data("list",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
51 data.matrix(estimated_music_props),
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
52 data.matrix(estimated_nnls_props)),
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
53 method.name = methods, title = "Jitter plot of Est Proportions",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
54 size = 2, alpha = 0.7) + theme_minimal(), maxyscale)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
55
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
56 ## Make a Plot
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
57 ## A more sophisticated jitter plot is provided as below. We separated
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
58 ## the T2D subjects and normal subjects by their disease factor levels.
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
59 m_prop <- sieve_data("rbind",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
60 estimated_music_props_flat,
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
61 estimated_nnls_props_flat)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
62 colnames(m_prop) <- c("Sub", "CellType", "Prop")
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
63
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
64 if (is.null(celltypes)) {
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
65 celltypes <- levels(m_prop$CellType)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
66 message("No celltypes declared, using:")
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
67 message(celltypes)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
68 }
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
69
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
70 if (is.null(phenotype_factors)) {
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
71 phenotype_factors <- colnames(pData(bulk_eset))
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
72 }
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
73 ## filter out unwanted factors like "sampleID" and "subjectName"
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
74 phenotype_factors <- phenotype_factors[
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
75 !(phenotype_factors %in% phenotype_factors_always_exclude)]
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
76 message("Phenotype Factors to use:")
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
77 message(paste0(phenotype_factors, collapse = ", "))
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
78
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
79 m_prop$CellType <- factor(m_prop$CellType, levels = celltypes) # nolint
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
80 m_prop$Method <- factor(rep(methods, each = nrow(estimated_music_props_flat)), # nolint
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
81 levels = methods)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
82
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
83 if (use_disease_factor) {
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
84
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
85 if (phenotype_target_threshold == -99) {
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
86 phenotype_target_threshold <- -Inf
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
87 message("phenotype target threshold set to -Inf")
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
88 }
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
89 ## the "2" here is to do with the sample groups, not number of methods
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
90 m_prop$Disease_factor <- rep(bulk_eset[[phenotype_target]], 2 * length(celltypes)) # nolint
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
91 m_prop <- m_prop[!is.na(m_prop$Disease_factor), ]
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
92 ## Generate a TRUE/FALSE table of Normal == 1 and Disease == 2
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
93 sample_groups <- c("Normal", sample_disease_group)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
94 m_prop$Disease <- factor(sample_groups[(m_prop$Disease_factor > phenotype_target_threshold) + 1], # nolint
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
95 levels = sample_groups)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
96
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
97 ## Binary to scale: e.g. TRUE / 5 = 0.2
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
98 m_prop$D <- (m_prop$Disease == # nolint
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
99 sample_disease_group) / sample_disease_group_scale
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
100 ## NA's are not included in the comparison below
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
101 m_prop <- rbind(subset(m_prop, Disease != sample_disease_group),
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
102 subset(m_prop, Disease == sample_disease_group))
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
103
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
104 jitter_new <- scale_yaxes(
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
105 ggplot(m_prop, aes(Method, Prop)) +
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
106 geom_point(aes(fill = Method, color = Disease,
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
107 stroke = D, shape = Disease),
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
108 size = 2, alpha = 0.7,
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
109 position = position_jitter(width = 0.25, height = 0)) +
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
110 facet_wrap(~ CellType, scales = "free") +
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
111 scale_colour_manual(values = c("white", "gray20")) +
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
112 scale_shape_manual(values = c(21, 24)) + theme_minimal(), maxyscale)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
113
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
114 }
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
115
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
116 if (use_disease_factor) {
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
117
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
118 ## Plot to compare method effectiveness
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
119 ## Create dataframe for beta cell proportions and Disease_factor levels
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
120 ## - Ugly code. Essentially, doubles the cell type proportions for each
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
121 ## set of MuSiC and NNLS methods
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
122 m_prop_ana <- data.frame(
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
123 pData(bulk_eset)[rep(1:nrow(estimated_music_props), length(methods)), #nolint
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
124 phenotype_factors],
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
125 ## get proportions of target cell type
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
126 ct.prop = sieve_data("c",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
127 estimated_music_props[, phenotype_scrna_target],
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
128 estimated_nnls_props[, phenotype_scrna_target]),
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
129 ##
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
130 Method = factor(rep(methods,
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
131 each = nrow(estimated_music_props)),
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
132 levels = methods))
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
133 ## - fix headers
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
134 colnames(m_prop_ana)[1:length(phenotype_factors)] <- phenotype_factors #nolint
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
135 ## - drop NA for target phenotype (e.g. hba1c)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
136 m_prop_ana <- subset(m_prop_ana, !is.na(m_prop_ana[phenotype_target]))
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
137 m_prop_ana$Disease <- factor( # nolint
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
138 ## - Here we set Normal/Disease assignments across the methods
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
139 sample_groups[(
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
140 m_prop_ana[phenotype_target] > phenotype_target_threshold) + 1
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
141 ],
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
142 sample_groups)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
143 ## - Then we scale this binary assignment to a plotable factor
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
144 m_prop_ana$D <- (m_prop_ana$Disease == # nolint
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
145 sample_disease_group) / sample_disease_group_scale
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
146
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
147 jitt_compare <- scale_yaxes(
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
148 ggplot(m_prop_ana, aes_string(phenotype_target, "ct.prop")) +
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
149 geom_smooth(method = "lm", se = FALSE, col = "black", lwd = 0.25) +
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
150 geom_point(aes(fill = Method, color = Disease,
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
151 stroke = D, shape = Disease),
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
152 size = 2, alpha = 0.7) + facet_wrap(~ Method) +
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
153 ggtitle(paste0(toupper(phenotype_target), " vs. ",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
154 toupper(phenotype_scrna_target),
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
155 " Cell Type Proportion")) +
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
156 theme_minimal() +
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
157 ylab(paste0("Proportion of ",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
158 phenotype_scrna_target, " cells")) +
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
159 xlab(paste0("Level of bulk factor (", phenotype_target, ")")) +
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
160 scale_colour_manual(values = c("white", "gray20")) +
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
161 scale_shape_manual(values = c(21, 24)), maxyscale)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
162 }
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
163
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
164 ## BoxPlot
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
165 plot_box <- scale_yaxes(Boxplot_Est(
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
166 sieve_data("list",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
167 data.matrix(estimated_music_props),
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
168 data.matrix(estimated_nnls_props)),
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
169 method.name = methods) +
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
170 theme(axis.text.x = element_text(angle = -90),
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
171 axis.text.y = element_text(size = 8)) +
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
172 ggtitle(element_blank()) + theme_minimal(), maxyscale)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
173
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
174 ## Heatmap
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
175 plot_hmap <- Prop_heat_Est(
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
176 sieve_data(
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
177 "list",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
178 data.matrix(estimated_music_props),
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
179 data.matrix(estimated_nnls_props)),
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
180 method.name = methods) +
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
181 theme(axis.text.x = element_text(angle = -90),
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
182 axis.text.y = element_text(size = 6))
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
183
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
184 pdf(file = outfile_pdf, width = 8, height = 8)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
185 if (length(celltypes) <= 8) {
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
186 plot_grid(jitter_fig, plot_box, labels = "auto", ncol = 1, nrow = 2)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
187 } else {
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
188 print(jitter_fig)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
189 plot_box
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
190 }
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
191 if (use_disease_factor) {
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
192 plot_grid(jitter_new, jitt_compare, labels = "auto", ncol = 1, nrow = 2)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
193 }
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
194 plot_hmap
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
195 message(dev.off())
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
196
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
197 writable <- function(obj, prefix, title) {
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
198 write.table(obj,
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
199 file = paste0("report_data/", prefix, "_",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
200 title, ".tabular"),
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
201 quote = F, sep = "\t", col.names = NA)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
202 }
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
203
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
204 ## Output Proportions
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
205 if ("NNLS" %in% methods) {
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
206 writable(est_prop$Est.prop.allgene, "prop",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
207 "NNLS Estimated Proportions of Cell Types")
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
208 }
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
209
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
210 if ("MuSiC" %in% methods) {
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
211 writable(est_prop$Est.prop.weighted, "prop",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
212 "Music Estimated Proportions of Cell Types")
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
213 writable(est_prop$Weight.gene, "weightgene",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
214 "Music Estimated Proportions of Cell Types (by Gene)")
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
215 writable(est_prop$r.squared.full, "rsquared",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
216 "Music R-sqr Estimated Proportions of Each Subject")
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
217 writable(est_prop$Var.prop, "varprop",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
218 "Matrix of Variance of MuSiC Estimates")
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
219 }
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
220
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
221
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
222 <<<<<<< HEAD
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
223 =======
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
224 write.table(est_prop$Est.prop.weighted,
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
225 file = paste0("report_data/prop_",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
226 "Music Estimated Proportions of Cell Types",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
227 ".tabular"),
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
228 quote = F, sep = "\t", col.names = NA)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
229 write.table(est_prop$Est.prop.allgene,
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
230 file = paste0("report_data/prop_",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
231 "NNLS Estimated Proportions of Cell Types",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
232 ".tabular"),
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
233 quote = F, sep = "\t", col.names = NA)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
234 write.table(est_prop$Weight.gene,
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
235 file = paste0("report_data/weightgene_",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
236 "Music Estimated Proportions of Cell Types (by Gene)",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
237 ".tabular"),
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
238 quote = F, sep = "\t", col.names = NA)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
239 write.table(est_prop$r.squared.full,
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
240 file = paste0("report_data/rsquared_",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
241 "Music R-sqr Estimated Proportions of Each Subject",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
242 ".tabular"),
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
243 quote = F, sep = "\t", col.names = NA)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
244 write.table(est_prop$Var.prop,
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
245 file = paste0("report_data/varprop_",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
246 "Matrix of Variance of MuSiC Estimates",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
247 ".tabular"),
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
248 quote = F, sep = "\t", col.names = NA)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
249
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
250
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
251 >>>>>>> 7a416140 (fitting summaries only apply when disease factor is used)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
252 if (use_disease_factor) {
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
253 ## Summary table of linear regressions of disease factors
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
254 for (meth in methods) {
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
255 ##lm_beta_meth = lm(ct.prop ~ age + bmi + hba1c + gender, data =
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
256 sub_data <- subset(m_prop_ana, Method == meth)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
257
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
258 ## We can only do regression where there are more than 1 factors
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
259 ## so we must find and exclude the ones which are not
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
260 gt1_facts <- sapply(phenotype_factors, function(facname) {
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
261 return(length(unique(sort(sub_data[[facname]]))) == 1)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
262 })
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
263 form_factors <- phenotype_factors
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
264 exclude_facts <- names(gt1_facts)[gt1_facts]
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
265 if (length(exclude_facts) > 0) {
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
266 message("Factors with only one level will be excluded:")
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
267 message(exclude_facts)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
268 form_factors <- phenotype_factors[
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
269 !(phenotype_factors %in% exclude_facts)]
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
270 }
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
271 lm_beta_meth <- lm(as.formula(
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
272 paste("ct.prop", paste(form_factors, collapse = " + "),
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
273 sep = " ~ ")), data = sub_data)
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
274 message(paste0("Summary: ", meth))
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
275 capture.output(summary(lm_beta_meth),
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
276 file = paste0("report_data/summ_Log of ",
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
277 meth,
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
278 " fitting.txt"))
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
279 }
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
280 }
2b58bd6becc7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/music/ commit d5c7ca22af1d4f0eaa7a607886554bebb95e8c50
bgruening
parents:
diff changeset
281