x=1:5*10 x=1:5*10 y=100+10*x+rnorm(length(x),5) plot(x,y) rl=lm(y~x) abline(rl) x=1:5*10 y=100+10*x+rnorm(length(x),10) plot(x,y) rl=lm(y~x) abline(rl) x=1:5*10 y=100+10*x+rnorm(length(x),50) plot(x,y) rl=lm(y~x) abline(rl) x=1:5*10 y=100+10*x+rnorm(length(x),500) plot(x,y) rl=lm(y~x) abline(rl) x=1:5*10 y=100+10*x+rnorm(length(x),0,500) plot(x,y) rl=lm(y~x) abline(rl) x=seq(10,50,length=25) y=100+10*x+rnorm(length(x),0,500) plot(x,y) rl=lm(y~x) abline(rl) cor(x,y) cor.test(x,y) cov(x,y)/var(x) cov(x,y)/sd(x)/sd(y) x=seq(10,50,length=25) y=100+10*x+rnorm(length(x),0,500) plot(x,y) rl=lm(y~x) abline(rl) x=seq(10,50,length=25) y=100+10*x+rnorm(length(x),0,100) plot(x,y) rl=lm(y~x) abline(rl) cor.test(x,y) cor.test(x,y,confidence=0.01) ?cor.test cor.test(x,y,conf.level = =0.01) cor.test(x,y,conf.level = 0.01) cor.test(x,y,conf.level = 0.99) cor.test(x,y,conf.level = 0.99999) cor.test(x,y,conf.level = 0.999999999) cor.test(x,y,conf.level = 0.99) d=read.table("clipboard",h=T) d head(d) table(d) table(d$det) length(d$det) xdd=mean(d$time) xdd SST=sum((d$time-xdd)^2) SST aov(d$time~d$det) summary(aov(d$time~d$det)) 187+6015 46.67/158.29 1-pf(0.2948386,4,38) d=read.table("clipboard",h=T) summary(aov(d$time~d$det)) d=read.table("clipboard",h=T) summary(aov(d$time~d$det)) savehistory("N:/321_su19_data/201906124_in_class.Rhistory")