similar to: Contrasts for 2x4 interaction in mixed effects model

Hello everyone, I am testing a model in which I have a two-level factor (let's call it First [1, 2]) nested under a four-level factor (let's call it Second [A, B, C, D]). I have used the following model to get coefficients representing whether, for each level of Second, there is a significant difference (in the outcome variable, Latency) between the levels of First: test <- lmer(

Hello everyone, Does anyone know what would cause the skewness() function (from e1071), as well as skew() from psych, to return a value of NaN? I have a vector of positively-skewed data (https://docs.google.com/file/d/0B6-m45Jvl3ZmYzlHRVRHRURzbVk/edit?usp=sharing) which these functions return a value for like normal: > skewness( data ) # returns 1.400405 but when I instead give those

Right, that's my point. The help page mentions a `title.cex`, like I said; saying that `cex` sets the default `title.cex` sure implies to me (and presumably to the other people whose discussion I linked) that a `title.cex` parameter exists. Since no such parameter exists, this bit in the documentation is misleading (suggesting that there is a `title.cex` parameter which can be set, when there

Hello everyone, I've conducted a Type III repeated-measures ANOVA using Anova() from the car package, based on the suggestions at http://blog.gribblelab.org/2009/03/09/repeated-measures-anova-using-r/(option 3) and http://languagescience.umd.edu/wiki/EEG#ERP_ANOVA_in_R. My ANOVA has two factors: Condition (3 levels) and Region (6 levels) and their interaction. Below is code to run the Anova

To whom it may concern: The help page for ?legend refers to a `title.cex` parameter, which suggests that the function has such a parameter. As far as I can tell, though, it doesn't; here's an example: > plot(1,1) > legend("topright",pch=1, legend="something", title="my legend", title.cex=2) Error in legend("topright", pch = 1, legend =

Hello, duplicated() does not seem to work for a long vector. For example, if you download the data from https://docs.google.com/open?id=0B6-m45Jvl3ZmNmpaSlJWMXo5bmc (a vector with about 12,000 numbers) and then run the following code which does duplicated() over the whole vector but just shows the last 30 elements: data.frame( tail(verylong, 30), tail(duplicated(verylong), 30) ) you'll see

Hello, I have a list of data with multiple elements, and each element in the list has multiple variables in it. Here's an example: ### Make the fake data dv <- c(1,3,4,2,2,3,2,5,6,3,4,4,3,5,6) subject <- factor(c("s1","s1","s1","s2","s2","s2","s3","s3","s3",

>>>>> POLITZER-AHLES, Stephen [CBS] <stephen.politzerahles at polyu.edu.hk> >>>>> on Sat, 25 Mar 2017 13:25:32 +0000 writes: > Right, that's my point. The help page mentions a > `title.cex`, like I said; saying that `cex` sets the > default `title.cex` sure implies to me (and presumably to > the other people whose discussion

> On 25 Mar 2017, at 00:39 , POLITZER-AHLES, Stephen [CBS] <stephen.politzerahles at polyu.edu.hk> wrote: > > To whom it may concern: > > > The help page for ?legend refers to a `title.cex` parameter, which suggests that the function has such a parameter. No it does not. All arguments are listed and documented, none of them is title.cex, and there's no

>>>>> Jan van der Laan <rhelp at eoos.dds.nl> >>>>> on Fri, 6 Oct 2017 12:13:39 +0200 writes: > It is actually model.matrix that crashes, not glm. Same > crash occurs with e.g. lm. > model.matrix(dob_mon ~ dob_day*dob_mon, data = tab) > also crashes R. Yes, segmentation fault. It only happens when these are *logical*

Dear all, I studied in tank prey fish behavior. Using the design described below (and R code), I want to test the effects of both habitat and predator (and interaction) on prey fish's vertical distribution, which was recorded (with repeated measures) as a categorical variable. I found that package mlogit might fit to my need but I don't know how to specify my complex design in the

The following code crashes R (I know I shouldn't try to estimate such a model; this was a bug in some code of mine). I also tried with R-devel; same result. tab <- structure(list(dob_day = c(FALSE, FALSE, FALSE, FALSE, TRUE, TRUE, TRUE, TRUE), dob_mon = c(FALSE, FALSE, TRUE, TRUE, FALSE, FALSE, TRUE, TRUE), dob_year = c(FALSE, TRUE, FALSE, TRUE, FALSE, TRUE, FALSE, TRUE), n =

Hello Amelie, I don't have an answer to your question, but I just wanted to point out this page I noticed recently ( http://hlplab.wordpress.com/2009/05/07/multinomial-random-effects-models-in-r/), which might be helpful. I'm also interested in figuring out how to do a multinomial glmm, so if you find out anything I'd be happy to hear more about it! Based on what I've found so

hi, folks: i need to transpose the following data: gene tissue patient1 patient2 patient3..... --------------------------------------------- gene1 breast 10 100 1 gene2 breast 20 200 4 gene3 breast 30 50 5 gene4 breast 40 400 9 ................................ to the

Dear all, I have following data called "paitent" day patient1 patient4 patient5 patient6 0 -0.27842688 -0.04080808 -0.41948398 -0.04508318 56 -0.22275425 -0.01767067 -0.30977249 -0.03168185 112 -0.08217659 -0.26209243 -0.29141451 -0.09876170 252 0.08044537 -0.26701769 0.05727087 -0.09663701 where each patient have response values at four time points. I want to

I am making some changes to the permax library (so that it will accept NA's). This function performs a permutation analysis to identify discriminating attributes distinguishing two groups of observations. It takes the form (at its most simplistic): permax(data, ig1) where ig1 is one group of interest. The other group (if not specified) is assumed to be the remaining observations, namely,

hi, everyone: i have a question on the reshape function. i have the following dataset : gene tissue patient1 patient2 patient3............. _________________________________________________ gene1 breast 10 20 50 gene2 breast 20 40 60 gene3 breast 100 200 300 which i hope to convert to the following format: gene patientID

Dear R users, I have a reproducible data and try to create new variable "clo" is 1 if know variable is equal to "very well" or "fairly well" and getalong is 4 or 5 otherwise it is 0. rep_data<- read.table(header=TRUE, text=" id1 id2 know getalong 100000016_a1 100000016_a2 very well 4 100000035_a1 100000035_a2 fairly

Dear r-list, I've got a data base: > HData[1:10,] NumTree Site Species Date Age DBH H IdentTree 1 1 Queige Spruce 2002 184 49 33.5 Queige 1 2 2 Queige Fir 2002 NA 5 4.6 Queige 2 3 3 Queige Fir 2002 25 8 6.6 Queige 3 4 4 Queige Spruce 2002 198 47 32.5 Queige 4 5 5 Queige Fir 2002 200 59 35.3 Queige 5 6 6 Queige

Hi Helpers: I have a factor object that has 314k entries of 39 land cover types. (This object can be coerced to characters neatly should that be easier to work with.) > length(foo) [1] 314482 > foo[1:10] [1] Montane Chaparral Barren Red Fir Red Fir [5] Red Fir Red Fir Red Fir Red Fir [9] Red Fir Red Fir 39 Levels: