Displaying 6 results from an estimated 6 matches for "lowerci".
Did you mean:
lower
2009 Jan 15
1
Confidence Intervals for Poisson
Hi folks!
I run the following code to get a CI for a Poisson with lambda=12.73
library(MASS)
set.seed(125)
x <- rpois(100,12.73)
lambda_hat<-fitdistr(x, dpois, list(lambda=12))$estimate
#Confidence Intervals - Normal Approx.
alpha<-c(.05,.025,.01)
for(n in 1:length(alpha)) {
LowerCI<-mean(x)-(qnorm(1-alpha[n]/2, mean = 0, sd =
1)*sqrt(var(x)/length(x)))
UpperCI<-mean(x)+(qnorm(1-alpha[n]/2, mean = 0, sd =
1)*sqrt(var(x)/length(x)))
cat("For
Alpha=",alpha[n],"LowerCI=",LowerCI,"<","Lambda=",mean(x),"<","UpperC...
2010 Jan 03
3
R2HTML Report number format, or Better Way?
...to produce a report and here's where I get stuck:
How do I get R2HTML to produce the same number format?
Particularly remove the decimal places for Par and Sal.
Are there better methods to produce this type of report?
Thanks,
L.A.
R version 2.10.0 XP
srtype<-cbind(Par,Sal,Median,COD,PRD,LowerCI,UpperCI)
srtype
Par Sal Median COD PRD LowerCI UpperCI
RES I 9683 578 0.9533 29.69 1.191 0.9382 0.9582
RES V 4003 155 0.9763 16.51 1.091 0.9499 0.9943
OTHER 1542 10 0.8367 82.35 1.253 0.4759 2.2293
COM I...
2000 Jul 05
2
couldn't find function "..."
...I created my function named "CIfunc" and saved into the current
directory using the "save" command.
function(Mean, SD, N, Sig) {
if (N > 30) { Zscore <- pnorm(Sig)} else {Zscore <- pt(Sig,N-1)}
Zscore
UpperCI <- Mean + Zscore * SD / sqrt(N)
UpperCI
LowerCI <- Mean - Zscore * SD / sqrt(N)
LowerCI
}
...with no error.
Next I write the function call "testcall":
testcall <- call('CIfunc', prescan10$F532Mean, prescan10$F532SD,
prescan10$FPixels, 0.95)
...with no error.
testcall CIfunc(c(189, 238, 363, 869, 585, 350, 3624, 2...
2000 Jul 06
0
Fwd: Re: couldn't find function ...
...save" command.
> >
> > function(Mean, SD, N, Sig) {
> >
> > if (N > 30) { Zscore <- pnorm(Sig)} else {Zscore <- pt(Sig,N-1)}
> > Zscore
> >
> > UpperCI <- Mean + Zscore * SD / sqrt(N)
> > UpperCI
> >
> > LowerCI <- Mean - Zscore * SD / sqrt(N)
> > LowerCI
> >
> > }
> >
> > ...with no error.
> >
> > Next I write the function call "testcall":
> >
> > testcall <- call('CIfunc', prescan10$F532Mean, prescan10$F532SD,
> >...
2009 Jul 16
0
how to get means and confidence limits after glmmPQL or lmer
...imating the CIs.
variance<-mm$sigma #IS THIS THE RIGHT VARIANCE?
speciessize<- tapply(species, burnagecat1 ,length)
speciesmeans<- (tapply(predict(mm, type="response"), burnagecat1 ,mean))
upperci<- exp(log(speciesmeans) +
qt(0.95,residdf)*sqrt(variance/speciessize))
lowerci<- exp(log(speciesmeans)- qt(0.95,residdf)*sqrt(variance/speciessize))
result<-cbind(lowerci, speciesmeans, upperci)
result
I also explored lmer, but haven't found a way to extract means and
CIs on the original scale for the levels of treatment.
library(lme4)
dlm<-lmer(response~ tr...
2002 Apr 22
2
lattice x(y)lab and expression
Another question about lattice:
Is there a way to use plotmath in the labels and strips?
E.g.
xyplot(y~x|group,data=test.df,
panel=function(x,y,subscripts,df) {
llines(x,y,col="black",lwd=2)
llines(x,df[subscripts,"lowerCI"],lty=2,col="#aaaaaa")
llines(x,df[subscripts,"upperCI"],lty=2,col="#aaaaaa")
},
strip=function(...) {
# not sure what to put here
# but I would like something like
# substitute(gamma[r]==2*pi*r,list(r=group))...