similar to: anova grouping of factors in lme4 / lmer

Hello, I’m trying to fit the following model: Dependent variable: MAXDEPTH (the maximum depth reached by a penguin during a given dive) Fixed effects: SUCCESSMN (an index of the “individual quality” of a bird), STUDYDAY (the day of the study, from -5 to 20, with 0=Dec 20), and the interaction SUCCESSMN*STUDYDAY Random effect: BIRD (the bird id, as each bird is performing several dives)

Dear All, I want to test whether a method exists for given object. For example, whether a function "deviance" is defined for an object of the "lm" class. My imperfect understanding leads me to think something like hasMethod("deviance", object) hasMethod("deviance", "lm") existsMethod("deviance", signature(class="lm"))

Dear All, I have been trying to simulate data from a fitted glm using the simulate() function (version details at the bottom). This works for lm() fits and even for lmer() fits (in lme4). However, for glm() fits its output does not make sense to me -- am I missing something or is this a bug? Consider the following count data, modelled as gaussian, poisson and binomial responses: counts

Dear All, Is there is a fast way of refitting lm() when the design matrix stays constant but the response is different? For example, y1 ~ X y2 ~ X y3 ~ X ...etc. where y1 is the 1st instance of the response vector. Calling lm() every time seems rather wasteful since the QR-decomposition of X needs to be calculated only once. It would be nice if qr() was called only once and then the same

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

Christoph -- There was just a thread on this earlier this week. You can search in the archives for the title: "refitting lm() with same x, different y". (Actually, it doesn't turn up in the R site search yet, at least for me. But if you just go to the archive of recent messages, available through CRAN, you can search on refitting and find it. The original post was from William

Dear R help, I've encountered a problem with survreg's anova(). I am currently writing general code to fit a variety of models using different fitting functions. Here's a simple example of what I'm trying to do: ---begin code--- # general function to analyse data analyse.data <- function(formula, FUN, data, ...) { fit <- FUN(formula, data=data, ...) anova(fit)

Hi. I am not that good at R but I was wondering if there is either a tool or a strategy for testing many different models in R in a batch. I have used something in Weka called the Experimenter interface which helps with doing this kind of thing. Thank you. Andy [[alternative HTML version deleted]]

Hi, guys: These days I keep using R and Weka to do data mining. I think my next step is open the source codes so that I can "customrize them" and make them better server my purpose. But now I kinda hesitate to do so b/c I am really not sure which is better to start with. You know, both require some time and I cannot clone myself to work on both:) If here are some persons who used both

Hi everyone, My purpose is to read a .CEL file into R. The .CEL file was created from a .CAB by using DTT software found on Affymetrix website I read the .CEL file in R using ReadAffy as follows: > d2=ReadAffy(widget=T) and I complete the fields as required. It does not complain. For example I could find the description: > description(d2) Experimenter name: BB Laboratory: FFL Contact

I'm working with a genomic data-set with ~31k end-points and have performed an F-test across 5 groups for each end-point. The QA measurments on the individual micro-arrays all look good. One of the first things I do in my work-flow is take a look at the p-valued distribution. it is my understanding that, if the findings are due to chance alone, the p-value distribution should be uniform. In

Peter Dalgaard writes (in response to a question about 2-way ANOVA with imbalance): > ... There are various > boneheaded ways in which people try to use to assign some kind of > SumSq to main effects in the presence of interaction, and they are all > wrong - although maybe not very wrong if the unbalance is slight. People keep saying this

Hi everyone, I have a question about calculating r-squared in R. I have tried searching the archives and couldn't find what I was looking for - but apologies if there is somewhere I can find this... I carried out a droughting experiment to test plant competition under limited water. I had: - 7 different levels of watering treatment (1 -7 - from most watered to least watered/) - 15

Bonjour. Ne pouvant pas afficher toutes les pages web avec html window, est-ce qu''il est possible à partir d''un bouton de coller l''url choisi dans la barre de recherche d''IE7 par exemple et d''ouvrir IE7 automatiquement à la page recherchée? Est-ce que vous pouvez me donnez quelques explications si cela est possible? merci sebastien

Hi, A have small technical question about the calculation of R-squared using lm(). In a study case with experimental values, it seems more logical to force the regression line to pass through origin with lm(y ~ x +0). However, R-squared values are higher in this case than when I compute the linear regression with lm(y ~ x). It seems to be surprising to me: is this result normal ? Is there

Good morning. I wish to plot two data on the same axis. I tried plot(x,y, type = "l") for the first and tried to use lines or points(x,y, lty = 2, col = 4) to add or plot the second data on alongside the first. However, what I got was not encouraging. I have attached the two data and would be pleased if anybody could be of help. Thank you Best regards Ogbos -------------- next part

Hi, I was looking at the thread on the 'null' cipher and was wondering if you have done anything one it, if you need any help (or just encouragement ;)). I'm a radio experimenter (radio ham) and currently I'm using ssh illegally over the air, as I'm not legally allowed to encrypt sessions over the Radio Link. Being able to have strong authentication with no

http://www.thirdlane.com/screenshots.htm (Asterisk PBX Manager from Thirdlane) looks like a great program for "eye candy configuration" of Asterisk. However it costs lost of $, and I'm currently only an "experimenter" so to speak. Anyone advice of a decent alternative that is similar?? Currently, we only have VOIP connections, but will have a couple of Digium

http://www.thirdlane.com/screenshots.htm (Asterisk PBX Manager from Thirdlane) looks like a great program for "eye candy configuration" of Asterisk. However it costs lost of $, and I'm currently only an "experimenter" so to speak. Anyone advice of a decent alternative that is similar?? Currently, we only have VOIP connections, but will have a couple of Digium

[sorry if this arrives in duplo: it doesn't show up in the archives and it seems that the address I posted this from originally is no longer functional] Hello, I think I could do with some suggestions concerning the following problem. I have data from a set of experiments on motion sickness where for each subject, I have 1) personal data like age and gender 2) a subjective rating of