similar to: Mixed Model nested ANOVA (lme help)

Hello again R help. This is a simple question perhaps(laughing as I type) with a simple answer. Is the multcomp package appropriate for using with a lme built under the nlme package. I know i can get it to work (i.e report p-values for a tukey test) but i am not unsure if this is appropriate for a linear mixed effects model. Below is the code i used lmm1 <- lme(asinh(density_recruit) ~

The multcomp package has code in multcomp:::expression2coef that attaches the 'coef' attribute to symbols. Since there is only one symbol object in a session with a given name, this means that this attaching has a global effect. Should this be quietly allowed or should there be a warning or an error? E.g., str(quote(Education)) # symbol Education lmod <- stats::lm(Fertility ~ .,

Hello, is it possible to obtain type III sums of squares for a nested model as in the following: lmod <- lm(resp ~ A * B + (C %in% A), mydata)) I have tried library(car) Anova(lmod, type="III") but this gives me an error (and I also understand from the documentation of Anova as well as from a previous request (http://finzi.psych.upenn.edu/R/Rhelp02a/archive/64477.html) that it is

Here is a simpler example: > ex <- as.name("a") > attr(ex, "test") <- 1 > quote(a) a attr(,"test") [1] 1 Torsten On Thu, 6 Jul 2017, William Dunlap wrote: > The multcomp package has code in multcomp:::expression2coef that attaches the 'coef' attribute to > symbols.? Since there is only one symbol object in a session with a given name,

Thanks for spotting this issue. The short answer is yes, adding attributes to a symbol is a bad idea and will be turned into a runtime error soon. Maintainers of packages that add attributes to symbols have been notified and some have already fixed their code. At least in one case the package is not working properly, even in isolation, because of the global effect of adding an attribute to a

Hey, I have a data set of 1.20MB. I used read.table( ) to read in the data and then tried to convert it to geodata type. But it failed and I got "Error in vector("double", length) : cannot allocate vector of length 1351974000" message. All the commands I used were: rawdata<-read.table(file=" ") datag<-as.geodata(rawdata) BTW, what is the biggest origin

Hey there! I got a pseudo code and don't know how to apply it to R, maybe someone can help me: Input: A dataset X, kmax: maximum number of clusters, num_subsamples: number of subsamples. Output: S(i; k) - a distribution of similarities between partitions into k clusters of a reference clustering and clustering of subsamples; i = 1 to num_subsamples Requires: T = cluster(X): A hierarchical

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 The lmer() function in the lme4 package has some code of the form mc <- match.call() mc[[1]] <- as.name("lFormula") lmod <- eval(mc, parent.frame(1L)) this is a fairly common idiom in R, found e.g. in lm(), used when one wants to pass all of the arguments of a function to a different function (in the case of lm() it's

I tried to use calcMin with a function that uses a number of ... arguments (all args from resid on) besides the vector of parameters being fit. Same idea as optim, nlm, nlminb for which this form of ... syntax works. But with calcMin I get an error regarding unused arguments. No partial matches to previous arguments that I can see. Anybody know the reason or fix for this?

Full_Name: Tom Short Version: 2.4.0 OS: Windows XP Submission from: (NULL) (68.236.159.227) It looks like there's a bug in acosh with complex number in windows: > acosh(2) [1] 1.316958 > acosh(2+0i) [1] 0+NaNi This happens for me on Windows XP with the following versions: R version 2.2.0, 2005-10-06, i386-pc-mingw32 R version 2.4.0 (2006-10-03) i386-pc-mingw32 It works fine with

Hi, I'm currently in the process of writing an R-installation SOP for my company. As part of that process I'm using the recommendations from the 'R Installation and Administration' document, section 3.2, "Testing an installation". This is done on an XP machine, using the latest binary of 2.11.0. The binary is downloaded and then installed from the installer. I then

I am very new to R and programming and thank you in advance for your patience and help with a complete novice! I am working with a large multivariate data set that has 10 explanatory environmental variables (e.g. temp, depth) and over 60 response variables (each is a separate species). My data frame is set up like the simplified version below: JulianDay Temperature Salinity Depth Copepod

Full_Name: Brendan McMahon Version: 1.3 OS: SunOS Submission from: (NULL) (12.42.50.51) Probably just a configuration problem but ... fails to compile. configure completed without a problem. Any help appreciated!! SunOS flexdev1 5.7 Generic_106541-11 sun4u sparc SUNW,Ultra-5_10 `Makedeps' is up to date. gcc -I. -I../../src/include -I../../src/include -I/usr/local/include -DHAVE_CONFIG_H

Hi I am having difficulty with setMethod(). I have a "brob" class of objects whose representation has two slots: "x" and "positive". Slot "x" (double) holds the log of a number and slot "positive" (logical) its sign. The idea is that large numbers can be handled. I'm trying to implement a log() method using an analogue of the

RSQLite.extfuns provides SQLite extension functions for use with RSQLite. The package is a wrapper of extension functions written by Liam Healy and made available through the SQLite website (http://www.sqlite.org/contrib). You can make the extension functions available on a per db connection basis like this: library("RSQLite.extfuns") db <- dbConnect(SQLite(), dbname =

RSQLite.extfuns provides SQLite extension functions for use with RSQLite. The package is a wrapper of extension functions written by Liam Healy and made available through the SQLite website (http://www.sqlite.org/contrib). You can make the extension functions available on a per db connection basis like this: library("RSQLite.extfuns") db <- dbConnect(SQLite(), dbname =

Hi, I am running R 2.13.0 on a Windows 7 machine. I ran the script: testInstalledBasic('devel') and received the following warning message: running tests of consistency of as/is.* creating ?isas-tests.R? running code in ?isas-tests.R? comparing ?isas-tests.Rout? to ?isas-tests.Rout.save? ...running tests of random deviate generation -- fails occasionally running code in

I haven't been able to build yesterday's R-devel (or, as a test, R-devel_2005-03-11) on Windows XP. The error message I'm getting is to do with 'arithmetic.c', as shown below. Current R-patched seems to build fine, and I've built previous R-devels from source fairly often in the past-- though I have no understanding whatsoever of what's going on. FWIW the only changes

I am trying to install R-1.4.1 on the Solaris 2.7 OS by compiling the source code. Unfortunately, I obtain a fatal error during make (apparently while using gcc to compile arithmetic.c). The configure command seemed to work acceptably (end of output follows). --(snip)-- R is now configured for sparc-sun-solaris2.7 Source directory: . Installation directory: /usr/local C

Hi there, I'm using the lattice package to create an xy plot of abundance vs. depth for 5 stages of barnacle larvae from 5 species. Each panel of the plot represents a different stage, while different loess smoothers within each panel should represent different species. However, I would like depth to be on the y-axis and abundance to be on the x-axis, because this is more intuitive as an