similar to: studentized CIs for a correlation using package boot

Hi, I'm encountering a very odd problem with calls to anova.mlm() from within a function. Consider the following code (data.n is a matrix of numeric values): mlmfit <- lm(data.n ~ 1) mlmfit0 <- lm(data.n ~ 0) print(mlmfit) anova(mlmfit,mlmfit0,test="Spherical") If I run it just like this from the console, it works just fine. If, however, I call it from within a function,

Hi all, I have been worriyng and googling a lot but I can't find my mistake. I am trying to regiter an X-Lite Softphone to Asterisk, but I am getting a SIP/2.0 403 Forbidden response: SEND TIME: 10157385 SEND >> 10.100.249.12:5060 REGISTER sip:10.100.249.12 SIP/2.0 Via: SIP/2.0/UDP 10.100.249.86:5060;rport;branch=z9hG4bKFAC1B6F2B5414EE9855696A09A83FB22 From: Tester

Hi, yet another anova.mlm problem - it doesn't seem to end. This time, I have a setup with a few within-subject factors and a between-subject factor (SGROUP). Consider the most simple case with only one within-factor (apo): > mlmfit0 <- lm(data.n ~ 0 + SGROUP) > mlmfit1 <- lm(data.n ~ 1 + SGROUP) > anova(mlmfit1,mlmfit0,test="Spherical",M=~hemi,X=~1) Analysis of

Hello, i have a problem with the function nls(). This are my data in "k": V1 V2 [1,] 0 0.367 [2,] 85 0.296 [3,] 122 0.260 [4,] 192 0.244 [5,] 275 0.175 [6,] 421 0.140 [7,] 603 0.093 [8,] 831 0.068 [9,] 1140 0.043 With the nls()-function i want to fit following formula whereas a,b, and c are variables: y~1/(a*x^2+b*x+c) With the standardalgorithm

I would have expected the function getAnywhere to have behaved differently in the following: > search() [1] ".GlobalEnv" "file:C:/R/Rdata/miya/.Rdata" [3] "package:boot" "package:methods" [5] "package:ctest" "package:mva" [7] "package:modreg"

Dear Friends and Colleagues, I'm puzzling over how to interpret or use some bootstrap intervals. I think that I know what I should do, but I want to check with knowledgeable people first! I'm using a studentized non-parametric bootstrap to estimate 95% confidence intervals for three parameters. I estimate the variance of the bootstrap replicates using another bootstrap. The script

I am trying to find studentized bootstrap intervals for the skewness of a data set. I have tried the following (nerve.dat is a set of data containing observations on one variable) (using Windows XP): x <- scan("e:/Flashbackup2009/Nonparametrics/nerve.dat") n <- length(x) library(e1071) skewness(x) library(boot) sampleskew <- function(x,d) {return(skewness(x[d]))} bb <-

Dear friends, In s-plus, lm() generates the the studentized residuals automatically for us, and In R, it seems don't have the results: After i fitted lm(), i use attibutes() to see the objects and didn't find studentized residuals . How to get the the studentized residuals in lm(),have i missed something? thanks very much! -- Kind Regards, Zhi Jie,Zhang ,PHD Department of Epidemiology

Hi all, I'm trying to apply loess regression to my data and then use the fitted model to get the *standardized/studentized residuals. I understood that for linear regression (lm) there are functions to do that:* * * fit1 = lm(y~x) stdres.fit1 = rstandard(fit1) studres.fit1 = rstudent(fit1) I was wondering if there is an equally simple way to get the standardized/studentized residuals for a

Has anyone written a function to find the quantiles for the studentized maximum distribution, and the studentized maximum modulus distribution? Thanks for any help. Kostas --------------------------------- [[alternative HTML version deleted]]

Dear R-Help, I am using the studres function from the MASS package to compute studentized deleted residuals in a oneway anova. I'm having trouble interpreting the results in situations where a factor level has only one observation. Sometimes studres yields an NaN and sometimes it produces a numeric value for cases where a factor level has only one observation. I would think it should

I'm sorry to waste bandwidth. I re-read the console message for the tenth time and finally noticed R was looking in rw2001pat for libraries. Looking at my Env vars I see I set R_LIBS to look there. Changed R_LIBS, fixed problem. DOH! So now if I want to use several versions of R simultaneously, what do I do. I set R_LIBS so it would look also in C:/R/extra for some added packages. Thanks and

Hi, I must be doing something silly here, because I can't get the studentised and standardised residuals from r output of a linear model to agree with what I think they should be from equation form. Thanks in advance, Jennifer x = seq(1,10) y = x + rnorm(10) mod = lm(y~x) rstandard(mod) residuals(mod)/(summary(mod)$sigma) rstudent(mod)

Can anyone tell me how to calculate quantiles of the studentized maximum modulus distribution using R? Thanks! __________________ Rob Kabacoff, Ph.D. www.statmethods.net [[alternative HTML version deleted]]

Hi, I am wondering if there is a way to extract studentized residues from a nls() object? I searched archive, someone posted a similar question before (http://tolstoy.newcastle.edu.au/R/e6/help/09/04/10845.html), but didn't get an answer. Thanks for any suggestions, John

Dear colleagues, hi. I am using the function gls in R-package nlme and I want to compute either the leverage values or the studentized residuals. Can any body tell me how can do that? With regards ___________________________________________________________ 围蟻畏蟽喂渭慰蟺慰喂蔚委蟍[elided Yahoo spam]] 螔伪蟻蔚胃萎魏伪蟿蔚 蟿伪 蔚谓慰蠂位畏蟿喂魏维 渭畏 未喂伪胃苇蟿蔚喂 蟿畏谓 魏伪位蠉蟿蔚蟻畏 未蠀谓伪蟿萎 蟺蟻慰蟽蟿伪蟽委伪 魏伪蟿维 蟿蠅谓 蔚谓慰蠂位畏蟿喂魏蠋谓 渭畏谓蠀渭维蟿蠅谓

Hi there, I hope I can get advice regarding the calculation of leverage values or studentized residual values of a non-linear regression model. It seems like rstudent() does not work on a nls object. Many thanks in advance! Best regards, Xingli

Well, it is clear enough that the problem is in interpreting the documentation. However, when you claim you tested something, and found it inconsistent with tables, it would be advisable to back it up with examples! The description in the help files and in the sources is admittedly confusing. The original paper has this, rather more clear, description in the abstract: "We consider the

We cannot help you understand what you are doing if you do not show us what you are doing. Here are some discussions about how to communicate questions about R [1][2][3]. [1] http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example [2] http://adv-r.had.co.nz/Reproducibility.html [3] https://cran.r-project.org/web/packages/reprex/index.html -- Sent from my phone.

Hi, Can anyone tell me why I am not getting the correct intervals for fixed effect terms for the following generalized linear mixed model from HPDinterval: > sessionInfo() R version 2.4.1 (2006-12-18) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252