Statistics R command

. Read a file
height <- read.table("H:/public folder/hzhang/R/height.txt",header=T) . Print summary summary(height) . Set the factor level height$Treatment <- relevel(height$Treatment, ref="kitaake") . Variance analysis summary(lm((height$Height~height$Treatment)))