similar to: mblm, confint problem

Hello, I am working on a simple non-parametric (Theil) regression function and and am following Hollander and Wolfe 1999 text. I would like some help making my function faster. I have compared with pre-packaged version from "MBLM", which isnt very fast either, but it appears mine is faster with N = 1000 (see results below). I plan on running this function repeatedly, and I generally

I haven't seen this cryptic warning before: > update.views('Robust') Warning message: In update.views("Robust") : The following packages are not available: covRobust, distr, FRB, MASS, mblm, multinomRob, mvoutlier, quantreg, RandVar, rgam, RobAStBase, robfilter, RobLox, RobRex, robust, RobustAFT, robustbase, ROptEst, ROptRegTS, rrcov, sandwich, wle >

Hi All... I haven?t found mention of this error anywhere. I'm trying to draw spline curves using grid graphics. Most of the time, I have no problems, but I have some data sets that give the error in the subject line. I'm not sure which end points are identical, but the end points passed to the function are definitely not identical. Any assistance appreciated! Bryan tst <-

Hello, I am analyzing some climate time series data using the Mann Kendall package and was wondering if there was a way to calculate the trend using Sen's nonparametric estimator slope in R? Thank you in advance, Barry _________________________ Barry Baker, Ph.D. Global Climate Change Initiative The Nature Conservancy 2424 Spruce St., Suite 100 Boulder, CO 80302 Tel: (303)-541-0322 Fax:

Hi All Please see the Rprofile file which i have modified as follows and after that when I start R then I see that R says to me "TRUE" for all the packages implying that all loaded at once. But when i try to use commands as simple as help("lm"), it doesnt work nor any of the menu "Packages" is not working. Although the regression using lm ( Y ~ X ) is working

>>>>> 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

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

I don't know if this has already come up in the list or elsewhere - a quick search did't show anything relevant - but I think it's worth of mention. The confint.lm function in package MASS doesn't work correctly when called on a subset of parameters. The bug, easy to fix, is that confidence intervals are computed for all parameters anyway, and then assigned to a matrix which is too

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 sample program > set.seed(7882) > x<-rep(c(0,1),c(20,20)) >

On July 12, 2004 Spencer Graves wrote an email describing essentially the same issue that I would like help on: calling the confint function from within another homemade function. Because he provided many good examples of the problem, I will not reproduce them here but will instead refer readers to the original posting: http://www.mail-archive.com/r-help at stat.math.ethz.ch/msg23826.html It is

Dear colleagues, I have used statistical software "R" in academic research, and I find very interesting. So now I have started my studies at advanced level. I have conducted several statistical models of the type "*nls*". In respect of procedures for models of the type nls, I'd like you answer only one question that I have had doubts, related to the formula

Hi all, I've got a list (Reg3Lst) with 1000 nls regression results in it. I'd like to get the confidence interval of the parameters obtained with the nonlinear regressions. Thus I've used this: for (i in 1:1000) { foo<-list(foo,confint(Reg3Lst[[i]])) } For some regressions the confidence interval is not estimated because of a singular gradient result. Then it

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 find any way to turn off 'Waiting for. .." I tried

Hi there, I had a problem when I hoped to get confidence intervals for the parameters I got using mle() of stats4 package. This problem would not appear if ``fixed'' option was not used. The following mini-example will demo the problem: x <- c(100, 56, 32, 18, 10, 1) r <- c(18, 17, 10, 6, 4, 3) n <- c(18, 22, 17, 21, 23, 20) loglik.1 <- function(alpha, beta, c) { x

confint.glm prints a message "Waiting for profiling to be done..." I could have sworn that there used to be an option (quiet = TRUE?) to turn this message off without resorting to suppressMessages() (finer/more specific control is always preferable ...) -- but on the basis of looking back at archived versions of MASS, and at this Stack Overflow post:

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() version, if any such help page exists. Is there a general way of turning of output

I fit a proportional odds model with the polr-function of the MASS package from Venables and Ripley Applying the confint method to calculate confidence intervals for the parameters I get the following error message Waiting for profiling to be done... Re-fitting to get Hessian Error in X[, -i, drop = FALSE] : incorrect number of dimensions Can someone explain the error-message? (The

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 a vector of numbers or

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

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, family = binomial, data = DF, weights = DF[, "N"]) Coefficients: