[R] Accessing lme source code

Dear all; This an FAQ. I tried to access lme source script so I can step into it to debug the problems resulting from a lme() call. I used getAnywhere(lme) or nlme:::lme, both produced only the function definition and UseMethod(lme). Any idea how to list the source code?

[R] Variance-covariance by factor

Dear all, I have a data frame with one factor and four numeric variables and wish to obtain the var-cor matrix separately by factor. I tried by() and sapply() but getting nowhere. I understand this can be done by subsetting the dataframe, but there should have some sleek ways of doing

[R] Lattice bwplot error

Dear all; Searching the R site for answers to my problem, but found none. Here is the run and error: bwplot(dev ~ Dbhcl | Period, data = DbhValid2, font = 2, + main=list( ), axis.font =2, + ylab = list(label = Residual (cm), font = 2), + xlab = list(label = Dbh class (cm), font =

RE: [R] Lattice bwplot error

Thank, Deepayan Sarkar and Sundar Dorai-Raj for their quick responses. After updating lattice, the problem resolved. Richard -Original Message- Subject: [R] Lattice bwplot error Dear all; Searching the R site for answers to my problem, but found none. Here is the run and

[R] gnls(0 error: invalid variable type

Dear R-helpers; While using gnls() to fit a function Gbht0t.gnls - gnls(h2 ~ Rht(b0, b1, b2, h1,t1, t2), data=gbht10, + params=list(b0 + b1 + b2 ~ Sisp -1), start=c(strssb0,strssb1,strssb2)) I encountered an error: Error in model.frame(formula, rownames, variables,

[R] nls() error: Object not found

Dear R-helpers; Using nls() to fit a function,Rdum, defined below I stumbled on an error: Error in eval(expr, envir, enclos): Object s0 not found. The function Rdum is defined as Rdum - deriv(~ h1 * (s0 + sl0*sl + sm0*sm + sp01*sp1 + sp02*sp2 + sp03*sp3+sp04*sp4) *

[R] Plot.lme error

Dear All; Attempting to reproduce Figure 4.15 of MEMSS (p. 171) in R using plot(Wafer, outer = ~ Wafer) yields an error: plot(Wafer, outer = ~ Wafer) Error in order(na.last, decreasing, ...) : Argument lengths differ The plot() produces the figure without

[R] NLME augPred Error

Dear all; This error has been posted to the list a couple of times, but no answer / walk around solution is available. An example on P. 377 MEMSS plot(augPred(fm5CO2.nlme, levels=0:1), layout =c(6,2)) generates an error: Error in predict.nlme(object, value[1:(nrow(value)/nL), , drop =

RE: [R] [OFF] Dataset for extra Crawley Chapter

Ronaldo; You can download all datasets used in Crawley's book from his web site. -- Richard -Original Message- From: Ronaldo Reis Jr. [mailto:[EMAIL PROTECTED] Sent: Tuesday, October 14, 2003 6:47 AM To: R-Help Subject: [R] [OFF] Dataset for extra Crawley Chapter Hi,

[R] RE: Installing GLMMGibbs problems

Thank Andy Liaw for his response to my enquiry. He suggested creating the tar.gz source package at the parent directory of the fixed code using R CMD build GLMMGibbs and then R CMD INSTALL GLMMGibbs to install the package. The package is now installed and works fine. Thanks Andy, --

[R] Installing GLMMGibbs problems

Dear all; Installing the GLMMGibbs package to my Solaris Unix box, I got an compiling error: ars.c:497:10: missing terminating character ars.c: In function `dump_arse': ars.c:498: error: parse error before mylesj . The

[R] Rgui access violation

Dear All; While using EMclust() in the mclust package, I frequently encountered a program error. A message window popped up with the message Rgui.exe has generated errors and will be closed by Windows. You will need to restart the program. An error log is be created. version

[R] Trouble with optim

Dear All; Searching on the achieve, many questions on optim() have been asked, but I haven't seen the following. The question began with my original inquiry on Optimization failed in fitting mixture 3-parameter Weibul l distribution using fitdistr() which I posted on Jul. 28,

[R] Optimization failed in fitting mixture 3-parameter Weibulldistri bution using fitdistr()

Dear All; I tried to use fitdistr() in the MASS library to fit a mixture distribution of the 3-parameter Weibull, but the optimization failed. Looking at the source code, it seems to indicate the error occurs at if (res$convergence 0) stop(optimization failed).

[R] Winedt and R on Windows XP

Dear All; I have used Winedt editor in conjunction with R on the Win2K platform for more than a year without any problem. Recently I purchased a P4 machine with Windows XP. Following the installation of R1.62 and WinEdt 5.3 to separate directories, I copied R-WinEdt to the Plugins

[R] Mtext and xyplot

Dear all; I wish to create a graphic object combing an xyplot() and an mtext(). My code looks like following, gmv - { trellis.device(windows, bg=white, width = 7, height = 7) xyplot(Mvol ~ Age | Nl * Th , data = Hft1, prepanel = function(x, y) prepanel.loess(x, y,

[R] Multiple xyplot on a page

Dear all; I have a set of 5 xyplots each with 8 panels and wish to place 2 on a letter-sized page. I tried using pdf() to specify paper size and plot size: pdf(file = fig1b.pdf, paper = letter, onefile = FALSE, width = 5, height = 5, family=Helvetica) par(mfrow=c(1,2)) but it did not