Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
2012 4.35 2012 12.7 2012 18.1 2012 17.85 2012 16.6 2012 12.6 2012 17.1 2012 19.1 2012 16.1 2012 13.35 2012 18.4 2012 14.7 2012 10.6 2012 12.6 2012 16.2 2012 13.6 2012 18.9 2012 14.1 2012 14.5 2012 16.15 2012 14.75 2012 14.8 2012 12.45 2012 12.65 2012 17.35 2012 8.6 2012 18.4 2012 16.1 2012 11.6 2012 17.75 2012 15.25 2012 17.65 2012 15.6 2012 16.35 2012 17.65 2012 13.6 2012 11.7 2012 14.35 2012 14.75 2012 18.25 2012 9.9 2012 16 2012 18.25 2012 16.85 2012 14.6 2012 13.85 2012 18.95 2012 15.6 2012 14.85 2012 11.75 2012 18.45 2012 15.9 2012 17.1 2012 16.1 2012 19.9 2012 10.95 2012 18.45 2012 15.1 2012 15 2012 11.35 2012 15.95 2012 18.1 2012 14.6 2012 15.4 2012 15.4 2012 17.6 2012 13.35 2012 19.1 2012 15.35 2012 7.6 2012 13.4 2012 13.9 2012 19.1 2012 15.25 2012 12.9 2012 16.1 2012 17.35 2012 13.15 2012 12.15 2012 12.6 2012 10.35 2012 15.4 2012 9.6 2012 18.2 2012 13.6 2012 14.85 2012 14.75 2012 14.1 2012 14.9 2012 16.25 2012 19.25 2012 13.6 2012 13.6 2012 15.65 2012 12.75 2012 14.6 2012 9.85 2012 12.65 2012 11.9 2012 19.2 2012 16.6 2012 11.2 2012 15.25 2012 11.9 2012 13.2 2012 16.35 2012 12.4 2012 15.85 2012 14.35 2012 18.15 2012 11.15 2012 15.65 2012 17.75 2012 7.65 2012 12.35 2012 15.6 2012 19.3 2012 15.2 2012 17.1 2012 15.6 2012 18.4 2012 19.05 2012 18.55 2012 19.1 2012 13.1 2012 12.85 2012 9.5 2012 4.5 2012 11.85 2012 13.6 2012 11.7 2012 12.4 2012 13.35 2012 11.4 2012 14.9 2012 19.9 2012 17.75 2012 11.2 2012 14.6 2012 17.6 2012 14.05 2012 16.1 2012 13.35 2012 11.85 2012 11.95 2012 14.75 2012 15.15 2012 13.2 2012 16.85 2012 7.85 2012 7.7 2012 12.6 2012 7.85 2012 10.95 2012 12.35 2012 9.95 2012 14.9 2012 16.65 2012 13.4 2012 13.95 2012 15.7 2012 16.85 2012 10.95 2012 15.35 2012 12.2 2012 15.1 2012 17.75 2012 15.2 2012 14.6 2012 16.65 2012 8.1 2011 12.9 2011 7.4 2011 12.2 2011 12.8 2011 7.4 2011 6.7 2011 12.6 2011 14.8 2011 13.3 2011 11.1 2011 8.2 2011 11.4 2011 6.4 2011 10.6 2011 12 2011 6.3 2011 11.3 2011 11.9 2011 9.3 2011 9.6 2011 10 2011 6.4 2011 13.8 2011 10.8 2011 13.8 2011 11.7 2011 10.9 2011 16.1 2011 13.4 2011 9.9 2011 11.5 2011 8.3 2011 11.7 2011 6.1 2011 9 2011 9.7 2011 10.8 2011 10.3 2011 10.4 2011 12.7 2011 9.3 2011 11.8 2011 5.9 2011 11.4 2011 13 2011 10.8 2011 12.3 2011 11.3 2011 11.8 2011 7.9 2011 12.7 2011 12.3 2011 11.6 2011 6.7 2011 10.9 2011 12.1 2011 13.3 2011 10.1 2011 5.7 2011 14.3 2011 8 2011 13.3 2011 9.3 2011 12.5 2011 7.6 2011 15.9 2011 9.2 2011 9.1 2011 11.1 2011 13 2011 14.5 2011 12.2 2011 12.3 2011 11.4 2011 8.8 2011 14.6 2011 7.3 2011 12.6 2011 NA 2011 13 2011 12.6 2011 13.2 2011 9.9 2011 7.7 2011 10.5 2011 13.4 2011 10.9 2011 4.3 2011 10.3 2011 11.8 2011 11.2 2011 11.4 2011 8.6 2011 13.2 2011 12.6 2011 5.6 2011 9.9 2011 8.8 2011 7.7 2011 9 2011 7.3 2011 11.4 2011 13.6 2011 7.9 2011 10.7 2011 10.3 2011 8.3 2011 9.6 2011 14.2 2011 8.5 2011 13.5 2011 4.9 2011 6.4 2011 9.6 2011 11.6 2011 11.1
Names of X columns:
Jaar Totale_ptn
Response Variable (column number)
Factor Variable (column number)
Include Intercept Term ?
TRUE
TRUE
FALSE
Chart options
Title:
Label y-axis:
Label x-axis:
R Code
cat1 <- as.numeric(par1) # cat2<- as.numeric(par2) # intercept<-as.logical(par3) x <- t(x) x1<-as.numeric(x[,cat1]) f1<-as.character(x[,cat2]) xdf<-data.frame(x1,f1) (V1<-dimnames(y)[[1]][cat1]) (V2<-dimnames(y)[[1]][cat2]) names(xdf)<-c('Response', 'Treatment') if(intercept == FALSE) (lmxdf<-lm(Response ~ Treatment - 1, data = xdf) ) else (lmxdf<-lm(Response ~ Treatment, data = xdf) ) (aov.xdf<-aov(lmxdf) ) (anova.xdf<-anova(lmxdf) ) load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'ANOVA Model', length(lmxdf$coefficients)+1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, paste(V1, ' ~ ', V2), length(lmxdf$coefficients)+1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, 'means',,TRUE) for(i in 1:length(lmxdf$coefficients)){ a<-table.element(a, round(lmxdf$coefficients[i], digits=3),,FALSE) } a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab') a<-table.start() a<-table.row.start(a) a<-table.element(a,'ANOVA Statistics', 5+1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, ' ',,TRUE) a<-table.element(a, 'Df',,FALSE) a<-table.element(a, 'Sum Sq',,FALSE) a<-table.element(a, 'Mean Sq',,FALSE) a<-table.element(a, 'F value',,FALSE) a<-table.element(a, 'Pr(>F)',,FALSE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, V2,,TRUE) a<-table.element(a, anova.xdf$Df[1],,FALSE) a<-table.element(a, round(anova.xdf$'Sum Sq'[1], digits=3),,FALSE) a<-table.element(a, round(anova.xdf$'Mean Sq'[1], digits=3),,FALSE) a<-table.element(a, round(anova.xdf$'F value'[1], digits=3),,FALSE) a<-table.element(a, round(anova.xdf$'Pr(>F)'[1], digits=3),,FALSE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, 'Residuals',,TRUE) a<-table.element(a, anova.xdf$Df[2],,FALSE) a<-table.element(a, round(anova.xdf$'Sum Sq'[2], digits=3),,FALSE) a<-table.element(a, round(anova.xdf$'Mean Sq'[2], digits=3),,FALSE) a<-table.element(a, ' ',,FALSE) a<-table.element(a, ' ',,FALSE) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable1.tab') bitmap(file='anovaplot.png') boxplot(Response ~ Treatment, data=xdf, xlab=V2, ylab=V1) dev.off() if(intercept==TRUE){ 'Tukey Plot' thsd<-TukeyHSD(aov.xdf) bitmap(file='TukeyHSDPlot.png') plot(thsd) dev.off() } if(intercept==TRUE){ a<-table.start() a<-table.row.start(a) a<-table.element(a,'Tukey Honest Significant Difference Comparisons', 5,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, ' ', 1, TRUE) for(i in 1:4){ a<-table.element(a,colnames(thsd[[1]])[i], 1, TRUE) } a<-table.row.end(a) for(i in 1:length(rownames(thsd[[1]]))){ a<-table.row.start(a) a<-table.element(a,rownames(thsd[[1]])[i], 1, TRUE) for(j in 1:4){ a<-table.element(a,round(thsd[[1]][i,j], digits=3), 1, FALSE) } a<-table.row.end(a) } a<-table.end(a) table.save(a,file='mytable2.tab') } if(intercept==FALSE){ a<-table.start() a<-table.row.start(a) a<-table.element(a,'TukeyHSD Message', 1,TRUE) a<-table.row.end(a) a<-table.start() a<-table.row.start(a) a<-table.element(a,'Must Include Intercept to use Tukey Test ', 1, FALSE) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable2.tab') } library(car) lt.lmxdf<-leveneTest(lmxdf) a<-table.start() a<-table.row.start(a) a<-table.element(a,'Levenes Test for Homogeneity of Variance', 4,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,' ', 1, TRUE) for (i in 1:3){ a<-table.element(a,names(lt.lmxdf)[i], 1, FALSE) } a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Group', 1, TRUE) for (i in 1:3){ a<-table.element(a,round(lt.lmxdf[[i]][1], digits=3), 1, FALSE) } a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,' ', 1, TRUE) a<-table.element(a,lt.lmxdf[[1]][2], 1, FALSE) a<-table.element(a,' ', 1, FALSE) a<-table.element(a,' ', 1, FALSE) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable3.tab')
Compute
Summary of computational transaction
Raw Input
view raw input (R code)
Raw Output
view raw output of R engine
Computing time
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation