Home » date » 2010 » Dec » 17 »

Histogram Happiness

*The author of this computation has been verified*

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

Title produced by software: Histogram

Date of computation: Fri, 17 Dec 2010 15:52:25 +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/t1292601039bdp3e6t5blnx2d5.htm/, Retrieved Thu, 23 May 2013 09:09:05 +0000

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords:

System-generated keywords (parent):

(pk = 0)

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	1 seconds
R Server	'George Udny Yule' @ 72.249.76.132

Frequency Table (Histogram)
Bins	Midpoint	Abs. Frequency	Rel. Frequency	Cumul. Rel. Freq.	Density
[7,8[	7.5	1	0.006897	0.006897	0.006897
[8,9[	8.5	1	0.006897	0.013793	0.006897
[9,10[	9.5	5	0.034483	0.048276	0.034483
[10,11[	10.5	4	0.027586	0.075862	0.027586
[11,12[	11.5	11	0.075862	0.151724	0.075862
[12,13[	12.5	12	0.082759	0.234483	0.082759
[13,14[	13.5	17	0.117241	0.351724	0.117241
[14,15[	14.5	24	0.165517	0.517241	0.165517
[15,16[	15.5	29	0.2	0.717241	0.2
[16,17[	16.5	21	0.144828	0.862069	0.144828
[17,18[	17.5	12	0.082759	0.944828	0.082759
[18,19]	18.5	8	0.055172	1	0.055172

Charts produced by software:

http://www.freestatistics.org/blog/date/2010/Dec/17/t1292601039bdp3e6t5blnx2d5/1018a1292601143.png (opens in new window)

http://www.freestatistics.org/blog/date/2010/Dec/17/t1292601039bdp3e6t5blnx2d5/1018a1292601143.ps (opens in new window)

Click here to open pdf file.

Parameters (Session):

par2 = grey ; par3 = FALSE ; par4 = Unknown ;

Parameters (R input):

par2 = grey ; par3 = FALSE ; par4 = Unknown ;

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

par1 <- as.numeric(par1)
if (par3 == 'TRUE') par3 <- TRUE
if (par3 == 'FALSE') par3 <- FALSE
if (par4 == 'Unknown') par1 <- as.numeric(par1)
if (par4 == 'Interval/Ratio') par1 <- as.numeric(par1)
if (par4 == '3-point Likert') par1 <- c(1:3 - 0.5, 3.5)
if (par4 == '4-point Likert') par1 <- c(1:4 - 0.5, 4.5)
if (par4 == '5-point Likert') par1 <- c(1:5 - 0.5, 5.5)
if (par4 == '6-point Likert') par1 <- c(1:6 - 0.5, 6.5)
if (par4 == '7-point Likert') par1 <- c(1:7 - 0.5, 7.5)
if (par4 == '8-point Likert') par1 <- c(1:8 - 0.5, 8.5)
if (par4 == '9-point Likert') par1 <- c(1:9 - 0.5, 9.5)
if (par4 == '10-point Likert') par1 <- c(1:10 - 0.5, 10.5)
bitmap(file='test1.png')
if(is.numeric(x[1])) {
if (is.na(par1)) {
myhist<-hist(x,col=par2,main=main,xlab=xlab,right=par3)
} else {
if (par1 < 0) par1 <- 3
if (par1 > 50) par1 <- 50
myhist<-hist(x,breaks=par1,col=par2,main=main,xlab=xlab,right=par3)
}
} else {
plot(mytab <- table(x),col=par2,main='Frequency Plot',xlab=xlab,ylab='Absolute Frequency')
}
dev.off()
if(is.numeric(x[1])) {
myhist
n <- length(x)
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,hyperlink('http://www.xycoon.com/histogram.htm','Frequency Table (Histogram)',''),6,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Bins',header=TRUE)
a<-table.element(a,'Midpoint',header=TRUE)
a<-table.element(a,'Abs. Frequency',header=TRUE)
a<-table.element(a,'Rel. Frequency',header=TRUE)
a<-table.element(a,'Cumul. Rel. Freq.',header=TRUE)
a<-table.element(a,'Density',header=TRUE)
a<-table.row.end(a)
crf <- 0
if (par3 == FALSE) mybracket <- '[' else mybracket <- ']'
mynumrows <- (length(myhist$breaks)-1)
for (i in 1:mynumrows) {
a<-table.row.start(a)
if (i == 1)
dum <- paste('[',myhist$breaks[i],sep='')
else
dum <- paste(mybracket,myhist$breaks[i],sep='')
dum <- paste(dum,myhist$breaks[i+1],sep=',')
if (i==mynumrows)
dum <- paste(dum,']',sep='')
else
dum <- paste(dum,mybracket,sep='')
a<-table.element(a,dum,header=TRUE)
a<-table.element(a,myhist$mids[i])
a<-table.element(a,myhist$counts[i])
rf <- myhist$counts[i]/n
crf <- crf + rf
a<-table.element(a,round(rf,6))
a<-table.element(a,round(crf,6))
a<-table.element(a,round(myhist$density[i],6))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable.tab')
} else {
mytab
reltab <- mytab / sum(mytab)
n <- length(mytab)
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Frequency Table (Categorical Data)',3,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Category',header=TRUE)
a<-table.element(a,'Abs. Frequency',header=TRUE)
a<-table.element(a,'Rel. Frequency',header=TRUE)
a<-table.row.end(a)
for (i in 1:n) {
a<-table.row.start(a)
a<-table.element(a,labels(mytab)$x[i],header=TRUE)
a<-table.element(a,mytab[i])
a<-table.element(a,round(reltab[i],4))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable1.tab')
}