similar to: polynomials REML and ML in nlme

I have data that looks like Year,Recruit,Spawner,Mtempcv 1958,4532,775,0.24125 1959,22996,2310,0.16319 1960,628,2990,0.46056 1961,879,1400,0.33028 1962,14747,1130,0.22618 1963,13205,790,0.20596 1964,31793,1195,0.19229 1965,10621,981,0.20363 1966,22271,870,0.3452 1967,8736,1104,0.27511 1968,8761,883,0.10884 1969,18885,1421,0.17799 1970,10098,1198,0.2106 1971,3394,760,0.22098 1972,1697,1354,0.39461

Hello Everyone - I want to print a number of results from lme function objects out to a txt file. How could I do this more efficiently than what you see here: out2 <- capture.output(summary(mod2a)) out3 <- capture.output(summary(mod3)) out4 <- capture.output(summary(mod5)) out5 <- capture.output(summary(mod6)) out6 <- capture.output(summary(mod7))

1. I'm attempting to test for Random Effects. I've grouped the data on subject (grid) but want to use lme to build the model without subject as a RE then add it and do anova between the 2 models. This is the result I get and it appears it's adding Random Effects. tmp.dat4 <- groupedData(Trials ~ 1 | grid, data = tmp.dat4) mod2a <- lme(Trials ~ factor(group_id) + reversal,

Hi R-Helpers: I?m trying to use the function GLMM from lme4 package, (R-1.8.1, Windows 98),and I get the following error: > pd5 = GLMM(nplant~sitio+ + fert+ + remo+ + sitio:fert+ + remo:sitio+ + remo:fert+ + remo:fert:sitio + data=datos, + family=binomial, +

Could anyone help with what should be a simple task? I have data as a fixed format (fortran) table. I have no trouble getting it into R using read.table. Each column is separated by a space, including the first column that begins with a space, and aligned. It reads into R as if separated by tabs. However I want to manipulate two columns of data then write the results out into exactly the

Hello, My question is the following. Has anyone put together an "idiots" checklist of the steps needed to build a (small personal) package under Windows, or documented their own experiences of building packages under Windows? I ask this as last week there were some very interesting comments regarding the advantages of putting personally written functions together as a package

Does anyone have a way of producing solid shapes other than spheres in rgl? I am using rgl to produce a simple visualisation of a forest model results using "lollipops". Its just a bit of fun, but as many of the trees are pines I would like to depict their crowns as cones. If there is a solution I need it to work under windows. Here is the example. library(rgl) library(misc3d)

I have found that the "schema.table" syntax used in Postgresql (and Oracle) does not work directly with RODBC. This works library(RODBC) con<-odbcConnect("mydb") d<-sqlQuery(con,"select * from meso.trees") However this does not. d<-sqlFetch(con,"meso.trees") Error in odbcTableExists(channel, sqtable) : ?meso.trees?: table not found on channel

Hello, I've seen this question asked in the archives but no clear reply or solution provided. So, just to be sure it is not possible in R: Can I replace backslashes with slashes in a string ? I am writing a GUI for R with the Rpad library. I have a "browse" button for data loading and Windows return a path string with backslashes. I need to convert them into slashes to use the

I am running R on Ubuntu Heron. I recently updated to R 2.8.0 from the default (2.6) in the Ubuntu heron repository by adding deb http://cran.stat.sfu.ca/bin/linux/ubuntu hardy/ I updated all packages installed in /usr/lib/R/site-library and those in my home directory /home/duncan/R/i486-pc-linux-gnu-library/2.6 without obvious issues. I installed adehabitat and gpclib, which should be good

Dear List-members, Hopefully someone will help through my confusion: In order to get the same coefficients as we get from the following ## require (MASS) summary ( lm(Gas ~ Insul/Temp - 1, data = whiteside) ) ...................... we need to do the following (if we use model.matrix to specify the model) ## summary ( lm(Gas ~ model.matrix(~ Insul/Temp - 1) - 1, data = whiteside) )