similar to: regression parms var-cov matrix

Hi list, can I extract the cov.unscaled ("the unscaled covariance matrix") from a gls fit (package nlme), like with summary.lm? Background: In a fixed effect meta analysis regression the standard errors of the coefficients can be computed as sqrt(diag(cov.unscaled)) where cov.unscaled is (X'WX). I try do do this with a gls-fit. Thanks, Sven

R 1.6.2 for Windows, Win2k: I have fitted a weighted least squares model using the code "wls.out <- gls(y ~ x1 + x2 + x3 + x4 + x5 + x6 - 1, data = foo.frame, weights = varConstPower(form = ~ fitted(.), fixed = list(power = 0.5), const = 1))" The data has 62 rows and the response is zero when the covariates are zero. The purpose of the model was to account for the the fact that

Hi, I am running a simple linear model with (say) 5 independent variables. Is there a simple way of getting the variance-covariance matrix of the coeffcient estimates? None of the values of the lm() seem to provide this. Thanks in advance, Ritwik Sinha rsinha@darwin.cwru.edu Grad Student Case Western Reserve University [[alternative HTML version deleted]]

Hello, Is there a way to get the Var-Cov matrix from the LMER function? Thanks [[alternative HTML version deleted]]

Dear all, I'm trying to get the variance/covarince matrix after fitting a quantile regression model (either linear or non linear), in order to get the variance of my predictions and be able to calculate the median squared error. The commands working for the lm models (corr=T or vcov=T) do not seem to work for the rq models. Could you advise me a way of getting it? Best regards

This is my first post to this list so I suppose a quick intro is in order. I've been using SPLUS 2000 and R1.6.2 for just a couple of days, and love S already. I'm reading MASS and also John Fox's book - both have been very useful. My background in stat software was mainly SPSS (which I've never much liked - thanks heavens I've found S!), and Perl is my tool of choice for

Folks I have modified an existing function to calculate 'ec/ld/lc' 50 values and their associated Fieller's confidence limits. It is based on EC50.calc (writtien by John Bailer) - but also borrows from the dose.p (MASS) function. My goal was to make the original EC50.calc function flexible with respect to 1) probability at which to calculate the expected dose, and 2) the link

Dear list, I'd like to create a visual plot of a clmm() I've fitted using the 'ordinal' package in R. It's possible to do this with a glm() by using the 'effects' package. For example: library(effects) data(BEPS) mod <- lm(political.knowledge ~ age + gender + vote, data=BEPS) eff <- effect("age", mod, default.levels=100) plot(eff,

I am running a 3-parameter nonlinear fit using the default Gauss-Newton method of nls. initialValues.L = list(b=4,d=0.04,t=180); fit.nls.L = nls( myModel.nlm , fData.L, start = initialValues.L, control = nls.control(warnOnly = TRUE), trace=T ); summary.nls.L = summary(fit.nls.L); I run the same analysis in SAS proc NLIN. proc nlin data=apples outest=a; parms b=4 d=.04 t=180; model Y =

Dear friends, After fitting a generalized linear models ,i hope to get the variance of alpha,variance of beta and their covariance, that is , the variance-covariance matrix/information of alpha and beta , suppose *B* is the object of GLMs, i use attributes(B) to look for the options ,but can't find it, anybody knows how to get it? > attributes(B) $names [1] "coefficients"

I have had a visit from our windows buddy the BSOD twice in one day when running R1.6.2 under ESS5.1.24 on a win2000 platform. I wasn't doing anything special, just moving or resizing the emacs window, although I did have 4 or 5 emacs windows open at the time. Has anyone else observed this problem? Thanks, Darryl.

I've seen that there's an history for graphics. How can I save it manually? Can I save more than 4 plots? I'm using R1.6.2 under win98 [[alternative HTML version deleted]]

hi, for working with the simseg/acm approach i need multilm and it seems that the last version is 0.1-4.tar ? What's wrong in the description file, if i attempt install the package ? c:\DataMining\rw1062\bin>Rcmd Install multilm Malformed DCF file (file multilm/DESCRIPTION, line 6) [R1.6.2 /W2K] P.S. ..hmm a solution might be import the functions from multilm in my .Rprofile, but i

Hi, I have installed R1.6.2 on a Sun Sparc workstatioin 60 under Solaris 8. I' m using csh, tcsh or bash- shell. Unfortunately, it is impossible to use the "arrow-up" key on Sun's German keyboard to repeat a R command by using this key. Always following syntax error occurs: ^[[A Does anybody have a hint ? Thanks Klaus-P. Pleissner -- Dr. Klaus-Peter Pleissner Max Planck

Dear all, <BR><BR>I was used to work on R1.6 and I have now passed on R1.6.2 but I can''t read my <BR>file (and that is a big problem!!).<BR>I made a data sheet with some <BR>spreadsheet in Excell, and save it as separeted by tab .txt.<BR>I write in R<BR>read.table ("file.txt",h=T,sep="/t",dec=",")<BR>But R

Hi all, Here's a curosity I ran into since upgrading to 1.6.2 (precompiled for Windows). When using the chm help I get the following warning. I saw a recent post regarding this as a new warning (``dyn.load warning message in R1.6.2 on Windows XP'' dated 1/28/03), but not in the context of the help system. The warning only appears once and does not prevent the chm file from

Dear Experts, In glim() of S Plus, one of the returned values is "var", the estimated variance matrix of coefficients. However, in glm() of R (there is no glim() in R), "var" is not one of the returned values. Anyone know what could I get the varience matrix of coefficients in glm() in R? As a novice in R and S+, I'd appreciate your help Sincerely, Charlie Liu

Dear r-bugs, I looked over the FAQ. Hope I'm reporting this correctly. I ran this on both solaris and windows. I've provided terminal snapshots which include how R was called from the command line, and the result of version at the R prompt. I have attached the .r file, and the data file and the output snapshots. Below also find everything except only a few lines of the data file. Note

I am trying to do logistic regression with a categorical predictor variable with the glm() function, family=binomial. Using glm() I would like to be able to calculate the confidence intervals of all three possible odds ratios for a factor (the factor has three categories). Three categories imply two columns of 0's and 1's in the design matrix, and two parameter estimates with their

Hi, I would like to be able to access the calculated results from the nls package. Using the example in R, fm3DNase1 we can reurn certain parts of the calculations: > coef(fm3DNase1) Asym xmid scal 2.345179 1.483089 1.041454 > resid(fm3DNase1) [1] -0.0136806237 -0.0126806237 0.0089488569 0.0119488569 -0.0025803222 [6] 0.0064196778 0.0026723396 -0.0003276604