2018 Hakyung R command

From Crop Genomics Lab.
Jump to: navigation, search

library(ggplot2) ggplot(summary, aes(x=summary$no, y = summary$m)) + geom_point() + geom_errorbar(aes(ymin=summary$m-summary$se, ymax=summarWarning message: width=.1)

> pdf("RWC_distribution.pdf") > hist(xvalues, xlim=c(0,1), xlab="Relative Water Contents", ylab = "Density", col="lightgray", main= "RWC distribution after> curve(dnorm(x, mean= mean(xvalues), sd = sd(xvalues)), add = TRUE, col = "red", lwd = 2) dev.off()


ggplot(summary, aes(x=summary$no, y = summary$m)) + geom_point() + geom_errorbar(aes(ymin=summary$m-summary$se, ymax=summary$m+summary$se)

data.frame()

ref1 <- unlist(ref1, use.names = FALSE)

shapiro.test(ref1)

m <- read.table("./tem/transpiration_initial_transpiration",header = T) res_transpiration_initial = lmm(initial_transpiration~genotype + 1|rep, data = m, method = "reml") res_transpiration_initial


m <- read.table("initial_transpiration_fixed_sort") pdf("RWC_distribution.pdf") hist(m$V2, xlab="fitted initial transpiraiton", ylab = "Density", col="lightgray", main= "initial transpiraion", freq = FALSE) curve(dnorm(x, mean= mean(m$V2), sd = sd(m$V2)), add = TRUE, col = "red", lwd = 2) dev.off()