similar to: aov bug? (PR#1930)

I'm using R 1.7.0 on linux. With this version of R the package modreg is automatically loaded at start of session. However attempting to use predict.smooth.spline() produces Error: couldn't find function predict.smooth.spline. The function smooth.spline() is OK. What am I missing? ====================================== I.White ICAPB, University of Edinburgh Ashworth Laboratories, West

Apologies if this is this too obscure for R-help. In package splines, ns(x,,knots,intercept=TRUE) produces an n by K+2 matrix N, the values of K+2 basis functions for the natural splines with K (internal) knots, evaluated at x. It does this by first generating an n by K+4 matrix B of unconstrained splines, then postmultiplying B by H, a K+4 by K+2 representation of the nullspace of C (2 by K+4),

Can anyone explain why with latest version of R (2.0.1) on FC3, installed from R-2.0.1-0.fdr.2.fc3.i386.rpm, update.packages() produces the message /usr/lib/R/bin/Rcmd exec: INSTALL: not found. Indeed /usr/lib/R/bin seems to lack various shell scripts (INSTALL, REMOVE, etc). ====================================== I.White University of Edinburgh Ashworth Laboratories, West Mains Road Edinburgh

# Your mailer is set to "none" (default on Windows), # hence we cannot send the bug report directly from R. # Please copy the bug report (after finishing it) to # your favorite email program and send it to # # r-bugs@r-project.org # ###################################################### # The last line of following code produces a segmentation fault: x <- 1:10 f <- gl(5,2)

Hi, You should look at Pinheiro and Bates (2000) Mixed-effects models in S and S-Plus. It describes how to format the correlation matrix to pass to functions lme and gls. Basically, the correlation matrix has to be one of the corStruct classes, probably corSymm for your example. So in the call to lme (or gls if you really have no random effects), use something like:

Could anyone help with a linear mixed model fitting problem ? The model is : Y= Xp + Zu + e where X, Z are known design matrix, p is fixed effect factor, u is random effect, u~ (0, G) , e~(0,R) The main problem is , I want to fix the covariance matrix G to be a constant times a known covariance matrix A, G = c*A (c is positive constant, A is a predefined matrix with values manually set by

I am seeking for a method to calculate, given a fitted lme model and some data for a subject, the random effects predictors for this subject. I can only find predictors for the subjects used in creating the fit. Of course I could just add the subject and redo the fit. But I want to avoid just this refitting. Thanks for help wbk

After a one-way anova, se.contrast computes the standard error of a contrast, but not the value of the contrast itself. Wouldn't this be useful? Am I missing something? ====================================== I.White ICAPB, University of Edinburgh Ashworth Laboratories, West Mains Road Edinburgh EH9 3JT Fax: 0131 650 6564 Tel: 0131 650 5490 E-mail: iwhite at staffmail.ed.ac.uk

Hello everyone, I hope this is a fair enough question, but I don’t have access to a copy of Bates and Pinheiro. It is probably quite obvious but the answer might be of general interest. If I fit a fixed effect with an added quadratic term and then do it as an orthogonal polynomial using maximum likelihood I get the expected result- they have the same logLik.

Hi, I'm having trouble with glmmPQL. I'm fitting a 2 level random intercept model, with 90,000 cases and about 330 groups. I'm unable to get any results on the full data set. I can get it to work if I sample down to about 30,000 cases. But for models with N's much larger than that I get the following warning message:

Dear list, I use Sweave almost exclusively for writing papers, and I have become quite spoiled by the excellent xtable export facilities. Has anybody written an xtable method for the Anova function in CAR, or has anybody used a different set of functions to import Anova results into a table in an Sweave document? If not, any handy hints on how to write a good homebrew based on the output of Anova

Anybody have suggestions on having a 56K dialpool and VOIP connections with an Asterisk box over the same set of PRIs? We've done the PM3 with PRIs for just dialup, but are looking for a way to integrate our Asterisk box and move our voice calls onto the same PRIs. Ian -- Ian White Victoria Free-Net Association email: iwhite@victoria.tc.ca http://victoria.tc.ca/

On recent releases of the snom200 firmware, the MWI indicator will turn on, but won't turn off when the message has been checked. It works on firmware 2.03o, but not in 2.04g or newer. I filed a bug report with snom, but they're claiming it is an asterisk issue and that it should have been resolved. They suggested that I ask on the list. "Anyway, Asterisk had a bug where it

I can't get LaTeX to recognize the output from the pictex device driver. Are these commands for some special latex package which I don't know about? ************************************************ * I.White * * ICAPB, University of Edinburgh * * Ashworth Laboratories, West Mains Road * * Edinburgh EH9 3JT

I'm running R 1.2.2. The help information for manova says that the result is "A list with components SS: A names list of sums of squares and product matrices. Eigenvalues: A matrix of eigenvalues, stats: A matrix of the statistics, approximate F value and degrees of freedom." However, when I run the example, with fit <- manova(Y ~ rate*additive), I find that fit$SS is NULL.

Hi, I would like to repsond to a few questions raised in a reply to a question posted by another user on March 17, 2004. The entire message is copied at the end of this email. The relevant questions and statements are as follows: What did you not understand about help(memory.size)? This is also in the rw-FAQ: what in that did you not understand? ... Yes, so try a machine with 2Gb RAM. I

BUG IN SPLINE()? Version R-1.0.1, system i486,linux If the spline(x,y,method="natural") function is given values outside the range of the data, it does not give a warning. Moreover, the extrapolated value reported is not the ordinate of the natural spline defined by (x,y). Example. Let x <- c(2,5,8,10) and y <- c(1.2266,-1.7606,-0.5051,1.0390). Then interpolate/extrapolate with

#I am trying to understand how R fits models for contrasts in a #simple one-way anova. This is an example, I am not stupid enough to want #to simultaneously apply all of these contrasts to real data. With a few #exceptions, the tests that I would compute by hand (or by other software) #will give the same t or F statistics. It is the contrast estimates that R produces #that I can't seem to

https://bugzilla.mindrot.org/show_bug.cgi?id=1930 Damien Miller <djm at mindrot.org> changed: What |Removed |Added ---------------------------------------------------------------------------- Status|NEW |RESOLVED Resolution|--- |FIXED --- Comment #1 from Damien Miller <djm at

CentOS Errata and Security Advisory 2015:1930 Important Upstream details at : https://rhn.redhat.com/errata/RHSA-2015-1930.html The following updated files have been uploaded and are currently syncing to the mirrors: ( sha256sum Filename ) i386: 3cd9b3811298e411dfa65eb67a74b61bdee315640dfc11d1f02c28c8953af6aa ntp-4.2.6p5-5.el6.centos.2.i686.rpm