Home » date » 2010 » Dec » 17 »

workshop9

*The author of this computation has been verified*

R Software Module: /rwasp_smp.wasp (opens new window with default values)

Title produced by software: Standard Deviation-Mean Plot

Date of computation: Fri, 17 Dec 2010 20:13:38 +0000

Cite this page as follows:

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2010/Dec/17/t1292616732fyfawn08kvhdbig.htm/, Retrieved Fri, 17 Dec 2010 21:12:14 +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/2010/Dec/17/t1292616732fyfawn08kvhdbig.htm/},
    year = {2010},
}
@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 = {2010},
    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 «

40.7819 39.5915 38.8859 39.9068 41.47 41.5613 41.6005 41.4113 41.84 42.2892 43.1521 43.5998 43.116 42.4185 42.3687 42.2975 42.8528 43.535 44.7265 45.7293 45.7585 46.1685 46.5075 46.527 46.601 46.4607 46.7135 46.4113 45.55 44.6081 44.4395 44.9847 45.7558 45.3942 45.697 45.5664 46.0205 45.9195 45.8005 45.535 45.4977 45.5782 45.7697 45.2445 45.0615 45.2865 44.791 44.7625 44.7644 44.9973 44.7265 45.1465 44.7465 45.1795 45.6515 45.492 45.2775 45.2115 45.411 45.4005 44.7692 44.8913 45.032 44.879 44.833 44.8257 44.7815 44.479 44.6317 44.5043 44.3217 44.1005 44.047 43.6835 43.7864 44.1807 43.9595 43.937 43.991 43.865 43.671 43.93 43.863 43.7095 43.9435 43.736 43.6295 43.598 43.8726 43.8935 43.5957 43.7155 43.528 43.3415 43.3374 43.332 43.3869 43.5016 43.4875 43.6023 43.3886 43.3105 43.4455 43.5185 43.5755 43.6217 43.644 43.5789 43.5215 43.5033 43.632 43.263 43.3717 43.2745 43.2647 43.324 43.4455 43.409 etc...

Output produced by software:

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	3 seconds
R Server	'George Udny Yule' @ 72.249.76.132

Standard Deviation-Mean Plot
Section	Mean	Standard Deviation	Range
1	40.12722	1.01342491927128	2.5841
2	41.74046	0.343168039012961	0.877900000000004
3	42.93102	0.526602598740265	1.23110000000000
4	43.82822	1.39656166960145	3.4318
5	46.3125	0.351752576394261	0.842500000000001
6	45.94872	0.86850904543361	2.10540000000000
7	45.25424	0.548552807849892	1.31630000000000
8	45.76838	0.213714522201932	0.485500000000002
9	45.43032	0.279349515839924	0.708199999999998
10	44.92034	0.226911641393736	0.523999999999994
11	45.0901	0.379673939058239	0.924999999999997
12	45.3585	0.112374930478286	0.280499999999996
13	44.8809	0.097100823889398	0.262799999999999
14	44.64444	0.157172319445887	0.346699999999998
15	43.98782	0.255370117672369	0.638199999999998
16	43.98664	0.117978061519929	0.315700000000000
17	43.8234	0.126050287583964	0.272500000000001
18	43.74592	0.135425281982354	0.295500000000004
19	43.50362	0.164224106025880	0.378099999999996
20	43.46206	0.105406939999224	0.270299999999999
21	43.44772	0.104483213962819	0.265000000000001
22	43.57388	0.0611432089442456	0.140699999999995
23	43.36118	0.158046882284972	0.369
24	43.50386	0.242388011254683	0.606000000000002
25	43.78538	0.139616159523173	0.341000000000001
26	43.86982	0.114427037888779	0.274799999999999
27	45.0249	0.70376597672806	1.765
28	44.54592	0.462245072445344	0.955100000000002
29	44.0922	0.226226225270193	0.517499999999998
30	44.5103	0.224176716007707	0.625
31	45.21262	0.379760479776398	0.874400000000001
32	45.4414	0.399843469372706	0.985000000000007
33	45.12938	0.176941436639358	0.450099999999999
34	45.30092	0.249225965741935	0.637
35	45.26668	0.139517192488955	0.291600000000003
36	44.99462	0.155694836137878	0.398500000000006
37	44.62458	0.173039394358627	0.427100000000003
38	44.15918	0.194474193660752	0.525599999999997
39	43.89526	0.0981817600168177	0.250599999999999
40	43.70722	0.111230625279191	0.299999999999997
41	43.73062	0.0675053479362911	0.182599999999994
42	43.4974	0.110854296263159	0.231500000000004
43	43.1261	0.0956401327895363	0.225999999999999
44	43.12716	0.0921451192413367	0.226700000000001
45	42.95376	0.103973135953475	0.280499999999996
46	43.0443	0.113046273711254	0.2577
47	43.79158	0.206421975574306	0.529999999999994
48	44.30164	0.432317531451131	1.04
49	44.42894	0.421978841886651	0.901200000000003
50	43.92952	0.235748227989099	0.5869
51	43.09346	0.349838602787056	0.789999999999999
52	42.9763	0.237336154009456	0.542499999999997
53	42.70778	0.200714578444117	0.539999999999999
54	42.25776	0.303126480202569	0.676199999999994
55	42.52038	0.211694525200817	0.552099999999996
56	42.0367	0.270623262119131	0.655000000000001
57	41.6316	0.160696297406009	0.391000000000005
58	41.07754	0.148847599913468	0.366
59	40.714	0.110038629580705	0.262999999999998
60	40.63842	0.099697577703775	0.216000000000001
61	40.4127	0.157336661334859	0.371500000000005
62	40.18062	0.145367197125072	0.355499999999999
63	39.73466	0.176119641153395	0.412500000000001
64	39.6601	0.108580269846783	0.289500000000004
65	39.23766	0.100055199764928	0.251099999999994
66	39.3759	0.109799476319333	0.234999999999999
67	38.95914	0.183024378157665	0.483499999999999
68	38.79744	0.0323951076553235	0.0812999999999988
69	38.68662	0.184017450259480	0.387000000000000
70	37.89804	0.310675204997117	0.758000000000003
71	38.28758	0.591900652136826	1.44290000000000
72	38.1085	0.286782583153162	0.710000000000001
73	37.95776	0.124670317237105	0.318000000000005
74	38.13634	0.168261486977857	0.423999999999999
75	38.3131	0.136839504529940	0.342499999999994
76	38.17328	0.134517162473790	0.302500000000002
77	38.62236	0.43094049821292	0.950099999999999
78	39.03716	0.158797443304356	0.383800000000001
79	39.14092	0.369592298350493	0.8185
80	39.15734	0.0855640286569063	0.235799999999998
81	39.39024	0.0633032621592299	0.154400000000003
82	39.40022	0.137493607124113	0.338999999999999
83	39.25644	0.156470182462984	0.353999999999999
84	39.0399	0.115549686282571	0.248400000000004
85	39.03162	0.0496847260231951	0.1235
86	39.33416	0.132234443319431	0.335999999999999
87	39.27492	0.0381661106218600	0.102400000000003
88	40.2559	0.436378906456303	1.03050000000000
89	41.24094	0.262605489280025	0.660500000000006
90	41.56386	0.290716912132747	0.654499999999999
91	41.5485	0.108700275988611	0.280000000000001
92	42.18536	0.346158667376682	0.884999999999998
93	42.5818	0.276033231695027	0.712800000000001
94	42.56126	0.267116047814427	0.611799999999995
95	43.1798	0.23092060323843	0.615000000000002
96	42.25398	0.283607179034666	0.779700000000005
97	42.38106	0.0620769119721652	0.145699999999998
98	41.6929	0.418318419388868	1.14650000000000

Regression: S.E.(k) = alpha + beta * Mean(k)
alpha	-0.388060213957644
beta	0.0148742742195947
S.D.	0.0085834395887362
T-STAT	1.73290369971425
p-value	0.0863239046359169

Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha	-12.0961543921023
beta	2.78611752886801
S.D.	1.18221991665320
T-STAT	2.3566829568862
p-value	0.0204704805611829
Lambda	-1.78611752886801

Charts produced by software:

http://www.freestatistics.org/blog/date/2010/Dec/17/t1292616732fyfawn08kvhdbig/16htz1292616815.png (open in new window)

http://www.freestatistics.org/blog/date/2010/Dec/17/t1292616732fyfawn08kvhdbig/16htz1292616815.ps (open in new window)

http://www.freestatistics.org/blog/date/2010/Dec/17/t1292616732fyfawn08kvhdbig/26htz1292616815.png (open in new window)

http://www.freestatistics.org/blog/date/2010/Dec/17/t1292616732fyfawn08kvhdbig/26htz1292616815.ps (open in new window)

Parameters (Session):

par1 = 4 ;

Parameters (R input):

par1 = 4 ;

R code (references can be found in the software module):

par1 <- 5

(n <- length(x))

(np <- floor(n / par1))

arr <- array(NA,dim=c(par1,np))

j <- 0

k <- 1

for (i in 1:(np*par1))

{

j = j + 1

arr[j,k] <- x[i]

if (j == par1) {

j = 0

k=k+1

}

}

arr

arr.mean <- array(NA,dim=np)

arr.sd <- array(NA,dim=np)

arr.range <- array(NA,dim=np)

for (j in 1:np)

{

arr.mean[j] <- mean(arr[,j],na.rm=TRUE)

arr.sd[j] <- sd(arr[,j],na.rm=TRUE)

arr.range[j] <- max(arr[,j],na.rm=TRUE) - min(arr[,j],na.rm=TRUE)

}

arr.mean

arr.sd

arr.range

(lm1 <- lm(arr.sd~arr.mean))

(lnlm1 <- lm(log(arr.sd)~log(arr.mean)))

(lm2 <- lm(arr.range~arr.mean))

bitmap(file='test1.png')

plot(arr.mean,arr.sd,main='Standard Deviation-Mean Plot',xlab='mean',ylab='standard deviation')

dev.off()

bitmap(file='test2.png')

plot(arr.mean,arr.range,main='Range-Mean Plot',xlab='mean',ylab='range')

dev.off()

load(file='createtable')

a<-table.start()

a<-table.row.start(a)

a<-table.element(a,'Standard Deviation-Mean Plot',4,TRUE)

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,'Section',header=TRUE)

a<-table.element(a,'Mean',header=TRUE)

a<-table.element(a,'Standard Deviation',header=TRUE)

a<-table.element(a,'Range',header=TRUE)

a<-table.row.end(a)

for (j in 1:np) {

a<-table.row.start(a)

a<-table.element(a,j,header=TRUE)

a<-table.element(a,arr.mean[j])

a<-table.element(a,arr.sd[j] )

a<-table.element(a,arr.range[j] )

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,'Regression: S.E.(k) = alpha + beta * Mean(k)',2,TRUE)

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,'alpha',header=TRUE)

a<-table.element(a,lm1$coefficients[[1]])

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,'beta',header=TRUE)

a<-table.element(a,lm1$coefficients[[2]])

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,'S.D.',header=TRUE)

a<-table.element(a,summary(lm1)$coefficients[2,2])

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,'T-STAT',header=TRUE)

a<-table.element(a,summary(lm1)$coefficients[2,3])

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,'p-value',header=TRUE)

a<-table.element(a,summary(lm1)$coefficients[2,4])

a<-table.row.end(a)

a<-table.end(a)

table.save(a,file='mytable1.tab')

a<-table.start()

a<-table.row.start(a)

a<-table.element(a,'Regression: ln S.E.(k) = alpha + beta * ln Mean(k)',2,TRUE)

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,'alpha',header=TRUE)

a<-table.element(a,lnlm1$coefficients[[1]])

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,'beta',header=TRUE)

a<-table.element(a,lnlm1$coefficients[[2]])

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,'S.D.',header=TRUE)

a<-table.element(a,summary(lnlm1)$coefficients[2,2])

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,'T-STAT',header=TRUE)

a<-table.element(a,summary(lnlm1)$coefficients[2,3])

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,'p-value',header=TRUE)

a<-table.element(a,summary(lnlm1)$coefficients[2,4])

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,'Lambda',header=TRUE)

a<-table.element(a,1-lnlm1$coefficients[[2]])

a<-table.row.end(a)

a<-table.end(a)

table.save(a,file='mytable2.tab')

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.

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