scan("c://julia//MAE5704//ref.txt",dec="," ) read.table("c://JULIA//MAE5704//tab1.txt",dec=",",col.names=c("metodo","yield")) x <- rnorm(10) x print(x) x y <- 2*x +1 + rnorm(10,sd=0.01) lm(y ~ x) z <- lm(y ~ x) z summary(z) anova(z) plot(z) idades <- scan() idades fix(idades) idades x y c(y,x) w <- c(y,x) w boxexpr <- read.table("c://Julia//MAE5704//tab1.txt",dec=",",header=T) boxexpr boxexpr$metodo boxexpr$yield attach(boxexpr) ls(2) str(boxexpr) metodo yield ?boxplot help.search("linear model") ?lm lm metodo yield mean(yield) tapply(yield,metodo,mean) tapply(yield,metodo,median) tapply(yield,metodo,sd) boxplot(yield ~ metodo) boxplot(yield ~ metodo,main="Boxplot") diff(tapply(yield,metodo,mean)) x x[1] x[10] x[c(1,10)] x[-c(1,10)] tapply(yield,metodo,mean)[1] tapply(yield,metodo,mean)[2] tapply(yield,metodo,mean)[2]-tapply(yield,metodo,mean)[1] diff(tapply(yield,metodo,mean)) dotchart(yield) ordem <- 1:20 ordem dotchart(yield,labels=ordem,groups=metodo) library(lattice) dotplot(yield,groups=metodo) refs <- scan("c://julia//mae5704//ref.txt",dec=",") refs meds <- filter(refs,filter=rep(1/10,10)) meds meds <- as.vector(meds) meds meds[1] meds[5] meds[201] meds[203] meds[208] meds[c(1,203)] meds[-c(1,203)] is.na(meds) meds[is.na(meds)] meds[!is.na(meds)] meds <- meds[!is.na(meds)] meds diff(meds) diff(meds,lag=10) diff(meds,lag=10)[1] difs <- diff(meds,lag=10) round(difs,1) stripchart(round(difs,1),method='stack',pch=1) abline(v=1.3) mtext(1.3,at=1.5) save.image("C:/Program Files/R/rw2001/.RData")