> > rawData<-read.table("/Users/cshallah/Desktop/11g_St.txt",header=FALSE) > attach(rawData) > ds1<-V1 > dends1<-density(ds1) > > > ds<-ds1 > fivenum(ds) [1] 0.004676033 0.004708505 0.004725338 0.004736550 0.004778740 > summary(ds) Min. 1st Qu. Median Mean 3rd Qu. Max. 0.004676 0.004709 0.004725 0.004725 0.004736 0.004779 > sd(ds) [1] 2.190288e-05 > shapiro.test(ds) Shapiro-Wilk normality test data: ds W = 0.9896, p-value = 0.7029 > hist(ds) > plot(density(ds)) > > rawData<-read.table("/Users/cshallah/Desktop/12c_St.txt",header=FALSE) > attach(rawData) The following object(s) are masked from 'rawData (position 3)': V1 > ds2<-V1 > dends2<-density(ds2) > > > ds<-ds2 > fivenum(ds) [1] 0.004359095 0.004393360 0.004407711 0.004437973 0.004947699 > summary(ds) Min. 1st Qu. Median Mean 3rd Qu. Max. 0.004359 0.004393 0.004408 0.004420 0.004438 0.004948 > sd(ds) [1] 6.3159e-05 > shapiro.test(ds) Shapiro-Wilk normality test data: ds W = 0.4578, p-value < 2.2e-16 > hist(ds) > plot(density(ds)) > > t.test(ds1,ds2) Welch Two Sample t-test data: ds1 and ds2 t = 43.3707, df = 110.102, p-value < 2.2e-16 alternative hypothesis: true difference in means is not equal to 0 95 percent confidence interval: 0.0002916475 0.0003195763 sample estimates: mean of x mean of y 0.004725306 0.004419694 > > wilcox.test(ds1,ds2) Wilcoxon rank sum test with continuity correction data: ds1 and ds2 W = 8010, p-value < 2.2e-16 alternative hypothesis: true location shift is not equal to 0 > > > plot(range(dends1$x, dends2$x), range(dends1$y, dends2$y), type = "n", xlab = "Sample Values", ylab = "Probability", main="Probability Density Curve\n11g:red, 12c:blue") > lines(dends1, col = "red") > lines(dends2, col = "blue") >