Displaying 2 results from an estimated 2 matches for "p1n".
Did you mean:
p1
2010 Jul 18
2
loop troubles
...The data below is dummy data, I used a starting
mu of 1, and starting precision of 0.
bayes.analysis.treat<-function(mu0,p0){
n1 = 5
n2 = 10
n3 = 15
ybar1 = 12
ybar2 = 13
ybar3 = 14
sd1 = 2
sd2 = 3
sd3 = 4
#posterior 1
var1 = sd1^2 #sample variance
p1 = n1/var1 #sample precision
p1n = p0+p1
mu1 = ((p0)/(p1n)*mu0)+((p1)/(p1n))*ybar1
sigma1 = 1/sqrt(p1n)
#posterior 2
var2 = sd2^2 #sample variance
p2 = n2/var2 #sample precision
p2n = p1n+p2
mu2 = ((p1n)/(p2n)*mu1)+((p2)/(p2n))*ybar2
sigma2 = 1/sqrt(p2n)
#posterior 3
var3 = sd3^2 #sample variance
p3 = n3/var3 #sample pr...
2007 Jun 08
1
Need Help with robustbase package: fitnorm2 and plotnorm2
...oading required package: RColorBrewer
Loading required package: grid
Loading required package: geneplotter
Loading required package: annotate
KernSmooth 2.22 installed
Copyright M. P. Wand 1997
[1] "/usr/local/lib/R/site-library/prada/extdata/fas-Bcl2-plate323-04-04.A01"
[1] 2115 8
$P1N $P2N $P3N $P4N $P5N $P6N $P7N $P8N
"FSC-H" "SSC-H" "FL1-H" "FL2-H" "FL3-H" "FL2-A" "FL4-H" "Time"
Loading required package: rrcov
Error in fitNorm2(fdat[, "FSC-H"], fdat[, "SSC-H&quo...