Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_samplenorm.wasp
Title produced by softwareMinimum Sample Size - Testing Mean
Date of computationSun, 09 Dec 2012 14:06:58 -0500
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2012/Dec/09/t1355080175wcd7h6sr7uvtqu7.htm/, Retrieved Thu, 31 Oct 2024 23:52:32 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=198036, Retrieved Thu, 31 Oct 2024 23:52:32 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact91
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Minimum Sample Size - Testing Mean] [Paper deel 2 test...] [2012-12-09 19:06:58] [fd3c35a156f52433b5d6e23e16a12a78] [Current]
Feedback Forum

Post a new message




\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input & view raw input (R code)  \tabularnewline
Raw Output & view raw output of R engine  \tabularnewline
Computing time & 1 seconds \tabularnewline
R Server & 'Sir Maurice George Kendall' @ kendall.wessa.net \tabularnewline
R Engine error message & 
Error in plot.window(...) : need finite 'ylim' values
Calls: plot -> plot.default -> localWindow -> plot.window
In addition: Warning messages:
1: In min(x) : no non-missing arguments to min; returning Inf
2: In max(x) : no non-missing arguments to max; returning -Inf
Execution halted
\tabularnewline \hline \end{tabular} %Source: https://freestatistics.org/blog/index.php?pk=198036&T=0

[TABLE]
[ROW][C]Summary of computational transaction[/C][/ROW]
[ROW][C]Raw Input[/C][C]view raw input (R code) [/C][/ROW]
[ROW][C]Raw Output[/C][C]view raw output of R engine [/C][/ROW]
[ROW][C]Computing time[/C][C]1 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Sir Maurice George Kendall' @ kendall.wessa.net[/C][/ROW]
[ROW][C]R Engine error message[/C][C]
Error in plot.window(...) : need finite 'ylim' values
Calls: plot -> plot.default -> localWindow -> plot.window
In addition: Warning messages:
1: In min(x) : no non-missing arguments to min; returning Inf
2: In max(x) : no non-missing arguments to max; returning -Inf
Execution halted
[/C][/ROW] [/TABLE] Source: https://freestatistics.org/blog/index.php?pk=198036&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=198036&T=0

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'Sir Maurice George Kendall' @ kendall.wessa.net
R Engine error message
Error in plot.window(...) : need finite 'ylim' values
Calls: plot -> plot.default -> localWindow -> plot.window
In addition: Warning messages:
1: In min(x) : no non-missing arguments to min; returning Inf
2: In max(x) : no non-missing arguments to max; returning -Inf
Execution halted



Parameters (Session):
par1 = 20000 ; par2 = 0.05 ; par3 = 0.95 ; par4 = 0.36 ; par5 = 0,8 ;
Parameters (R input):
par1 = 20000 ; par2 = 0.05 ; par3 = 0.95 ; par4 = 0.36 ; par5 = 0,8 ;
R code (references can be found in the software module):
par1 <- as.numeric(par1)
par2 <- as.numeric(par2)
par3 <- as.numeric(par3)
par4 <- as.numeric(par4)
par5 <- as.numeric(par5)
(z <- abs(qnorm((1-par3)/2)) + abs(qnorm(1-par5)))
(z1 <- abs(qnorm(1-par3)) + abs(qnorm(1-par5)))
z2 <- z*z
z2one <- z1*z1
z24 <- z2 * par4
z24one <- z2one * par4
npop <- array(NA, 200)
ppop <- array(NA, 200)
for (i in 1:200)
{
ppop[i] <- i * 100
npop[i] <- ppop[i] * z24 / (z24 + (ppop[i] - 1) * par2*par2)
}
bitmap(file='pic1.png')
plot(ppop,npop, xlab='population size', ylab='sample size (2 sided test)', main = paste('Confidence',par3))
dumtext <- paste('Margin of error = ',par2)
dumtext <- paste(dumtext,' Population Var. = ')
dumtext <- paste(dumtext, par4)
mtext(dumtext)
grid()
dev.off()
par2sq <- par2 * par2
num <- par1 * z24
denom <- z24 + (par1 - 1) * par2sq
(n <- num/denom)
num1 <- par1 * z24one
denom1 <- z24one + (par1 - 1) * par2sq
(n1 <- num1/denom1)
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Minimum Sample Size',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Population Size',header=TRUE)
a<-table.element(a,par1)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Margin of Error',header=TRUE)
a<-table.element(a,par2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Confidence',header=TRUE)
a<-table.element(a,par3)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Power',header=TRUE)
a<-table.element(a,par5)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Population Variance',header=TRUE)
a<-table.element(a,par4)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'z(alpha/2) + z(beta)',header=TRUE)
a<-table.element(a,z)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'z(alpha) + z(beta)',header=TRUE)
a<-table.element(a,z1)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Minimum Sample Size (2 sided test)',header=TRUE)
a<-table.element(a,n)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Minimum Sample Size (1 sided test)',header=TRUE)
a<-table.element(a,n1)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')
(ni <- z24 / (par2sq))
(ni1 <- z24one / (par2sq))
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Minimum Sample Size (for Infinite Populations)',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Population Size',header=TRUE)
a<-table.element(a,'infinite')
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Margin of Error',header=TRUE)
a<-table.element(a,par2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Confidence',header=TRUE)
a<-table.element(a,par3)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Power',header=TRUE)
a<-table.element(a,par5)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Population Variance',header=TRUE)
a<-table.element(a,par4)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'z(alpha/2) + z(beta)',header=TRUE)
a<-table.element(a,z)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'z(alpha) + z(beta)',header=TRUE)
a<-table.element(a,z1)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Minimum Sample Size (2 sided test)',header=TRUE)
a<-table.element(a,ni)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Minimum Sample Size (1 sided test)',header=TRUE)
a<-table.element(a,ni1)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')
(z <- abs(qt((1-par3)/2,n-1)) + abs(qt(1-par5,n-1)))
(z1 <- abs(qt(1-par3,n1-1)) + abs(qt(1-par5,n1-1)))
z2 <- z*z
z2one <- z1*z1
z24 <- z2 * par4
z24one <- z2one * par4
par2sq <- par2 * par2
num <- par1 * z24
denom <- z24 + (par1 - 1) * par2sq
(n <- num/denom)
num1 <- par1 * z24one
denom1 <- z24one + (par1 - 1) * par2sq
(n1 <- num1/denom1)
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Minimum Sample Size (Unknown Population Variance)',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Population Size',header=TRUE)
a<-table.element(a,par1)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Margin of Error',header=TRUE)
a<-table.element(a,par2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Confidence',header=TRUE)
a<-table.element(a,par3)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Power',header=TRUE)
a<-table.element(a,par5)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Population Variance',header=TRUE)
a<-table.element(a,'unknown')
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'t(alpha/2) + t(beta)',header=TRUE)
a<-table.element(a,z)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'t(alpha) + t(beta)',header=TRUE)
a<-table.element(a,z1)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Minimum Sample Size (2 sided test)',header=TRUE)
a<-table.element(a,n)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Minimum Sample Size (1 sided test)',header=TRUE)
a<-table.element(a,n1)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')
(z <- abs(qt((1-par3)/2,ni-1)) + abs(qt(1-par5,ni-1)))
(z1 <- abs(qt(1-par3,ni1-1)) + abs(qt(1-par5,ni1-1)))
z2 <- z*z
z2one <- z1*z1
z24 <- z2 * par4
z24one <- z2one * par4
(ni <- z24 / (par2sq))
(ni1 <- z24one / (par2sq))
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Minimum Sample Size
(Infinite Population, Unknown Population Variance)',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Population Size',header=TRUE)
a<-table.element(a,'infinite')
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Margin of Error',header=TRUE)
a<-table.element(a,par2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Confidence',header=TRUE)
a<-table.element(a,par3)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Power',header=TRUE)
a<-table.element(a,par5)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Population Variance',header=TRUE)
a<-table.element(a,'unknown')
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'t(alpha/2) + t(beta)',header=TRUE)
a<-table.element(a,z)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'t(alpha) + t(beta)',header=TRUE)
a<-table.element(a,z1)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Minimum Sample Size (2 sided test)',header=TRUE)
a<-table.element(a,ni)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Minimum Sample Size (1 sided test)',header=TRUE)
a<-table.element(a,ni1)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')