similar to: Confidence intervals

Dear R graphics gurus, Another question about lattice graphics. This time I would like to plot means and confidence intervals by group factor in a lattice graph. I can not find any working lattice examples. Maybe a custom panel function is the answer, but that is a bit beyond me for now. The individual plots within the lattice graph could look like this: # Example with confidence intervals

Hi, Sorry for repeated question. I performed logistic regression using lrm and penalized it with pentrace function. I wanted to get confidence intervals of odds ratio of each predictor and summary(MyModel) gave them. I also tried to get bootstrapping standard errors in the logistic regression. bootcov function in rms package provided them. Then, I found that the confidence intervals provided by

Hi - I am using R version 2.13.0. I have run several GLMMs using the glmmPQL function to model the proportion of fish caught in one net to the total caught in both nets by length. I started with a polynomial regression full model with three length terms: l, l^2, and l^3 (l=length). The length terms and intercept were the fixed effects and the random effect was a paired haul (n=18).

You need: myCor <- function(data, index){ mine(data[index, ])$MIC[1, 2] } results=boot(data = cbind(C,D), statistic = myCor, R = 2000) boot.ci(results,type="all") Look at the differences between: mine(C, D) and mine(cbind(C, D)) The first returns a value, the second returns a symmetric matrix. Just like cor() David L. Carlson Department of Anthropology Texas A&M

Dear R-Experts, Thanks to Prof. Bonnett, I have got an R script working to calculate confidence intervals around the semipartial correlation coefficients. Now, I would like to calculate BCa bootstrap CIs using the boot library and the boot.ci(results, type="all") function. How could I modify my R script (here below reproducible example) to get the BCa bootstrap CIs ? CIsemipartcorr

Hi, I would be very grateful for advice on getting confidence intervals for the ordinary (non smoothed) parameter estimates from a gam. Motivation I am studying hospital outcomes in a large data set. The outcomes of interest to me are all binary variables. The one in the example here, Dead30d, is death within 30 days of admission. Sexf is gender (M or F), Age is age in years at the start

Hello, First of all, when I tried to use function mic I got an error. mic(cbind(C, D)) Error in mic(cbind(C, D)) : could not find function "mic" So I've changed your function myCor and all went well, with a warning relative to BCa intervals. myCor <- function(data, index){ mine(data[index, ])$MIC } results=boot(data = cbind(C,D), statistic = myCor, R = 2000)

Hi Rui, Many thanks. The R code works BUT the results I get are quite weird I guess ! MIC = 0.2650 Normal 95% CI = (0.9614, 1.0398) The MIC is not inside the confidence intervals ! Is there something wrong in the R code ? Here is the reproducible example : ########## C=c(2,4,5,6,3,4,5,7,8,7,6,5,6,7,7,8,5,4,3,2) D=c(3,5,4,6,7,2,3,1,2,4,5,4,6,4,5,4,3,2,8,9) library(minerva) mine(C,D)$MIC

This is a follow-up mail of "the effect of blocking on the size of confidence intervals - analysis using aov". In both mails I pursue the idea of using blocking factors in order to reduce the width of confidence intervals. My dataset comprises, a quantitative response variable, namely: "response", and three categorical eplanatory variables, namely: "method",

Dear R-Experts, Here below is my R code (reproducible example) to calculate the confidence intervals around the spearman coefficient. ########## C=c(2,4,5,6,3,4,5,7,8,7,6,5,6,7,7,8,5,4,3,2) D=c(3,5,4,6,7,2,3,1,2,4,5,4,6,4,5,4,3,2,8,9) cor(C,D,method= "spearman") library(boot) myCor=function(data,index){ cor(data[index, ])[1,2] } results=boot(data=cbind(C,D),statistic=myCor, R=2000)

Hello! I want to calculate simultaneous confidence intervals for a nominal variable with three categories: "yes", "no", "partially" and I expect that far more than 5 samples fall into each category. I have read that Glaz & Sison's method is only appropriate for variables with 7 or more categories. Therefore, the Goodman method seems like a good idea. I have

Hi Cristian, I don't know about a R routine for exact CIs, but I found a function called "ageadjust" using a gamma distribution approximation some time ago -take a look at http://medepi.org/epitools/rfunctions/index.html. That Webpage seems a bit outdated now, but the calculations given in the ageadjust functions should still work, I think. Maybe this helps; regards Heinrich.

Dear List I?m using the ?lmomRFA? package to fit different distributions to my data sample. To calculate the error bounds I used: regsimq(?) and sitequantbounds(?) So my questions are: Are error bounds and confidence intervals the same thing? And: Does regsimq(? boundprob = c(0.05, 0.95)) calculate the 90 or the 95% confidence interval? If error bounds and confidence intervals are

I would like to obtain Confidence Intervals for the estimates (unstandardized beta weights) of each predictor in a WLS regression: m1 = lm(x~ x1+x2+x3, weights=W, data=D) SPSS offers that output by default, and I am not able to find a way to do this in R. I read through predict.lm, but I do not find a way to get the CIs for multiple independent variables. Thank you Torvon [[alternative HTML

Hi, I'm trying to get confidence intervals to slopes from a linear model and I can't figure out how to get at them. As a cut 'n' paste example: ################# # dummy dataset - regression data for 3 treatments, each treatment with different (normal) variance x <- rep(1:10, length=30) y <- 10 - (rep(c(0.2,0.5,0.8), each=10)*x)+c(rnorm(10, sd=0.1), rnorm(10,

My earlier posting should have said "dotchart", not "barchart". 1)  How does one show confidence-intervals in a "dotchart" and use rownames for labels on the y-axes?  I have looked at "plotCI" in "gregmisc" package . But it does not seem to produce something like a dotchart.  The statistic, error, upper-bound, and lower-bound are in a

Sometimes one might like to obtain pointwise bootstrap bias-corrected, accelerated (BCA) confidence intervals for a large number of statistics computed from a single dataset. For instance, one might like to get (so as to plot graphically) bootstrap confidence bands for the fitted values in a regression model. (Example: Chiu S et al., Early Acceleration of Head Circumference in Children with

Hi R users, 1) How does one show confidence-intervals in a barchart and use rownames for labels on the y-axes? I have looked at "plotCI" in "gregmisc" package . But it does not seem to produce something like a barchart. The statistic, error, upper-bound, and lower-bound are in a dataframe. 2) How to show CI in a barchart either using the statistic and, either (a)

Hi all!! I am trying to plot several confidence intervals in a unique plot. That is, for each x, I have a confidence interval for a parameter related to x and I would like to plot them in the same plot, in order to compare them. The plot would look like some parallel vertical lines, each one corresponding to a x value. Their extrem points would be the confidence interval limits. I do not know if

Good Morning: I posted an initial question a few days ago and I received some good advice from two R experts. I have re-examined the Davison-Hinkley text paying close attention to the examples of the boot() and boot.ci() in that text and the single example of a similar process in the MASS book (not the MASS package manual as I initially misunderstood). I think I understand how the stratified