similar to: xyplot link missing

Dear List, This may be related to this email thread initiated by Ben Bolker last month: https://stat.ethz.ch/pipermail/r-devel/2011-July/061630.html In answering this Question on StackOverflow http://stackoverflow.com/q/7195628/429846 I noticed that `methods()` was not listing some S3 methods for `plot()` provided by the mgcv package. At the time I wanted to check the development version of R as

Hi all, I figured out some things about an earlier question I posted today; this is a new question: I have a for-loop in which I have the following lines for (i in list.files(pattern = "*.dmp")) { currentfile <- read.table(i) vecrawrts <- as.vector(currentfile$V3) ... } Now, vecrawrts is numeric, but if I run the script on all the files with extension .dmp, I get the

Hello everybody, I'm trying to subset a data frame "foo", which looks like this: sent cond raw V1 10 2 A 614 7 11 2 A 724 6 12 2 A 641 5 13 2 A 665 9 33 5 B 510 7 34 5 B 659 7 35 5 B 607 6 37 5 B 467 8 38 5 B 586 9 39 5 B 523 10 57 8 C 858 8 58 8 C 513 9 59 8 C

The Rows function which is called from plot.compareFits in the nlme package is not found. > plot(compareFits(coef(bp.model3),coef(bp.model3M))) Error in plot.compareFits(compareFits(coef(bp.model3), coef(bp.model3M))) : couldn't find function "Rows" > Can I find it elswhere? Have I missed a required package? Thanks Ross Darnell > library(help=nlme) nlme

Dear all, I'm trying to install lme4 (after having installed R 2.4.0 from source, and having installed the latest Matrix package). lme4 fails with the following message: pedigree.o definition of _lme4_xSym in section (__DATA,__common) pedigree.o definition of _lme4_ySym in section (__DATA,__common) make: *** [lme4.so] Error 1 ERROR: compilation failed for package 'lme4' ** Removing

Hi all, I'm building a plot of the values in tmeant (below) against positions 1 to 5, using matplot. tmeant looks like this: case1 case2 pos1 861.8466 818.5909 pos2 961.2841 976.3466 pos3 878.6080 1262.8523 pos4 950.8011 1129.6080 pos5 968.1080 1063.3920 I also have lower (object tl) and upper (object tu) bounds on the confidence intervals as follows: tl: pos1

Hi folks, I want to have a y-axis using matplot where the y axis is truncated as shown schematically below: (msecs) 3000 | x ~ ~ 800 | x | | x | ----|----|----|----------------- 1 2 3 Position The ~ is supposed to stand for a curved line showing a break. The issue is that I want to plot the points shown by the

Hi all, All the files related to the question below are in http://ling.ohio-state.edu/~vasishth/TempRdir/ I've been using a script (called newresiduals.R, the code is in the above directory; I wanted to avoid cluttering up this message) that is run in the following manner at the command prompt: R --vanilla "wordlen-2-01-g1" < newresiduals.R What it does is: it goes through a

Greetings, I have a groupedData (nmGroupedData) object created with the following syntax: Soil <- groupedData( ksat ~ conc | soil_id/sar/rep, data=soil.data, labels=list(x='Solution Concentration', y='Saturated Hydraulic Conductivity'), units=list(x='(cmol_c)', y='(cm/s)') ) the original data represents longitudinal observations in the form of:

Hi all, I have the following dataset, call it test.data (30 columns, and one row named "0"): ADVP ADVP AP AP CONJ CONJ CP CP DU DU INF INF MWU MWU NP NP PP PP PPRT PPRT REL REL SMN SMN SSB SSB SV1 SV1 TI TI 0 96.85 2.05 89.07 2.54 70.91 2.37 94.92 3.46 82.31 11.33 40.96 2.25 98.06 3.43 90.77 17.63 86.60 10.78 60.27 1.32 93.27 0.97 77.60

Hi all, I'm trying to figure out the example about mixed models in the Pinheiro and Bates book (Mixed-Effects Models in S and S-Plus, 2000, pp. 135-137). One thing I don't understand is: When I run the command fm1Orth.lm <- lm( distance ~ age, Orthodont ) followed by fm2Orth.lm <- update( fm1Orth.lm, formula = distance ~ Sex*age ) and then do summary(fm2Orth.lm)

Dear R users, Many of you are aware of the statistical package lme4 used in R for carrying out mixed-effects modeling. The lme4 R-Forge website [ https://r-forge.r-project.org/projects/lme4/] is going to maintain a list of all dissertations and articles (accepted for publication, in press, or published) that use lme4 or its precursor, nlme. It would be helpful if members of this list could

I'm trying to write a function that has a dataframe as an argument. Within the function, I call groupedData to create a grouped-data object from the dataframe argument. However, the groupedData function doesn't seem to see the data. I'm guessing it's getting the argument name rather than its value, but I'm not positive. My system: R 1.3.1 on WinNT 4.0 (Pentium II), nlme

Hi, My apologies for sending this message earlier to the wrong list (next time I'll wait till I'm fully awake). I'm a new R user. I want to do the following: I want to pass a filename to R repeatedly (about 400 times) in order to do a repetitive task (calculating regressions). The way I'm setting it up is to write a bash script that will create an R source file on the fly for

Dear all, I'm using the nlme package. I would like to introduce the plot of my groupedData object in my latex report using Sweave but it doesn't work. I have no error message in R (Sweave(my.file_Rnw)) and the plot is generated in R, but when I want to see the pdf (I'm using TexnicCenter) of my report, there is an error message and no pdf file is generated. It seems that the plot

Thanks for the advice, folks! Shravan Vasishth <vasishth at ling.ohio-state.edu> Try pico. There's no substitute for/comparison with emacs, though, IMHO. but just to be sure: you don't have to have ESS to use (x)emacs for editing R scripts. Sven Garbade <garbade at psy.uni-muenchen.de> ESS can be installed via rpm on Mandrake. baron at cattell.psych.upenn.edu (Jonathan

Dea-R community. I'd like to draw your attention to an issue I have recently encountered while doing my current data analysis. I've got an unexpected (to me) result from the command: > augPred(lmList(my.object)), 'my.object' being a grouped data frame of class: > class(my.object) [1] "nfnGroupedData" "nfGroupedData" "groupedData"

Hi all, This question is about using the linear-mixed effects model for a within-subjects experiment. The experiment description: Each subject was asked to read a series of sentences (in Hindi) on a computer screen, and pressed the space bar to see each word in the sentence successively. The amount of time spent on each word was recorded on the computer, in milliseconds; this is the dependent

Dear All, I wish to model random effects that have known between-group covariance structure using the lme() function from library nlme. However, I have yet to get even a simple example to work. No doubt this is because I am confusing my syntax, but I would appreciate any guidance as to how. I have studied Pinheiro & Bates carefully (though it's always possible I've missed

Dear All, I would be grateful if you can tell how can I suppress the legend (automatically created) when I plot a groupedData. suppose that I have data farme df_0 which contains the following columns: ID represents the clusters, t represents time of observation, Y: the variable of interest and X a binary covariate. Then: df_1 <- groupedData(Y ~ t | ID) plot(df_1,