search for: confint

...Error z value Pr(>|z|) (Intercept) -2.03693352 1.449574526 -1.405194 0.159963578 A 0.01093048 0.006446256 1.695633 0.089955471 N 0.41060119 0.224860819 1.826024 0.067846690 S -0.20651005 0.067698863 -3.050421 0.002285206 then I use confint(k.glm) to obtain a confidnece interval for the estimates. > confint(k.glm,level=0.97) Waiting for profiling to be done... 1.5 % 98.5 % (Intercept) -5.471345995 0.94716503 A -0.002340863 0.02631582 N -0.037028592 0.95590178 S -0.365570347...

I can't get confint.glm to work from within a function. Consider the following (using R 1.9.1, Windows 2000): # FIRST: SOMETHING THAT WORKS FROM A COMMAND PROMPT DF <- data.frame(y=.1, N=100) (fit <- glm(y~1, family=binomial, data=DF, weights=DF[,"N"])) Call: glm(formula = y ~ 1, famil...

...commented in detail here<https://github.com/tidymodels/broom/issues/714>): # object with multiple classes: the output of function `geepack::geeglm`. The output of `class(object)`: ``` [1] "geeglm" "gee" "glm" "lm" ``` The generic function: `stats::confint`. The S3 method chosen: `confint.glm`. It produces an error. The correct method in this case would be `broom:::confint.geeglm`. Thank you! Iago [[alternative HTML version deleted]]

It seems that confint.default returns an empty data.frame for objects of class mlm. For example: ``` nobs <- 20 set.seed(1234) # some fake data datf <- data.frame(x1=rnorm(nobs),x2=runif(nobs),y1=rnorm(nobs),y2=rnorm(nobs)) fitm <- lm(cbind(y1,y2) ~ x1 + x2,data=datf) confint(fitm) # returns: 2.5 % 97.5 %...

Hi, I want to view the code for the function confint in the stats package. Typing just confint doesn't work. How can I view the code? Thanks, Walt ________________________ Walter R. Paczkowski, Ph.D. Data Analytics Corp. 44 Hamilton Lane Plainsboro, NJ 08536 ________________________ (V) 609-936-8999 (F) 609-936-3733 walt at dataanalyticsco...

Dear all, I was trying to see how the function 'confint' is defined. Doing > confint function (object, parm, level = 0.95, ...) UseMethod("confint") <environment: namespace:stats> does not really enlighten me. How can I get to see the implementation (I guess it should be possible according to the general philosophy of the R pro...

I need confint() for glm() to supress the messages "Waiting for profiling to be done..." because they mess up the caching mechanism of pgfSweave (see https://github.com/cameronbracken/pgfSweave/issues/40). I have read the help page of confint(), but I do not know how to get the help page for the glm...

...summary(sib.glm) The coefficients in the output are Estimate Std. Error z value Pr(>|z|) (Intercept) -3.6826 0.5062 -7.274 3.48e-13 *** sex 2.2059 0.5380 4.100 4.13e-05 *** I have calculated the .95 confidencce interval for sex two ways: (1) confint(sib.glm) The result is 2.5 % 97.5 % (Intercept) -4.861153 -2.823206 sex 1.263976 3.428764 Using the usual confidence interval formula, (2) 2.2059 +/- 1.96*.538 = 1.15142. 3.26038 The results from (2) are identical to those from SPSS but do not agree with thos...

>>>>> steven pav >>>>> on Thu, 19 Jul 2018 21:51:07 -0700 writes: > It seems that confint.default returns an empty data.frame > for objects of class mlm. For example: > It seems that confint.default returns an empty data.frame for objects of > class mlm. Not quite: Note that 'mlm' objects are also 'lm' objects, and so it is confint.lm() which is called her...

Hi, I have a question about the computation of confidence intervals in the zyp package, in particular using the functions zyp.sen and confint.zyp, or zyp.yuepilon. (1) I'm a bit confused about the confidence intervals given by zyp.sen and confint.zyp. When I request a certain confidence interval in the function, the R output seems to deliver another confidence interval, e.g. when I set level=0.95 in the function, then the output is...

A non-text attachment was scrubbed... Name: not available Type: text/enriched Size: 3706 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/20010913/b4135232/attachment.bin

Dear R developers, I noticed a bug in the stats4 package, specifically in the confint method applied to ?mle? objects. In particular, when some ?fixed? parameters define the log likelihood, these parameters are stored within the mle object but they are not used by the ?confint" method, which retrieves their value from the global environment (whenever they still exist). Sample...

the function confint uses the profiling method of the function of the package MASS confint.glm even after the package has been detached! 1: might this be the intenden behavior? 2. How does the function remember its 'MASS' functionality after detaching the package? R: 1.8.0; Windows 2000 Here is a sampl...

Dear R-helpers, I am using confint() within a function, and I want to turn off the message it prints: x <- rnorm(100) y <- x^1.1+rnorm(100) nlsfit <- nls(y ~ g0*x^g1, start=list(g0=1,g1=1)) > confint(nlsfit) Waiting for profiling to be done... 2.5% 97.5% g0 0.4484198 1.143761 g1 1.0380479 2.370057 I cannot...

Dear R users, The nnet package includes the multinom method for the confint function. The R Help file (?confint) for the generic function in the stats package and the help files for the glm and nls methods in the MASS package indicate that one can use the "parm" argument as "a specification of which parameters are to be given confidence intervals, either...

Full_Name: Christian Lajaunie Version: 2.5.1 OS: Fedora fc6 Submission from: (NULL) (193.251.63.39) confint() does not use the appropriate variance term when the design matrix contains a zero column (which of course should not happen). Example: A 10x2 matrix with trivial column 1: > junk <- data.frame(x=rep(0,10), u=factor(sample(c("Y", "N"), 10, replace=T))) The response: &...

....103.21) The print-out of xtable in the following example is wrong; instead of yielding the correct ci's for the second model it repeats the ci's from the first model. require(xtable) require(MASS) data(cats) b1<-lm(Hwt~Sex,cats) b2<-lm(Hwt~Sex+Bwt,cats) cbind(c(coef(b1),NA),rbind(confint(b1),c(NA,NA)),coef(b2),confint(b2)) xtable(cbind(c(coef(b1),NA),rbind(confint(b1),c(NA,NA)),coef(b2),confint(b2)))

Hi, After installing the latest versions of lme4, Matrix, VR on a Debian Box, I have run into a problem. When I use library(lme4) on R Version 2.1.0 (2005-04-18) I get the following error message: > library(lme4) Loading required package: Matrix Error in autoloader(name = "confint", package = "MASS") : autoloader did not find 'confint' in 'MASS' An example of the complete session is given at the end of this email. A search in the docs. At present I have 3 directories where packages are kept, as I use apt-get install "r-cran-pack...

...tercept) -0.06106 0.29808 -0.205 0.8377 pSusSD 0.21184 0.36886 0.574 0.5658 pBenSD 1.20255 0.52271 2.301 0.0214 * pBarSD -0.08654 0.48749 -0.178 0.8591 pSevSD 0.99759 0.44795 2.227 0.0259 * And this is are the corresponding CI when I call the confint function: 2.5 % 97.5 % (Intercept) -0.6548023 0.5264357 pSusSD -0.4980888 0.9733975 pBenSD 0.2665235 2.3495259 pBarSD -1.0695945 0.8740359 pSevSD 0.1877044 1.9747499 Utilizing the formula I mentioned above, the correct CI for pSusSD would actually be: >...

Dear R users this may be a simple question - but i would appreciate any thoughts does anyone know how you would get one lower and one upper confidence interval for a set of data that consists of proportions. i.e. taking a usual confidence interval for normal data would result in the lower confidence interval being negative - which is not possible given the data (which is constrained between