similar to: parallel REML computation

Hi, [I'mm cc-ing r-help on this, since some others might be interested (Tony Rossini is, at least)] Yes, I saw your mail and was wondering about what would be the best approach to do it. It's definitely not trivial, mainly because Trellis is not designed to handle factors on the x-axis (chiefly for convenience, as factors typically have long names which are difficult to write on the

I have the same question that Eusebio had: Is there a function similar to "sink" that redirect also R code to a file that is: sink("R001") x <- c(2,-6,-4,8,5,4,1,3,4,-9,0,1) A <- matrix(x, ncol=3) A A.prima <- t(A) A.prima dim(A) dim(A.prima) sink() create a file "R001" with contents: ------------------------------------------ [,1] [,2] [,3] [1,]

Hi, I have just downloaded R (v. 0.90.1) and am trying to compile it on a Sun/UNIX system running Solaris 2.6. Needless to say, I have run into a problem, and cannot find any info on it through the CRAN site, or searching the old discussion files. When I do the "make', after configuring, I get the message, "The current directory must be set to the RSI directory. Change the

Hello All, The attached patch allows openssh to specify which pam service name to authenticate users against by specifying the PAMServiceName attribute in the sshd_config file. Because the parameter can be included in the Match directive sections, it allows different authentication based on the Match directive. In our case, we use it to allow different levels of authentication based on the

When I generate levelplot figures with Postscript, the lines from the underlying grid (i.e., the x and y in z ~ x* y) show up as faint white lines in the output. If I have a dense grid, the lines wash out the contour colors. I can avoid the problem if I use png instead, but that brings up other problems. Any ideas for how I can get Postscript to not show those lines around the cells? Scott

I wanted a simple function for bilinear interpolation on a 2-D grid, and interp.surface() in the fields package didn't quite suit my needs. In particular, it requires uniform spacing between grid points. It also didn't have the "visual" reference frame I was looking for. Here is an alternative function, followed by an example. # A function for bilinear interpolation on a 2-d

I can't get Emacs to automatically do syntax highlighting of Sweave files. I have followed Friedrich's suggestion for code to insert into my .emacs file. The complete section from my .emacs file is given below. When I load a *.Snw file, font is white until I press M-x, then the first code and document chunks get highlighted, but not the rest of the file. Latex and Noweb menus are

I searched CRAN for material on gant (or Gantt) charts, those schedule plots that look like horizontal barplots where the x-axis is time. Finding none, I wrote the following and am submitting it to the archive so that it may help someone in the future. plotfile <- "gant.eps" regular.color <- "lightgray" important.color <- "yellow" tasks <- list() #

Thanks to Ross Ihaka and Bob Wheeler for responding to my earlier question. I looked into the Johnson system functions in SuppDists package. For now, I want to stick with the density(...) estimator, and so still need the variate lookup function. As per Ross' suggestion, I just did a numerical integration on the density object and used approxfun/splinefun to "lookup" the variate

The handy function summary() doesn't work correctly with Date class objects: > R.version.string [1] "R version 1.9.1, 2004-06-21" > b <- as.Date(c("2002-12-26", "2002-12-27", "2002-12-28", "2002-12-29", "2002-12-30")) > b [1] "2002-12-26" "2002-12-27" "2002-12-28" "2002-12-29"

Hi folks: How do I get shell() or system() to run external programs in a particular directory (different from the current R working directory)? > version _ platform i386-pc-mingw32 arch x86 os Win32 system x86, Win32 status major 1 minor 3.0 year 2001 month 06

Dear All, I have a question concerning the computation of the value of the Restricted Maximum Likelihood (REML) function evaluated at a given set of parameter estimates from the Maximum likelihood (ML) value. Following the book of Fitzmaurice, Laird and Ware (2004) "Applied Longitudinal Analysis" pp101, the REML likelihood can be computed by multiplying the ML likleihood by the square

Hi folks: I used to load a DLL fine up until R 1.5.1. The same DLL dyn.loaded under R 1.6.2 prints the following warning message: > dyn.load("C:/Rajiv/Bin/rpdate.dll") Warning message: DLL attempted to change FPU control word from 8001f to 9001f > Any indication as to what the warning means, and whether I should be concerned about it? The DLL is made from a pure C code,

I'm making some 3D plots with contour3d from misc3d and wireframe from lattice. I want to view them from below; i.e. the negative z-axis. I can't figure out how to do so. I would like my point of view looking up from below, with the z, y, and x axes positive going away. Can anyone tell me the correct settings for screen to achieve this? Here is what I've found so far:

Hi, I am trying to use glmmPQL package for Generalized linear mixed models. This package works by repeated calls to lme. lme uses by default REML method for estimation. Then, does glmmmPQL use REML too? In contrast, how can I change it? I have tried it, writing : method="REML", but the program says: invalid method REML. If somebody can answer me....thanks, Sonja

Hi Greg, For quasi families I've used extended quasi-likelihood (see Mccullagh and Nelder, Generalized Linear Models 2nd ed, section 9.6) in place of the likelihood/quasi-likelihood in the expression for the (RE)ML score. I hadn't realised that this was possible before the paper was published. best, Simon ps. sorry for slow reply, the original message slipped through my filter for

Hi, I've been using gamm4 to build GAMMs for exploring environmental influences on genetic ancestry. Things have gone well and I have 2 very straightforward questions: 1. I've used method=REML. Am I correct that this is an alternative method for estimating the smooth functions in GAMMs rather than GCV that is often used for GAMs? I've read up on REML and it makes sense, but I'm

Hi, In glmmPQL in the MASS library, the function uses repeated calls to the function lme(), using ML. Does anyone know how you can change this to REML? I know that in lme(), the default is actually set to REML and you can also specify this as 'method=REML' or 'method'ML' but this isn't applicable to glmmPQL(). I'd appreciate any help or advice! Thanks, Emma

Hi, I know that is not possible make a stepwise procedure using REML in R, I can use ML for this. For nested design it may be very dangerous due the difference in variance structure, mainly in a splitplot design. ML make significative variables that REML dont make. I read an article that is made a stepwise procedure using GENSTAT. from article: "Terms were dropped from a model in a

Hi folks: This may be a stupid question, but I cannot seem to find a way to tell readBin that I want to read 2-byte integers from the connection. The input file is 150,720 bytes long containing 75,360 short (2-byte) integers. But specifying "integer" or "int" for what in readBin only returns me a vector of length 37680, leading me to believe that sizeof(integer) or