Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- dat <- readRDS("TCGA GLUT data.rds")
- str(dat)
- a <- round(summary(lm(SLC2A1 ~ Gender, data=dat))$coefficients, 5)
- b <- round(summary(lm(SLC2A1 ~ Race, data=dat))$coefficients, 5)
- c <- round(summary(lm(SLC2A1 ~ Age_Dx, data=dat))$coefficients, 5)
- d <- round(summary(lm(SLC2A1 ~ Recurrence, data=dat))$coefficients, 5)
- e <- round(summary(lm(SLC2A1 ~ Vital_Status, data=dat))$coefficients, 5)
- f <- round(summary(lm(SLC2A1 ~ Hashimoto, data=dat))$coefficients, 5)
- g <- round(summary(lm(SLC2A1 ~ Histologic_Dx, data=dat))$coefficients, 5)
- h <- round(summary(lm(SLC2A1 ~ Max_Size, data=dat))$coefficients, 5)
- i <- round(summary(lm(SLC2A1 ~ Metastatic_LN, data=dat))$coefficients, 5)
- j <- round(summary(lm(SLC2A1 ~ ETE, data=dat))$coefficients, 5)
- k <- round(summary(lm(SLC2A1 ~ T_stage, data=dat))$coefficients, 5)
- l <- round(summary(lm(SLC2A1 ~ N_stage, data=dat))$coefficients, 5)
- m <- round(summary(lm(SLC2A1 ~ Stage, data=dat))$coefficients, 5)
- n <- round(summary(lm(SLC2A1 ~ BRAF_V600E, data=dat))$coefficients, 5)
- SLC2A1.result <- rbind(a,b,c,d,e,f,g,h,i,j,k,l,m,n)
- SLC2A1.result
- result <- data.frame()
- for (i in 28:63){
- a <- summary(lm(dat[,i] ~ Gender, data=dat))$coefficients
- b <- summary(lm(dat[,i] ~ Race, data=dat))$coefficients
- c <- summary(lm(dat[,i] ~ Age_Dx, data=dat))$coefficients
- d <- summary(lm(dat[,i] ~ Recurrence, data=dat))$coefficients
- e <- summary(lm(dat[,i] ~ Vital_Status, data=dat))$coefficients
- f <- summary(lm(dat[,i] ~ Hashimoto, data=dat))$coefficients
- g <- summary(lm(dat[,i] ~ Histologic_Dx, data=dat))$coefficients
- h <- summary(lm(dat[,i] ~ Max_Size, data=dat))$coefficients
- i <- summary(lm(dat[,i] ~ Metastatic_LN, data=dat))$coefficients
- j <- summary(lm(dat[,i] ~ ETE, data=dat))$coefficients
- k <- summary(lm(dat[,i] ~ T_stage, data=dat))$coefficients
- l <- summary(lm(dat[,i] ~ N_stage, data=dat))$coefficients
- m <- summary(lm(dat[,i] ~ Stage, data=dat))$coefficients
- n <- summary(lm(dat[,i] ~ BRAF_V600E, data=dat))$coefficients
- result[i] <- rbind(a,b,c,d,e,f,g,h,i,j,k,l,m,n)
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement