Blog & Share Statistical Computations at FreeStatistics.org

Home » date » 2010 » Jun » 03 »

*Unverified author*

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

Title produced by software: Standard Deviation-Mean Plot

Date of computation: Thu, 03 Jun 2010 09:01:40 +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/Jun/03/t1275555745bhmwjhn5b3e2ifx.htm/, Retrieved Thu, 03 Jun 2010 11:02:28 +0200

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/Jun/03/t1275555745bhmwjhn5b3e2ifx.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:

KDGP2W81

Dataseries X:

» Textbox « » Textfile « » CSV «

20 25 15 15 25 25 25 21 30 25 20 40 13 30 25 20 25 20 25 20 20 15 15 12 20 5 20 15 25 22 20 22 25 20 20 35 30 25 20 20 20 25 25 15 20 35 25 25 30 23 10 22 25 25 22 30 20 25 25 22 25 25 25 22 25 12 18 20 20 22 30 25 22 20 50 30 25 20 30 22 25 30 22 25 22 22 25 25 25 20 22 15 20 30 20 25 30 35 22 12 30 15 10 30 9 25 20 20 35 25 35 30 12 25 15 25 25 20 20 6 15 40 20 40 25 25 20 15 15 22 24 22 20 25 25 25 35 40 20 22 22 20 25 25 18 25 20 25 30 20 22 35 22 25 25 25 25 22 23 35 15 25 18 22 25 25 28 30 20 25 25 30 22 30 10 10 25 20 22 25 25 15 22 25 25 28 22 30 25 20 25 25 20 30 20 30 50 19 20 28 20 25 35 25 25 15 16 20 20 25 30 20 25 25 25 20 20 25 25 30 22 20 25 25 18 18 20 25 25 30 25 20 25 20 20 20 22 18 22 20 15 25 25 20 25 15 22 25 25 15 12 25 30 22 15 22 25 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	2 seconds
R Server	'Gwilym Jenkins' @ 72.249.127.135

Standard Deviation-Mean Plot
Section	Mean	Standard Deviation	Range
1	23.8333333333333	6.73975091707716	25
2	20	5.54321862524587	18
3	20.75	6.91671230360994	30
4	23.75	5.27644853011086	20
5	23.25	5.1720402163943	20
6	22.4166666666667	4.58174903747731	18
7	26.75	8.20338067783366	30
8	22.5833333333333	3.77692355166894	15
9	21.5	8.72301032275608	26
10	22.75	8.76070773396762	29
11	23.5833333333333	8.50089121887532	25
12	25.3333333333333	6.15457454896664	20
13	24.3333333333333	4.63844278989744	17
14	24.4166666666667	5.23030215246315	20
15	22	6.46669790682863	20
16	23.9166666666667	3.84845502255205	15
17	26.8333333333333	8.94257775866092	31
18	22.1666666666667	4.36584540044721	15
19	22.75	3.64629326103298	12
20	22.25	3.38781238076620	12
21	20.75	5.1012476192327	13
22	24	7.38548945875996	28
23	23.75	4.33012701892219	15
24	28	11.1273618534755	31
25	21.9166666666667	2.77843426585856	7
26	24.0833333333333	10.7910597114913	42
27	23.4166666666667	3.60450055381106	12
28	20.3333333333333	3.77391365063169	15
29	28.9166666666667	10.7657565991209	35
30	23.6666666666667	7.41415517572279	30
31	16.9166666666667	5.51787728179684	21
32	21.5833333333333	5.33357007050340	18
33	21.5	9.27851869055126	30
34	21.1666666666667	6.04277680500873	20
35	22.3333333333333	6.77338821261471	25
36	26.6666666666667	8.34847109936722	30
37	20.5833333333333	4.62126179804397	13
38	23.6666666666667	10.4736321157930	38
39	19.0833333333333	5.23030215246315	18
40	22.9166666666667	3.44985726538290	12
41	16.5833333333333	5.66421515598881	20
42	20.8333333333333	4.260459500044	10
43	28.25	8.87411967464942	35
44	22.3333333333333	4.4585634321817	13
45	26	13.5579564027246	50
46	19.5	5.35129551009507	15
47	27.9166666666667	7.21687836487032	25
48	25.5833333333333	3.08834563931546	10
49	24.4166666666667	5.28218848125615	20
50	25.3333333333333	4.92365963917331	15
51	26.1666666666667	5.23681605789648	20
52	19.75	4.02548698352702	10
53	20.5	8.55463723262525	35
54	24.9166666666667	2.90636709604442	10
55	23.9166666666667	2.93747985163948	10
56	24.9166666666667	3.70401093001856	10
57	21.5	4.23191338629359	15
58	23.4166666666667	5.46822778891984	20
59	23.25	3.76889180722205	15
60	23.75	4.9012985292397	17
61	21	5.15222811465627	20
62	21.1666666666667	4.68718433280546	15
63	22.1666666666667	7.7440103931811	25
64	22.1666666666667	3.45972498469126	12
65	21.9166666666667	5.36755121639338	18
66	22.5833333333333	5.28218848125615	20
67	23.25	6.78400525299968	25
68	29.1666666666667	11.2478954260409	30
69	25.8333333333333	6.02771377334171	22
70	22.25	9.91899005122818	28
71	19	5.18739731552258	14
72	17.9166666666667	5.23030215246315	16
73	22.9166666666667	9.80221063502174	35
74	26.8333333333333	11.2882346903421	32
75	30	8.79049072991532	30

Regression: S.E.(k) = alpha + beta * Mean(k)
alpha	-2.67284598886476
beta	0.383249416624467
S.D.	0.095612159510457
T-STAT	4.00837527974201
p-value	0.000145896888059319

Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha	-1.69292497882997
beta	1.09874908689156
S.D.	0.354809686550163
T-STAT	3.09672798838941
p-value	0.00277588110166215
Lambda	-0.098749086891564

Charts produced by software:

http://www.freestatistics.org/blog/date/2010/Jun/03/t1275555745bhmwjhn5b3e2ifx/1enda1275555698.png (open in new window)

http://www.freestatistics.org/blog/date/2010/Jun/03/t1275555745bhmwjhn5b3e2ifx/1enda1275555698.ps (open in new window)

http://www.freestatistics.org/blog/date/2010/Jun/03/t1275555745bhmwjhn5b3e2ifx/2enda1275555698.png (open in new window)

http://www.freestatistics.org/blog/date/2010/Jun/03/t1275555745bhmwjhn5b3e2ifx/2enda1275555698.ps (open in new window)

Parameters (Session):

par1 = 12 ;

Parameters (R input):

par1 = 12 ;

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

par1 <- as.numeric(par1)

(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