Difference between revisions of "2018 Hakyung R command"
(Created page with "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("...") |
|||
| Line 14: | Line 14: | ||
shapiro.test(ref1) | 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() | ||
Latest revision as of 05:55, 7 February 2020
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()
