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

Hello, R users, I have two factors (treat, section) anova design experiment where there are 3 replicates. The objective of the experiment is to test if there is significant difference of yield between top (section 9 to 11) and bottom (section 9 to 11) of the fruit tree under treatment. I found that there are interaction between two factors. I wonder if I can contrast means from levels of

Dear Thomas Lumley, and R-help list members, I have read your article "Network meta-analysis for indirect treatment comparisons" (Statist Med, 2002) with great interest. I found it very helpful that you included the R code to replicate your analysis; however, I have had a problem replicating your example and wondered if you are able to give me a hint. When I use the code from the

Hello, R experts, I have a list called dataHP which has 30 elements (m1, m2, ..., m30). Each element is a 7x6 matrix holding yield data from two factors experimental design, with treatment in column, position in row. For instance, the element 20 is: dataHP[[20]] col1 col2 col3 trt1 trt2 trt3 [1,] 22.0 20.3 29.7 63.3 78.5 76.4 [2,]

I'm running an ANOVA on some data for respiration in a forest. I am having a problem with my degrees of freedom. For one of my variables I get one fewer degrees of freedom than I should. I have 12 plots and I therefore expected 11 degrees of freedom, but instead I got 10. Any ideas? I have some code and output below: > class(Combined.Plot) [1] "character" >

Hello, R experts, If I understand Ted's anwser correctly, then I can not contrast the mean yields between sections 1-8 and 9-11 under "Trt" but I can contrast mean yields for sections 1-3 and 6-11 because there exists significant interaction between two factors (Trt:section4, Trt:section5). Could I use the commands below to test the difference between sections 1-3 and 6-11 ?

Dear R-list, I have a partly comma separated partly underscore separated string that I am trying to parse into R. Furthermore I have a bunch of them, and they are quite long. I have now spent most of my Sunday trying to figure this out and thought I would try the list to see if someone here would be able to get me started. My data structure looks like this, (in a example.txt file) Subject

Dear, I have a dataset with 4 subjects (see ID in example), and 4 treatment (see TRT in example) which are tested on 2 locations and in 3 blocs. By using Lattice dotplot, I made a graph that shows the raw data per location and per bloc. In that graph, I would like to have a reference line per bloc that refers to the first treatment (T1). However, I can not find how to do that. I can make