similar to: Nested anovas in R not doing what they ought to...

Hi, there I am not a member on your mailing list, so could you please cc: a response to me? I am currently attempting to perform an ANOVA with both nested and normal factors. My problem is that R is treating my nested factors the exact same way as it would interaction terms. My output for the nested model is as follows: > nested <- anova(lm(ltotinv ~ habitat + fish/lake +

Hi there Please excuse this elementary question, but I have been fumbling with this for hours and can't seem to get it right. I have a nested anova, with random factor "lakefac" nested within factor "fishfac" (fixed), with an additional fixed factor "Habfac". If I consider everything as fixed effects, it's addmittedly not the correct model, but I can at

Dear all R-users, I wanted to fit a Garch(1,1) model to a dataset by: >garch1 = garch(na.omit(dat)) But I got a warning message while executing, which is: >Warning message: >NaNs produced in: sqrt(pred$e) The garch parameters that I got are: > garch1 Call: garch(x = na.omit(dat)) Coefficient(s): a0 a1 b1 1.212e-04 1.001e+00 1.111e-14 Can any one

Dear R users: When I create a table without Major Column headings, my *regular* column headings appear correct in the typeset latex file. The major row heading and row groups are as they should. w <- latex(mytab,title="",file="tab/my.tex",ctable=TRUE,caption="Descriptive statistics by

Greetings. I'm trying to recreate in R some regression models I've done in SAS, but I'm not getting the same results. My advisor suspects this may be due to differences in precision between R and SAS. Does anyone know where I can find specifications for R's type double? (It doesn't seem to be in the R Language Definition.) Thanks in advance for any help anyone can

Dear all R user, Please forgive me if my problem is too simple. Actually my problem is basically Statistical rather directly R related. Suppose I have return series ret with mean zero. And I want to fit a Garch(1,1) on this. my is r[t] = h[i]*z[t] h[t] = w + alpha*r[t-1]^2 + beta*h[t-1] I want to estimate the three parameters here; the R syntax is as follows: #

I am wondering how to customize a pairwise comparisons plot of a factorial ANOVA, without doing a lot of manual manipulation of a TukeyHSD object. The customizations I'd like are: 1. The aov used log-transformed response data, but I'd like to plot the intervals on their original, untransformed scales 2. Plot all the main and interaction effects together, rather than in a separate

I'm having an issue with occasionally slow-running calls to normalizePath. If the path is a non-existent UNC path, then normalizePath sometimes takes 6 or 7 seconds to run, rather than its usual few microseconds. My big problem is that I can't reliably reproduce this across machines. The example below generates one or two slow runs out of 10000 on my Windows machine. I haven't been

Hi all, I've started looking at the GlobalMerge pass, enabled by default on ARM and AArch64. I think we should reconsider that, at least for AArch64. As is, the pass just merges all globals together, in groups of 4KB (AArch64, 128B on ARM). At the time it was enabled, the general thinking was "it's almost free, it doesn't affect performance much, we might as well use it".

Hi, I have a data frame with 6563 observations. I can run a regression with loess using four explanatory variables. If I add a fifth, R crashes. There are no missings in the data, and if I run a regression with any four of the five explanatory variables, it works. Its only when I go from four to five that it crashes. This leads me to believe that it is not an obvious problem with the data,