Home » date » 2009 » Dec » 03 »

week 9 Q2 Part 2

*The author of this computation has been verified*
R Software Module: Ian.Holliday/rwasp_One Factor ANOVA.wasp (opens new window with default values)
Title produced by software: Chi Square Measure of Association- Free Statistics Software (Calculator)
Date of computation: Thu, 03 Dec 2009 11:38:09 -0700
 
Cite this page as follows:
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2009/Dec/03/t1259869868cuvknshmtr1nydk.htm/, Retrieved Thu, 03 Dec 2009 20:51:10 +0100
 
BibTeX entries for LaTeX users:
@Manual{KEY,
    author = {{YOUR NAME}},
    publisher = {Office for Research Development and Education},
    title = {Statistical Computations at FreeStatistics.org, URL http://www.freestatistics.org/blog/date/2009/Dec/03/t1259869868cuvknshmtr1nydk.htm/},
    year = {2009},
}
@Manual{R,
    title = {R: A Language and Environment for Statistical Computing},
    author = {{R Development Core Team}},
    organization = {R Foundation for Statistical Computing},
    address = {Vienna, Austria},
    year = {2009},
    note = {{ISBN} 3-900051-07-0},
    url = {http://www.R-project.org},
}
 
Original text written by user:
 
IsPrivate?
No (this computation is public)
 
User-defined keywords:
 
Dataseries X:
» Textbox « » Textfile « » CSV «
36 1 36 1 56 2 48 2 32 2 44 1 39 2 34 2 41 3 50 3 39 1 62 3 52 2 37 3 50 2 41 1 55 2 41 2 56 3 39 2 52 1 46 2 44 2 48 2 41 2 50 3 50 3 44 2 52 1 54 2 44 2 52 3 37 2 52 3 50 3 36 1 50 1 52 3 55 3 31 2 36 1 49 1 42 1 37 2 41 2 30 1 52 1 30 3 41 2 44 1 66 2 48 3 43 2 57 2 46 1 54 3 48 3 48 2 52 1 62 1 58 3 58 2 62 2 48 2 46 2 34 1 66 2 52 3 55 2 55 1 57 3 56 1 55 2 56 3 54 1 55 3 46 2 52 1 32 2 44 1 46 2 59 2 46 3 46 3 54 3 66 3 56 2 59 2 57 2 52 3 48 1 44 1 41 2 50 1 48 3 48 2 59 2 34 46 2 54 2 55 2 54 3 59 2 44 2 54 3 52 3 66 3 44 2 57 2 39 1 60 3 45 2 41 2 50 2 39 2 43 2 48 1 37 2 58 2 46 1 43 1 44 2 34 3 30 1 50 3 39 1 37 2 55 2 48 3 39 1 36 3 43 2 50 3 55 2 43 2 60 3 48 2 30 3 43 2 39 1 52 2 39 1 39 1 56 1 59 1 46 2 57 2 50 2 54 1 50 3 60 3 59 3 41 2 48 1 59 2 60 3 56 2 56 2 51 1
 
Output produced by software:

Enter (or paste) a matrix (table) containing all data (time) series. Every column represents a different variable and must be delimited by a space or Tab. Every row represents a period in time (or category) and must be delimited by hard returns. The easiest way to enter data is to copy and paste a block of spreadsheet cells. Please, do not use commas or spaces to seperate groups of digits!


Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time8 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135


ANOVA Model
childsviq30 ~ Maternalwarmth
means45.7782.1525.222-43.278-43.778-44.778-44.778-43.635-43.278-43.578-43.778-42.778-33.444-43.528-44.178-43.778-43.778-43.778-43.444-43.778-43.778-43.778-43.778-43.778-42.778


ANOVA Statistics
DfSum SqMean SqF valuePr(>F)
Maternalwarmth2479518.23313.25861.2770
Residuals1347245.39854.07


Tukey Honest Significant Difference Comparisons
difflwruprp adj
2-12.152-4.5948.8991
3-15.222-2.25512.70.611
30-1-43.278-63.412-23.1440
34-1-43.778-71.757-15.7990
36-1-44.778-64.912-24.6440
37-1-44.778-64.912-24.6440
39-1-43.635-55.288-31.9820
41-1-43.278-63.412-23.1440
43-1-43.578-56.954-30.2010
44-1-43.778-60.498-27.0570
45-1-42.778-70.757-14.7990
46-1-33.444-50.165-16.7240
48-1-43.528-58.248-28.8080
50-1-44.178-57.554-30.8010
51-1-43.778-71.757-15.7990
52-1-43.778-63.912-23.6440
54-1-43.778-58.498-29.0580
55-1-43.444-60.165-26.7240
56-1-43.778-60.498-27.0570
57-1-43.778-63.912-23.6440
58-1-43.778-71.757-15.7990
59-1-43.778-58.498-29.0580
60-1-43.778-58.498-29.0580
66-1-42.778-70.757-14.7990
3-23.07-3.6779.8160.992
30-2-45.43-65.304-25.5560
34-2-45.93-73.723-18.1380
36-2-46.93-66.804-27.0560
37-2-46.93-66.804-27.0560
39-2-45.787-56.985-34.590
41-2-45.43-65.304-25.5560
43-2-45.73-58.712-32.7480
44-2-45.93-62.337-29.5240
45-2-44.93-72.723-17.1380
46-2-35.597-52.003-19.190
48-2-45.68-60.042-31.3180
50-2-46.33-59.312-33.3480
51-2-45.93-73.723-18.1380
52-2-45.93-65.804-26.0560
54-2-45.93-60.292-31.5680
55-2-45.597-62.003-29.190
56-2-45.93-62.337-29.5240
57-2-45.93-65.804-26.0560
58-2-45.93-73.723-18.1380
59-2-45.93-60.292-31.5680
60-2-45.93-60.292-31.5680
66-2-44.93-72.723-17.1380
30-3-48.5-68.634-28.3660
34-3-49-76.979-21.0210
36-3-50-70.134-29.8660
37-3-50-70.134-29.8660
39-3-48.857-60.51-37.2040
41-3-48.5-68.634-28.3660
43-3-48.8-62.176-35.4240
44-3-49-65.721-32.2790
45-3-48-75.979-20.0210
46-3-38.667-55.387-21.9460
48-3-48.75-63.47-34.030
50-3-49.4-62.776-36.0240
51-3-49-76.979-21.0210
52-3-49-69.134-28.8660
54-3-49-63.72-34.280
55-3-48.667-65.387-31.9460
56-3-49-65.721-32.2790
57-3-49-69.134-28.8660
58-3-49-76.979-21.0210
59-3-49-63.72-34.280
60-3-49-63.72-34.280
66-3-48-75.979-20.0210
34-30-0.5-34.1533.151
36-30-1.5-28.97525.9751
37-30-1.5-28.97525.9751
39-30-0.357-22.38621.6721
41-300-27.47527.4751
43-30-0.3-23.28722.6871
44-30-0.5-25.58124.5811
45-300.5-33.1534.151
46-309.833-15.24834.9140.999
48-30-0.25-24.04423.5441
50-30-0.9-23.88722.0871
51-30-0.5-34.1533.151
52-30-0.5-27.97526.9751
54-30-0.5-24.29423.2941
55-30-0.167-25.24824.9141
56-30-0.5-25.58124.5811
57-30-0.5-27.97526.9751
58-30-0.5-34.1533.151
59-30-0.5-24.29423.2941
60-30-0.5-24.29423.2941
66-300.5-33.1534.151
36-34-1-34.6532.651
37-34-1-34.6532.651
39-340.143-29.22929.5151
41-340.5-33.1534.151
43-340.2-29.89730.2971
44-340-31.72531.7251
45-341-37.85539.8551
46-3410.333-21.39242.0581
48-340.25-30.46830.9681
50-34-0.4-30.49729.6971
51-340-38.85538.8551
52-340-33.6533.651
54-340-30.71830.7181
55-340.333-31.39232.0581
56-340-31.72531.7251
57-340-33.6533.651
58-340-38.85538.8551
59-340-30.71830.7181
60-340-30.71830.7181
66-341-37.85539.8551
37-360-27.47527.4751
39-361.143-20.88623.1721
41-361.5-25.97528.9751
43-361.2-21.78724.1871
44-361-24.08126.0811
45-362-31.6535.651
46-3611.333-13.74836.4140.993
48-361.25-22.54425.0441
50-360.6-22.38723.5871
51-361-32.6534.651
52-361-26.47528.4751
54-361-22.79424.7941
55-361.333-23.74826.4141
56-361-24.08126.0811
57-361-26.47528.4751
58-361-32.6534.651
59-361-22.79424.7941
60-361-22.79424.7941
66-362-31.6535.651
39-371.143-20.88623.1721
41-371.5-25.97528.9751
43-371.2-21.78724.1871
44-371-24.08126.0811
45-372-31.6535.651
46-3711.333-13.74836.4140.993
48-371.25-22.54425.0441
50-370.6-22.38723.5871
51-371-32.6534.651
52-371-26.47528.4751
54-371-22.79424.7941
55-371.333-23.74826.4141
56-371-24.08126.0811
57-371-26.47528.4751
58-371-32.6534.651
59-371-22.79424.7941
60-371-22.79424.7941
66-372-31.6535.651
41-390.357-21.67222.3861
43-390.057-16.0316.1451
44-39-0.143-19.10218.8171
45-390.857-28.51530.2291
46-3910.19-8.76929.150.946
48-390.107-17.11417.3281
50-39-0.543-16.6315.5451
51-39-0.143-29.51529.2291
52-39-0.143-22.17221.8861
54-39-0.143-17.36417.0781
55-390.19-18.76919.151
56-39-0.143-19.10218.8171
57-39-0.143-22.17221.8861
58-39-0.143-29.51529.2291
59-39-0.143-17.36417.0781
60-39-0.143-17.36417.0781
66-390.857-28.51530.2291
43-41-0.3-23.28722.6871
44-41-0.5-25.58124.5811
45-410.5-33.1534.151
46-419.833-15.24834.9140.999
48-41-0.25-24.04423.5441
50-41-0.9-23.88722.0871
51-41-0.5-34.1533.151
52-41-0.5-27.97526.9751
54-41-0.5-24.29423.2941
55-41-0.167-25.24824.9141
56-41-0.5-25.58124.5811
57-41-0.5-27.97526.9751
58-41-0.5-34.1533.151
59-41-0.5-24.29423.2941
60-41-0.5-24.29423.2941
66-410.5-33.1534.151
44-43-0.2-20.26519.8651
45-430.8-29.29730.8971
46-4310.133-9.93130.1980.972
48-430.05-18.38118.4811
50-43-0.6-17.97716.7771
51-43-0.2-30.29729.8971
52-43-0.2-23.18722.7871
54-43-0.2-18.63118.2311
55-430.133-19.93120.1981
56-43-0.2-20.26519.8651
57-43-0.2-23.18722.7871
58-43-0.2-30.29729.8971
59-43-0.2-18.63118.2311
60-43-0.2-18.63118.2311
66-430.8-29.29730.8971
45-441-30.72532.7251
46-4410.333-12.132.7660.991
48-440.25-20.73421.2341
50-44-0.4-20.46519.6651
51-440-31.72531.7251
52-440-25.08125.0811
54-440-20.98420.9841
55-440.333-22.122.7661
56-440-22.43322.4331
57-440-25.08125.0811
58-440-31.72531.7251
59-440-20.98420.9841
60-440-20.98420.9841
66-441-30.72532.7251
46-459.333-22.39241.0581
48-45-0.75-31.46829.9681
50-45-1.4-31.49728.6971
51-45-1-39.85537.8551
52-45-1-34.6532.651
54-45-1-31.71829.7181
55-45-0.667-32.39231.0581
56-45-1-32.72530.7251
57-45-1-34.6532.651
58-45-1-39.85537.8551
59-45-1-31.71829.7181
60-45-1-31.71829.7181
66-450-38.85538.8551
48-46-10.083-31.06810.9010.984
50-46-10.733-30.7989.3310.949
51-46-10.333-42.05821.3921
52-46-10.333-35.41414.7480.998
54-46-10.333-31.31810.6510.979
55-46-10-32.43312.4330.994
56-46-10.333-32.76612.10.991
57-46-10.333-35.41414.7480.998
58-46-10.333-42.05821.3921
59-46-10.333-31.31810.6510.979
60-46-10.333-31.31810.6510.979
66-46-9.333-41.05822.3921
50-48-0.65-19.08117.7811
51-48-0.25-30.96830.4681
52-48-0.25-24.04423.5441
54-48-0.25-19.67819.1781
55-480.083-20.90121.0681
56-48-0.25-21.23420.7341
57-48-0.25-24.04423.5441
58-48-0.25-30.96830.4681
59-48-0.25-19.67819.1781
60-48-0.25-19.67819.1781
66-480.75-29.96831.4681
51-500.4-29.69730.4971
52-500.4-22.58723.3871
54-500.4-18.03118.8311
55-500.733-19.33120.7981
56-500.4-19.66520.4651
57-500.4-22.58723.3871
58-500.4-29.69730.4971
59-500.4-18.03118.8311
60-500.4-18.03118.8311
66-501.4-28.69731.4971
52-510-33.6533.651
54-510-30.71830.7181
55-510.333-31.39232.0581
56-510-31.72531.7251
57-510-33.6533.651
58-510-38.85538.8551
59-510-30.71830.7181
60-510-30.71830.7181
66-511-37.85539.8551
54-520-23.79423.7941
55-520.333-24.74825.4141
56-520-25.08125.0811
57-520-27.47527.4751
58-520-33.6533.651
59-520-23.79423.7941
60-520-23.79423.7941
66-521-32.6534.651
55-540.333-20.65121.3181
56-540-20.98420.9841
57-540-23.79423.7941
58-540-30.71830.7181
59-540-19.42819.4281
60-540-19.42819.4281
66-541-29.71831.7181
56-55-0.333-22.76622.11
57-55-0.333-25.41424.7481
58-55-0.333-32.05831.3921
59-55-0.333-21.31820.6511
60-55-0.333-21.31820.6511
66-550.667-31.05832.3921
57-560-25.08125.0811
58-560-31.72531.7251
59-560-20.98420.9841
60-560-20.98420.9841
66-561-30.72532.7251
58-570-33.6533.651
59-570-23.79423.7941
60-570-23.79423.7941
66-571-32.6534.651
59-580-30.71830.7181
60-580-30.71830.7181
66-581-37.85539.8551
60-590-19.42819.4281
66-591-29.71831.7181
66-601-29.71831.7181


Levenes Test for Homogeneity of Variance
DfF valuePr(>F)
Group243.0720
134
 
Charts produced by software:
http://www.freestatistics.org/blog/date/2009/Dec/03/t1259869868cuvknshmtr1nydk/3ld6g1259865480.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/03/t1259869868cuvknshmtr1nydk/3ld6g1259865480.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Dec/03/t1259869868cuvknshmtr1nydk/41jdi1259865480.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/03/t1259869868cuvknshmtr1nydk/41jdi1259865480.ps (open in new window)


 
Parameters (Session):
par1 = 1 ; par2 = 2 ; par3 = TRUE ;
 
Parameters (R input):
par1 = 1 ; par2 = 2 ; par3 = TRUE ;
 
R code (references can be found in the software module):
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){
thsd<-TukeyHSD(aov.xdf)
bitmap(file='TukeyHSDPlot.png')
plot(thsd)
dev.off()
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<-levene.test(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')
 





Copyright

Creative Commons License

This work is licensed under a Creative Commons Attribution-Noncommercial-Share Alike 3.0 License.

Software written by Ed van Stee & Patrick Wessa


Disclaimer

Information provided on this web site is provided "AS IS" without warranty of any kind, either express or implied, including, without limitation, warranties of merchantability, fitness for a particular purpose, and noninfringement. We use reasonable efforts to include accurate and timely information and periodically update the information, and software without notice. However, we make no warranties or representations as to the accuracy or completeness of such information (or software), and we assume no liability or responsibility for errors or omissions in the content of this web site, or any software bugs in online applications. Your use of this web site is AT YOUR OWN RISK. Under no circumstances and under no legal theory shall we be liable to you or any other person for any direct, indirect, special, incidental, exemplary, or consequential damages arising from your access to, or use of, this web site.


Privacy Policy

We may request personal information to be submitted to our servers in order to be able to:

  • personalize online software applications according to your needs
  • enforce strict security rules with respect to the data that you upload (e.g. statistical data)
  • manage user sessions of online applications
  • alert you about important changes or upgrades in resources or applications

We NEVER allow other companies to directly offer registered users information about their products and services. Banner references and hyperlinks of third parties NEVER contain any personal data of the visitor.

We do NOT sell, nor transmit by any means, personal information, nor statistical data series uploaded by you to third parties.

We carefully protect your data from loss, misuse, alteration, and destruction. However, at any time, and under any circumstance you are solely responsible for managing your passwords, and keeping them secret.

We store a unique ANONYMOUS USER ID in the form of a small 'Cookie' on your computer. This allows us to track your progress when using this website which is necessary to create state-dependent features. The cookie is used for NO OTHER PURPOSE. At any time you may opt to disallow cookies from this website - this will not affect other features of this website.

We examine cookies that are used by third-parties (banner and online ads) very closely: abuse from third-parties automatically results in termination of the advertising contract without refund. We have very good reason to believe that the cookies that are produced by third parties (banner ads) do NOT cause any privacy or security risk.

FreeStatistics.org is safe. There is no need to download any software to use the applications and services contained in this website. Hence, your system's security is not compromised by their use, and your personal data - other than data you submit in the account application form, and the user-agent information that is transmitted by your browser - is never transmitted to our servers.

As a general rule, we do not log on-line behavior of individuals (other than normal logging of webserver 'hits'). However, in cases of abuse, hacking, unauthorized access, Denial of Service attacks, illegal copying, hotlinking, non-compliance with international webstandards (such as robots.txt), or any other harmful behavior, our system engineers are empowered to log, track, identify, publish, and ban misbehaving individuals - even if this leads to ban entire blocks of IP addresses, or disclosing user's identity.


FreeStatistics.org is powered by