R version 2.13.0 (2011-04-13)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> x <- array(list(3.7
+ ,6.6
+ ,0.8564
+ ,3
+ ,6.5
+ ,0.8973
+ ,2.7
+ ,6.3
+ ,0.9383
+ ,2.5
+ ,6.2
+ ,0.9217
+ ,2.2
+ ,6.3
+ ,0.9095
+ ,2.9
+ ,6.5
+ ,0.892
+ ,3.1
+ ,6.6
+ ,0.8742
+ ,3
+ ,6.4
+ ,0.8532
+ ,2.8
+ ,6.2
+ ,0.8607
+ ,2.5
+ ,6.3
+ ,0.9005
+ ,1.9
+ ,6.6
+ ,0.9111
+ ,1.9
+ ,7.1
+ ,0.9059
+ ,1.8
+ ,7.2
+ ,0.8883
+ ,2
+ ,7.3
+ ,0.8924
+ ,2.6
+ ,7.3
+ ,0.8833
+ ,2.5
+ ,7.3
+ ,0.87
+ ,2.5
+ ,7.4
+ ,0.8758
+ ,1.6
+ ,7.3
+ ,0.8858
+ ,1.4
+ ,7.4
+ ,0.917
+ ,0.8
+ ,7.4
+ ,0.9554
+ ,1.1
+ ,7.6
+ ,0.9922
+ ,1.3
+ ,7.6
+ ,0.9778
+ ,1.2
+ ,7.7
+ ,0.9808
+ ,1.3
+ ,7.7
+ ,0.9811
+ ,1.1
+ ,7.8
+ ,1.0014
+ ,1.3
+ ,7.8
+ ,1.0183
+ ,1.2
+ ,8
+ ,1.0622
+ ,1.6
+ ,8.1
+ ,1.0773
+ ,1.7
+ ,8.1
+ ,1.0807
+ ,1.5
+ ,8.2
+ ,1.0848
+ ,0.9
+ ,8.1
+ ,1.1582
+ ,1.5
+ ,8.1
+ ,1.1663
+ ,1.4
+ ,8.1
+ ,1.1372
+ ,1.6
+ ,8
+ ,1.1139
+ ,1.7
+ ,8.2
+ ,1.1222
+ ,1.4
+ ,8.4
+ ,1.1692
+ ,1.8
+ ,8.4
+ ,1.1702
+ ,1.7
+ ,8.5
+ ,1.2286
+ ,1.4
+ ,8.6
+ ,1.2613
+ ,1.2
+ ,8.5
+ ,1.2646
+ ,1
+ ,8.3
+ ,1.2262
+ ,1.7
+ ,7.8
+ ,1.1985
+ ,2.4
+ ,7.8
+ ,1.2007
+ ,2
+ ,8
+ ,1.2138
+ ,2.1
+ ,8.6
+ ,1.2266
+ ,2
+ ,8.9
+ ,1.2176
+ ,1.8
+ ,8.9
+ ,1.2218
+ ,2.7
+ ,8.6
+ ,1.249
+ ,2.3
+ ,8.3
+ ,1.2991
+ ,1.9
+ ,8.3
+ ,1.3408
+ ,2
+ ,8.3
+ ,1.3119
+ ,2.3
+ ,8.4
+ ,1.3014
+ ,2.8
+ ,8.5
+ ,1.3201
+ ,2.4
+ ,8.4
+ ,1.2938
+ ,2.3
+ ,8.6
+ ,1.2694
+ ,2.7
+ ,8.5
+ ,1.2165
+ ,2.7
+ ,8.5
+ ,1.2037
+ ,2.9
+ ,8.5
+ ,1.2292
+ ,3
+ ,8.5
+ ,1.2256
+ ,2.2
+ ,8.5
+ ,1.2015
+ ,2.3
+ ,8.5
+ ,1.1786
+ ,2.8
+ ,8.5
+ ,1.1856
+ ,2.8
+ ,8.5
+ ,1.2103
+ ,2.8
+ ,8.6
+ ,1.1938
+ ,2.2
+ ,8.6
+ ,1.202
+ ,2.6
+ ,8.6
+ ,1.2271
+ ,2.8
+ ,8.6
+ ,1.277
+ ,2.5
+ ,8.4
+ ,1.265
+ ,2.4
+ ,8
+ ,1.2684
+ ,2.3
+ ,7.9
+ ,1.2811
+ ,1.9
+ ,8
+ ,1.2727
+ ,1.7
+ ,8
+ ,1.2611
+ ,2
+ ,8
+ ,1.2881
+ ,2.1
+ ,8
+ ,1.3213
+ ,1.7
+ ,7.9
+ ,1.2999
+ ,1.8
+ ,7.9
+ ,1.3074
+ ,1.8
+ ,7.9
+ ,1.3242
+ ,1.8
+ ,8
+ ,1.3516
+ ,1.3
+ ,7.9
+ ,1.3511
+ ,1.3
+ ,7.5
+ ,1.3419
+ ,1.3
+ ,7.2
+ ,1.3716
+ ,1.2
+ ,7
+ ,1.3622
+ ,1.4
+ ,6.9
+ ,1.3896
+ ,2.2
+ ,7.1
+ ,1.4227
+ ,2.9
+ ,7.1
+ ,1.4684
+ ,3.1
+ ,7.2
+ ,1.457
+ ,3.5
+ ,7.1
+ ,1.4718
+ ,3.6
+ ,6.9
+ ,1.4748
+ ,4.4
+ ,6.8
+ ,1.5527
+ ,4.1
+ ,6.7
+ ,1.575
+ ,5.1
+ ,6.7
+ ,1.5557
+ ,5.8
+ ,6.9
+ ,1.5553
+ ,5.9
+ ,7.3
+ ,1.577
+ ,5.4
+ ,7.4
+ ,1.4975
+ ,5.5
+ ,7.3
+ ,1.437
+ ,4.8
+ ,7.1
+ ,1.3322
+ ,3.2
+ ,7
+ ,1.2732
+ ,2.7
+ ,7.1
+ ,1.3449
+ ,2.1
+ ,7.5
+ ,1.3239
+ ,1.9
+ ,7.7
+ ,1.2785
+ ,0.6
+ ,7.8
+ ,1.305
+ ,0.7
+ ,7.7
+ ,1.319
+ ,-0.2
+ ,7.7
+ ,1.365
+ ,-1
+ ,7.8
+ ,1.4016
+ ,-1.7
+ ,8
+ ,1.4088
+ ,-0.7
+ ,8.1
+ ,1.4268
+ ,-1
+ ,8.1
+ ,1.4562
+ ,-0.9
+ ,8
+ ,1.4816
+ ,0
+ ,8.1
+ ,1.4914
+ ,0.3
+ ,8.2
+ ,1.4614
+ ,0.8
+ ,8.3
+ ,1.4272
+ ,0.8
+ ,8.4
+ ,1.3686
+ ,1.9
+ ,8.5
+ ,1.3569
+ ,2.1
+ ,8.5
+ ,1.3406
+ ,2.5
+ ,8.5
+ ,1.2565
+ ,2.7
+ ,8.5
+ ,1.2208
+ ,2.4
+ ,8.5
+ ,1.277
+ ,2.4
+ ,8.3
+ ,1.2894
+ ,2.9
+ ,8.2
+ ,1.3067
+ ,3.1
+ ,8.1
+ ,1.3898
+ ,3
+ ,7.9
+ ,1.3661
+ ,3.4
+ ,7.6
+ ,1.322
+ ,3.7
+ ,7.3
+ ,1.336
+ ,3.5
+ ,7.1
+ ,1.3649
+ ,3.5
+ ,7
+ ,1.3999
+ ,3.3
+ ,7
+ ,1.4442
+ ,3.1
+ ,7
+ ,1.4349
+ ,3.4
+ ,7
+ ,1.4388
+ ,4
+ ,6.9
+ ,1.4264
+ ,3.4
+ ,6.8
+ ,1.4343
+ ,3.4
+ ,6.7
+ ,1.377
+ ,3.4
+ ,6.6
+ ,1.3706)
+ ,dim=c(3
+ ,132)
+ ,dimnames=list(c('Inflatie'
+ ,'Werkloosheid'
+ ,'Wisselkoers')
+ ,1:132))
>  y <- array(NA,dim=c(3,132),dimnames=list(c('Inflatie','Werkloosheid','Wisselkoers'),1:132))
>  for (i in 1:dim(x)[1])
+  {
+  	for (j in 1:dim(x)[2])
+  	{
+  		y[i,j] <- as.numeric(x[i,j])
+  	}
+  }
> par3 = 'TRUE'
> par2 = '2'
> par1 = '1'
> ylab = 'Y Variable Name'
> xlab = 'X Variable Name'
> main = 'Title Goes Here'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Dr. Ian E. Holliday
> #To cite this work: Ian E. Holliday, 2009, YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: 
> #Technical description: 
> 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])
[1] "Inflatie"
> (V2<-dimnames(y)[[1]][cat2])
[1] "Werkloosheid"
> names(xdf)<-c('Response', 'Treatment')
> if(intercept == FALSE) (lmxdf<-lm(Response ~ Treatment - 1, data = xdf) ) else (lmxdf<-lm(Response ~ Treatment, data = xdf) )

Call:
lm(formula = Response ~ Treatment, data = xdf)

Coefficients:
 (Intercept)  Treatment6.3  Treatment6.4  Treatment6.5  Treatment6.6  
      2.6500       -0.1833        0.3500        0.3000        0.3750  
Treatment6.7  Treatment6.8  Treatment6.9    Treatment7  Treatment7.1  
      1.5500        1.2500        1.0500        0.3000        0.4214  
Treatment7.2  Treatment7.3  Treatment7.4  Treatment7.5  Treatment7.6  
     -0.5833        0.7500       -0.1250       -0.9500       -0.7167  
Treatment7.7  Treatment7.8  Treatment7.9    Treatment8  Treatment8.1  
     -1.6700       -1.6333       -0.6667       -1.3682       -1.7056  
Treatment8.2  Treatment8.3  Treatment8.4  Treatment8.5  Treatment8.6  
     -1.0500       -0.9167       -0.7833       -0.2312       -0.2875  
Treatment8.9  
     -0.7500  

> (aov.xdf<-aov(lmxdf) )
Call:
   aov(formula = lmxdf)

Terms:
                Treatment Residuals
Sum of Squares   95.97532 113.62733
Deg. of Freedom        25       106

Residual standard error: 1.035353 
Estimated effects may be unbalanced
> (anova.xdf<-anova(lmxdf) )
Analysis of Variance Table

Response: Response
           Df  Sum Sq Mean Sq F value    Pr(>F)    
Treatment  25  95.975   3.839  3.5813 2.338e-06 ***
Residuals 106 113.627   1.072                      
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
> 
> #Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/wessaorg/rcomp/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="/var/wessaorg/rcomp/tmp/1gvz61324304744.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="/var/wessaorg/rcomp/tmp/2ew881324304744.tab") 
> postscript(file="/var/wessaorg/rcomp/tmp/3h9ve1324304744.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> boxplot(Response ~ Treatment, data=xdf, xlab=V2, ylab=V1)
> dev.off()
null device 
          1 
> if(intercept==TRUE){
+ thsd<-TukeyHSD(aov.xdf)
+ postscript(file="/var/wessaorg/rcomp/tmp/49xbw1324304744.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
+ plot(thsd)
+ dev.off()
+ }
null device 
          1 
> 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="/var/wessaorg/rcomp/tmp/513kq1324304744.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="/var/wessaorg/rcomp/tmp/6fl0z1324304744.tab") 
+ }
> library(car)
Loading required package: MASS
Loading required package: nnet
Loading required package: survival
Loading required package: splines
> lt.lmxdf<-levene.test(lmxdf)
Warning message:
'levene.test' is deprecated.
Use 'leveneTest' instead.
See help("Deprecated") and help("car-deprecated"). 
> 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="/var/wessaorg/rcomp/tmp/7vld11324304744.tab") 
> 
> try(system("convert tmp/3h9ve1324304744.ps tmp/3h9ve1324304744.png",intern=TRUE))
character(0)
> try(system("convert tmp/49xbw1324304744.ps tmp/49xbw1324304744.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  3.190   0.164   3.348