Home » date » 2009 » Oct » 21 »

WS 3 - Vraag 1 (3)

*The author of this computation has been verified*

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

Title produced by software: Percentiles

Date of computation: Wed, 21 Oct 2009 09:42:06 -0600

Cite this page as follows:

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2009/Oct/21/t1256139781elij451hnj808l7.htm/, Retrieved Sun, 26 May 2013 09:20:43 +0000

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords:

SHW WS 3 - Vraag 1 (3)

System-generated keywords (parent):

t1255509090iv1ijspnf5sojst (pk = 46388)

Estimated Impact

Dataseries X:

» Textfile « » CSV « » Stem and Leaf « » Histogram « » Kernel Density « » Harrell-Davis Quantiles « » Central Tendency « » Variability «

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	'Sir Ronald Aylmer Fisher' @ 193.190.124.24

Percentiles - Ungrouped Data
p	Weighted Average at Xnp	Weighted Average at X(n+1)p	Empirical Distribution Function	Empirical Distribution Function - Averaging	Empirical Distribution Function - Interpolation	Closest Observation	True Basic - Statistics Graphics Toolkit	MS Excel (old versions)
0.01	46.266	46.2666	46.29	46.29	46.3608	46.29	46.2534	46.29
0.02	46.776	46.7776	46.84	46.84	46.912	46.76	46.8224	46.76
0.03	47.652	47.6547	47.67	47.67	47.7239	47.67	47.5953	47.67
0.04	48.15	48.16	48.3	48.3	48.354	48.05	48.19	48.05
0.05	48.84	48.843	48.84	48.87	48.897	48.84	48.897	48.84
0.06	49.202	49.2062	49.23	49.23	49.3218	49.23	49.1838	49.23
0.07	49.968	49.9708	50	50	50	49.96	49.9892	49.96
0.08	50.498	50.5028	50.51	50.51	50.7836	50.51	50.4572	50.51
0.09	52.482	52.5522	52.95	52.95	52.9531	52.17	52.5678	52.17
0.1	53.31	53.311	53.31	53.315	53.319	53.31	53.319	53.31
0.11	53.718	53.7378	53.79	53.79	53.8341	53.79	53.6622	53.79
0.12	54.824	54.8624	55.08	55.08	55.092	54.76	54.9776	54.76
0.13	55.604	55.6144	55.62	55.62	55.7004	55.62	55.5456	55.62
0.14	55.874	55.8824	55.91	55.91	55.9828	55.85	55.8776	55.91
0.15	57.11	57.209	57.11	57.44	57.671	57.11	57.671	57.11
0.16	58.34	58.372	58.42	58.42	58.4772	58.42	58.268	58.42
0.17	58.668	58.6918	58.78	58.78	58.7854	58.64	58.7282	58.64
0.18	59.152	59.1682	59.17	59.17	59.883	59.17	59.0818	59.17
0.19	60.772	60.8252	60.94	60.94	61.0345	60.66	60.7748	60.94
0.2	61.85	61.854	61.85	61.86	61.866	61.85	61.866	61.85
0.21	62.472	62.5077	62.54	62.54	62.7935	62.54	62.4023	62.54
0.22	64.7	64.7	64.7	64.7	64.7	64.7	64.7	64.7
0.23	65.282	65.3171	65.3	65.3	65.6249	65.3	65.8529	65.3
0.24	66.118	66.1828	66.28	66.28	66.2848	66.01	66.1072	66.28
0.25	66.64	66.7175	66.64	66.795	66.8725	66.64	66.8725	66.64
0.26	67.322	67.3532	67.37	67.37	67.3938	67.37	67.2668	67.37
0.27	67.62	67.62	67.62	67.62	67.62	67.62	67.62	67.62
0.28	68.3	68.3808	68.38	68.38	68.3852	68.38	68.3892	68.38
0.29	68.496	68.5221	68.55	68.55	68.5544	68.46	68.4879	68.55
0.3	68.77	68.797	68.77	68.815	68.833	68.77	68.833	68.77
0.31	69.426	69.4601	69.47	69.47	69.5077	69.47	69.3699	69.47
0.32	69.784	69.7904	69.8	69.8	69.7976	69.78	69.7896	69.8
0.33	70	70.0165	70.01	70.01	70.0335	70.01	70.0535	70.01
0.34	70.362	70.3722	70.38	70.38	70.3854	70.35	70.3578	70.38
0.35	70.91	70.945	70.91	70.91	70.975	70.91	70.975	70.91
0.36	71.454	71.5224	71.53	71.53	71.5924	71.53	71.3476	71.53
0.37	71.982	72.0227	72.07	72.07	72.0513	71.96	72.0073	72.07
0.38	72.392	72.486	72.45	72.45	72.534	72.45	72.614	72.45
0.39	72.87	72.9675	73.02	73.02	73.0206	72.77	72.8225	73.02
0.4	73.29	73.378	73.29	73.4	73.422	73.29	73.422	73.29
0.41	74.234	74.3909	74.39	74.39	74.4071	74.39	74.4791	74.39
0.42	74.586	74.6196	74.65	74.65	74.6324	74.57	74.6004	74.65
0.43	74.9	74.9146	74.91	74.91	74.9174	74.91	74.9254	74.91
0.44	75.208	75.3928	75.46	75.46	75.4432	75.04	75.1072	75.46
0.45	75.69	75.7305	75.69	75.735	75.7395	75.69	75.7395	75.69
0.46	75.85	75.876	75.87	75.87	75.884	75.87	75.964	75.87
0.47	76.086	76.1236	76.15	76.15	76.1284	76.07	76.0964	76.15
0.48	76.198	76.2196	76.2	76.2	76.2224	76.2	76.2504	76.2
0.49	76.422	76.4612	76.47	76.47	76.4628	76.39	76.3988	76.47
0.5	76.65	76.775	76.65	76.775	76.775	76.65	76.775	76.775
0.51	77.108	77.1411	77.14	77.14	77.1409	77.14	77.1489	77.14
0.52	77.42	77.42	77.42	77.42	77.42	77.42	77.42	77.42
0.53	77.726	77.7833	77.78	77.78	77.7827	77.78	77.7867	77.78
0.54	78.172	78.1882	78.19	78.19	78.1858	78.16	78.1618	78.19
0.55	78.38	78.5395	78.67	78.67	78.5105	78.38	78.5105	78.67
0.56	78.836	78.868	78.86	78.86	78.862	78.86	78.902	78.86
0.57	79.246	79.2916	79.31	79.31	79.2804	79.23	79.2484	79.31
0.58	79.888	79.8976	79.89	79.89	79.8944	79.89	79.9024	79.89
0.59	80.012	80.0297	80.03	80.03	80.0243	80	80.0003	80.03
0.6	80.47	80.476	80.47	80.475	80.474	80.47	80.474	80.48
0.61	80.676	80.7	80.7	80.7	80.6994	80.7	80.7	80.7
0.62	81.15	81.181	81.19	81.19	81.169	81.14	81.149	81.19
0.63	81.578	81.6273	81.58	81.58	81.5987	81.58	81.6427	81.58
0.64	81.938	81.9572	81.95	81.95	81.9452	81.93	82.1228	81.95
0.65	82.39	82.4225	82.39	82.415	82.4075	82.39	82.4075	82.44
0.66	82.974	83.1116	83.07	83.07	83.0556	83.07	83.1884	83.07
0.67	83.51	83.6105	83.63	83.63	83.5595	83.48	83.4995	83.63
0.68	84.012	84.0636	84.03	84.03	84.0384	84.03	84.0664	84.03
0.69	84.38	84.4454	84.44	84.44	84.411	84.34	84.4946	84.44
0.7	84.67	84.901	84.67	84.67	84.769	84.67	84.769	85
0.71	85.214	85.2362	85.23	85.23	85.2256	85.23	85.2438	85.23
0.72	85.374	85.5684	85.59	85.59	85.4496	85.32	85.3416	85.59
0.73	86.196	86.3248	86.24	86.24	86.2512	86.24	86.3152	86.4
0.74	86.766	86.834	86.82	86.82	86.7894	86.73	86.906	86.82
0.75	86.98	87.0475	86.98	87.025	87.0025	86.98	87.0025	87.07
0.76	87.28	87.3052	87.28	87.28	87.28	87.28	87.3248	87.28
0.77	87.502	87.5867	87.59	87.59	87.5273	87.48	87.4833	87.59
0.78	87.894	87.9796	87.91	87.91	87.9124	87.91	87.9604	88.03
0.79	88.374	88.41	88.41	88.41	88.3866	88.35	88.41	88.41
0.8	89.27	89.334	89.27	89.31	89.286	89.27	89.286	89.35
0.81	89.55	89.7181	89.55	89.55	89.55	89.55	89.7919	89.55
0.82	90.292	90.3028	90.3	90.3	90.2938	90.29	90.4372	90.3
0.83	90.642	90.7089	90.69	90.69	90.6828	90.69	90.7011	90.72
0.84	91.51	91.5172	91.51	91.51	91.51	91.51	91.5328	91.51
0.85	91.74	92.9045	91.74	92.425	91.9455	91.74	91.9455	93.11
0.86	93.284	93.4964	93.34	93.34	93.3036	93.34	93.5236	93.34
0.87	94.404	94.4221	94.42	94.42	94.4066	94.4	94.4479	94.42
0.88	95.208	95.422	95.32	95.32	95.2752	95.32	95.368	95.47
0.89	96.958	97.4943	97.3	97.3	97.0207	96.73	97.7757	97.3
0.9	99.2	99.254	99.2	99.23	99.206	99.2	99.206	99.26
0.91	100.064	101.0317	100.18	100.18	100.0901	100.18	100.9983	101.85
0.92	102.264	102.5752	102.52	102.52	102.2896	102.2	102.9248	102.52
0.93	104.674	104.9709	104.73	104.73	104.6936	104.73	104.8191	105.06
0.94	106.368	106.9936	106.98	106.98	106.4292	105.96	107.0064	106.98
0.95	107.93	109.507	107.93	108.76	108.013	107.93	108.013	109.59
0.96	110.954	111.6844	111.27	111.27	110.9856	111.27	111.5956	112.01
0.97	113.842	115.8206	115.45	115.45	113.9023	113.44	117.2594	115.45
0.98	120.514	123.8866	121.18	121.18	120.5806	121.18	121.9434	124.65
0.99	138.578	151.2605	147.38	147.38	138.7247	132.71	153.4495	147.38

Charts produced by software:

http://www.freestatistics.org/blog/date/2009/Oct/21/t1256139781elij451hnj808l7/1zngx1256139721.png (opens in new window)

http://www.freestatistics.org/blog/date/2009/Oct/21/t1256139781elij451hnj808l7/1zngx1256139721.ps (opens in new window)

Click here to open pdf file.

Parameters (Session):

Parameters (R input):

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

x <-sort(x[!is.na(x)])
q1 <- function(data,n,p,i,f) {
np <- n*p;
i <<- floor(np)
f <<- np - i
qvalue <- (1-f)*data[i] + f*data[i+1]
}
q2 <- function(data,n,p,i,f) {
np <- (n+1)*p
i <<- floor(np)
f <<- np - i
qvalue <- (1-f)*data[i] + f*data[i+1]
}
q3 <- function(data,n,p,i,f) {
np <- n*p
i <<- floor(np)
f <<- np - i
if (f==0) {
qvalue <- data[i]
} else {
qvalue <- data[i+1]
}
}
q4 <- function(data,n,p,i,f) {
np <- n*p
i <<- floor(np)
f <<- np - i
if (f==0) {
qvalue <- (data[i]+data[i+1])/2
} else {
qvalue <- data[i+1]
}
}
q5 <- function(data,n,p,i,f) {
np <- (n-1)*p
i <<- floor(np)
f <<- np - i
if (f==0) {
qvalue <- data[i+1]
} else {
qvalue <- data[i+1] + f*(data[i+2]-data[i+1])
}
}
q6 <- function(data,n,p,i,f) {
np <- n*p+0.5
i <<- floor(np)
f <<- np - i
qvalue <- data[i]
}
q7 <- function(data,n,p,i,f) {
np <- (n+1)*p
i <<- floor(np)
f <<- np - i
if (f==0) {
qvalue <- data[i]
} else {
qvalue <- f*data[i] + (1-f)*data[i+1]
}
}
q8 <- function(data,n,p,i,f) {
np <- (n+1)*p
i <<- floor(np)
f <<- np - i
if (f==0) {
qvalue <- data[i]
} else {
if (f == 0.5) {
qvalue <- (data[i]+data[i+1])/2
} else {
if (f < 0.5) {
qvalue <- data[i]
} else {
qvalue <- data[i+1]
}
}
}
}
lx <- length(x)
qval <- array(NA,dim=c(99,8))
mystep <- 25
mystart <- 25
if (lx>10){
mystep=10
mystart=10
}
if (lx>20){
mystep=5
mystart=5
}
if (lx>50){
mystep=2
mystart=2
}
if (lx>=100){
mystep=1
mystart=1
}
for (perc in seq(mystart,99,mystep)) {
qval[perc,1] <- q1(x,lx,perc/100,i,f)
qval[perc,2] <- q2(x,lx,perc/100,i,f)
qval[perc,3] <- q3(x,lx,perc/100,i,f)
qval[perc,4] <- q4(x,lx,perc/100,i,f)
qval[perc,5] <- q5(x,lx,perc/100,i,f)
qval[perc,6] <- q6(x,lx,perc/100,i,f)
qval[perc,7] <- q7(x,lx,perc/100,i,f)
qval[perc,8] <- q8(x,lx,perc/100,i,f)
}
bitmap(file='test1.png')
myqqnorm <- qqnorm(x,col=2)
qqline(x)
grid()
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Percentiles - Ungrouped Data',9,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'p',1,TRUE)
a<-table.element(a,hyperlink('http://www.xycoon.com/method_1.htm', 'Weighted Average at Xnp',''),1,TRUE)
a<-table.element(a,hyperlink('http://www.xycoon.com/method_2.htm','Weighted Average at X(n+1)p',''),1,TRUE)
a<-table.element(a,hyperlink('http://www.xycoon.com/method_3.htm','Empirical Distribution Function',''),1,TRUE)
a<-table.element(a,hyperlink('http://www.xycoon.com/method_4.htm','Empirical Distribution Function - Averaging',''),1,TRUE)
a<-table.element(a,hyperlink('http://www.xycoon.com/method_5.htm','Empirical Distribution Function - Interpolation',''),1,TRUE)
a<-table.element(a,hyperlink('http://www.xycoon.com/method_6.htm','Closest Observation',''),1,TRUE)
a<-table.element(a,hyperlink('http://www.xycoon.com/method_7.htm','True Basic - Statistics Graphics Toolkit',''),1,TRUE)
a<-table.element(a,hyperlink('http://www.xycoon.com/method_8.htm','MS Excel (old versions)',''),1,TRUE)
a<-table.row.end(a)
for (perc in seq(mystart,99,mystep)) {
a<-table.row.start(a)
a<-table.element(a,round(perc/100,2),1,TRUE)
for (j in 1:8) {
a<-table.element(a,round(qval[perc,j],6))
}
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable.tab')