similar to: data files in R

Hi all, I would be grateful if anyone could help me with the following. I am using nlme library and I am trying to do a trellis plot with an outer factor, but I have an error message which I can't understand. Here is the code : > mydata <- groupedData(y ~ x | warren/rabbit, outer= ~ treatment, data=mydata) > plot(mydata) # I obtain a plot with all rabbits displayed individually and

Hi all, I would like to know how the BIC criterion is calculated for models estimated using gls( ) function. I read in Pinheiro & Bates (2000) p84 that BIC = -2logL + npar*log(N) (for the ML method), or BIC = -2logLR + npar*log(N-p) (for the REML method) but when I use any of these formulae I don't obtain the result given by R. Thanks in advance for any help. Eve CORDA Office national

Dear all, I wonder if it is possible to obtain marginal tests for effects in generalized linear models. Indeed, the anova function produces sequential tests and it doesn't have any "type" argument to specify that we would like marginal tests instead, as in the similar anova function for lme objects. Thanks a lot for your help! Eve CORDA Office national de la chasse et de la faune

A non-text attachment was scrubbed... Name: not available Type: text/enriched Size: 1284 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/20020911/e68e2ecc/attachment.bin

I'm looking at Pinheiro & Bates "Mixed-Effects Models in S and S-PLUS" at the moment. Several datasets are used, one of which is called "PBIB" (a partially balanced incomplete block design). All the other datasets can be found somewhere or other in R. However, I cannot locate PBIB, and it does not seem to be mentioned in the latest edition of the R Full Reference

The bio.infer package contains a data frame /usr/lib/R/library/bio.infer/data/itis.ttable.rda that needs to be modified. After loading the bio.infer package and attaching the data frame with the data() function, I wrote the data frame to a text file. After adding another row to the data frame I applied read.table() to create a data frame, but it's in my pwd, not the R library data

Hello, Is there an easy way to format the output of a result table that R generates from a regression? I like the table, but would like to limit the number of decimal points in the entries if possible. For instance I would like only 3 digits of precision for the Value, Std.Error. (And if it would be easy to get rid of scientific notation, that would be good to know too). So ideally keep the

Dear all, Lets say I have the following: > x <- c("Eve: Going to try something new today...", "Adam: Hey @Eve, how are you finding R? #rstats", "Eve: @Adam, It's awesome, so much better at statistics that #Excel ever was! @Cain & @Able disagree though :(", "Adam: @Eve I'm sure they'll sort it out :)", "blahblah") > x [1]

I have csv data that extend beyond the area I want for an existing map. I want using the boundaries of the polygon shape file as a cookie cutter so that I can overlay the csv data on map without including anything outside the map boundaries and create a dbf file or shapefile of the clipped data .  The reproducible example: ############################################### library(RColorBrewer)

Hey all I'm using Wine 1.1.22 and I'm having a hell of a time getting Eve-Online to work. I downloaded eve, installed it into Wine, everything seem to go well, no problems. However I when I run it it shows the prescreen and then seems to close itself and just shows the Wine Desktop. When I run it in form the Command Shell this is what I get: darkslayer at

I?d like to report the ZFS related crash/bug described below. How do I go about reporting the crash and what additional information is needed? I?m using my own very simple test app that creates numerous directories and files of randomly generated data. I have run the test app on two machines, both 64 bit. OpenSolaris crashes a few minutes after starting my test app. The crash has occurred on

A picture is worth a thousand words: http://dl.eve-files.com/media/corp/aurix/eve_premium_wine_0.9.61_1.png The black columns are supposed to be spot lights on the ship in the middle. I haven't tried eve premium before 0.9.60, however, people claim that it didn't have this problem in 0.9.58. Here is a link to the thread in eve online:

Hi all, I've an asterisk 1.2.5 running very well for about a 9 months, and suddenly i cannot dial extensions 4XXX from SIP Phones. Now comes the wired stuff... I can dial this extensions from IAX phones as well as from Analogue extensions connected to our legacy pbx, that is installed on front of asterisk. So : Zapata Calls to SIP extensions 4XXX - OK IAX to SIP 4XXX-OK SIP to SIP 4XXX -

R-helpers: I would like to measure the correlation coefficient between the repeated measures of a single variable that is measured over time and is unbalanced. As an example, consider the Orthodont dataset from package nlme, where the model is: fit <- lmer(distance ~ age + (1 | Subject), data=Orthodont) I would like to measure the correlation b/t the variable "distance" at

Consider the following plot that shows separate regression lines ~ age for each subject in the Pothoff-Roy Orthodont data, with separate panels by Sex: library(nlme) #plot(Orthodont) xyplot(distance ~ age|Sex, data=Orthodont, type='r', groups=Subject, col=gray(.50), main="Individual linear regressions ~ age") I'd like to also show in each panel the pooled OLS

I tried to build the gold plugin and receive the errors posted below. I checked out the gold plugin using "cvs -z 9 -d :pserver:anoncvs at sourceware.org:/cvs/src co src," as documented on llvm.org. Is that version bleeding edge and perhaps not stable? The first errors have to do with the libtool version. I have libtool 2.2.6 on my system, which is what the README-maintainer-mode

Hi everybody! I'm sure that I overlook something and feel quite stupid to ask, but I have not found an easy solution to the following problem: Take e.g. the Orthodont data from the nlme package: > head(Orthodont) Grouped Data: distance ~ age | Subject distance age Subject Sex 1 26.0 8 M01 Male 2 25.0 10 M01 Male 3 29.0 12 M01 Male 4 31.0 14 M01 Male

I know these topics exist, I have tried a multitude of them, but no show. My eve stops each time at the eula. screen is black after the splash screen until the eula tries to visit my screen. (for those not in the know, I just need to click yes to actually play the game) when running it through terminal, I get the following (... equals a lot of the line above it) bram at

I am using rw0641. In my continuing quest to understand repeated measures analysis, I again return to lme. I exported the Potthoff and Roy data Orthodont.dat from S-PLUS 4.5 to avoid capture errors and ran the examples in the R help. I imported the data.frame with data <- read.table("Orthodont.dat",header=T) attach(data) and created the objects Orthodont.fit1 <-

In a paper I'm writing using Sweave, I make use of lattice graphics, but don't want to explicitly show (or explain) in the article text the print() wrapper I need in code chunks for the graphs to appear. I can solve this by including each chunk twice, with different options, as in <<ortho-xyplot1-code, keep.source=TRUE, eval=FALSE>>= library(nlme) library(lattice)