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

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

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

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,

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

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

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

Hi, there might be a bug in read.delim (also read.table) in R.1.6.2. I included the correct input done in R1.6.1 first, then the wrong input in R1.6.2. The discrepancies are highlighted in red. Michael Man ####### correct input in R1.6.1 ########################## > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386,

I am running R1.6.2 under windo$e98. I have followed the instructions in the Win port faq concerning tck, but when I try 'library(tcltk)', I get Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load shared library "C:/PROGRA~1/R/RW1062/library/tcltk/libs/tcltk.dll": LoadLibrary failure: One of the library files needed to run this application

R1.6.2 winnt. file.choose() gives an error if the dialog is cancelled without choosing a file. This can be easily dealt with via try(), but I would have thought that a NULL or NA return would be the expected behavior. Cheers, Bert Gunter Biometrics Research RY 84-16 Merck & Company P.O. Box 2000 Rahway, NJ 07065-0900 Phone: (732) 594-7765 mailto: bert_gunter@merck.com "The business of