similar to: Anova problem with order of terms in model

Good morning, suppose the following: m <- round(matrix(rnorm(16), ncol=4), 3) a <- rev(c(0.01, 0.025, 0.05, 0.1)) rownames(m) <- a colnames(m) <- c("0.25,0.75", "0.4,0.6", "0.1,0.9", "0.4,0.9") m 0.25,0.75 0.4,0.6 0.1,0.9 0.4,0.9 0.1 1.034 -0.119 -1.213 0.619 0.05 0.035 1.074 0.525 1.671 0.025 -1.687 0.960

Dear R users I have a question related to the interpretation of results based on GAMMs using Simon Woods package gamm4. I have repeated measurements (hours24) of subjects (vpnr) and one factor with three levels (pred). The outcome (dv) is binary. In the first model I'd like to test for differences among factor levels (main effects only): gamm.11<-gamm4(dv ~ pred +s(hours24), random = ~

Hi all and thanks for your time in advance, I can't figure out why summary.lmrob complains when lmrob is used on a zooreg object. If the zooreg object is converted to vector before calling lmrob, no problems appear. Let me clarify this with an example: >library(robustbase) >library(zoo) >dad<-c(801.4625,527.2062,545.2250,608.2313,633.8875,575.9500,797.0500,706.4188,

Hello, I would like to perform a multivariate regression analysis to model the relationship between m responses Y1, ... Ym and a single set of predictor variables X1, ..., Xr. Each response is assumed to follow its own regression model, and the error terms in each model can be correlated. Based on my readings of the R help archives and R documentation, the function lm() should be able to

Hi, I am using ccf but I could not figure out how to calculate the actual lag in number of periods from the returned results. The documentation for ccf says:"The lag is returned and plotted in units of time". What does "units of time" mean? For example: > x=ldeaths > x1=lag(ldeaths,1) > results=ccf(x,x1) > results Autocorrelations of series 'X', by lag

Hi, I have the following dataset and I would like to create a dotplot with confidence limits: CAT1 CAT2 MEAN Lower Upper 1 1 1 0.619 0.392 0.845 2 1 10 1.774 1.030 2.518 3 1

On Wed, Jun 18, 2008 at 6:32 AM, Joost Hietbrink <joost at joopp.com> wrote: > Hi Sam Roberts, > First of all. Thanks for the Vpim gem. We use it at www.yelloyello.com and > it works great. I''m glad to hear that. > We''ve only encountered the following error: > # NoMethodError (undefined method `to_str'' for []:Array): > # >

Dear gracious R community, I would like to produce charts of phytoplankton biomass changes through time. Each species has a line, and the biomass varies in mirror form along the line for each species along the X time axis. Here is an example of what I'd like to do: http://www.er.uqam.ca/nobel/r30240/Succession.jpg Any suggestions? Thanks David Bird UQAM, Montreal

Dear R community members, I think I am asking a very simple question, but I really looked up in the faqs and manuals and found nothing helpful. I am trying to plot a data frame with the following structure (this is just a small extract): glo conc odor line series X1 X2 X3 X4 X5 X6 X7 X8 X9 X10 X11 X12 X13 1 0 AIR LN1 UP -0.488

http://surf.chimie.uqam.ca/~kruus/vorbis/ now contain branch_beta3 and trunk subdirs that are updated daily (until branch_beta3 becomes historical next week, like the other snapshots at that URL). to web-browse the source code, go to http://surf.chimie.uqam.ca/~kruus/vorbis/lxr/http/source erik (back after audio-video sync miniproject. Yay!) --- >8 ---- List archives:

Hola, me estoy crispando con un tema que no es la primera vez que he resuelto pero (extraamente) ahora no puedo resolver. Estoy instalando un paquete que no est disponible en la versin 4 de R: packageurl <- "http://cas.uqam.ca/pub/R/bin/macosx/el-capitan/contrib/3.6/CASdatasets_1.0-10.tgz" install.packages(packageurl, repos = NULL, dependencies = TRUE, type="source")

I have been a member of the mailing list for a couple months now. I have not contributed anything because it has been difficult to gain a toehold from which I could really understand the concepts used in Ogg Vorbis. I would like to help out, but it is hard to get an idea of the overall architecture with docs like http://www.xiph.org/ogg/vorbis/doc/vorbis.html , in which every link is dead.

I am trying to use t.test on the following data: date type INTERVAL nCASES MTF SDF MTO SDO nFST MF nOBS MO MB BIASCV BIASEV ME MAE RMSE CRCF 2001-06-15 avn GE1.00 4385 0.246 0.300 1.502 0.556 1367 1.373 4385 1.502 1.471 0.285 0.164 -1.256 1.266 1.399 0.056 2001-06-15 avn

I've looked at how to speed up the transpose function in R (ie, t(X)). The existing code does the work with loops like the following: for (i = 0; i < len; i++) REAL(r)[i] = REAL(a)[(i / ncol) + (i % ncol) * nrow]; It seems a bit optimistic to expect a compiler to produce good code from this. I've re-written these loops as follows: for (i = 0, j = 0; i<len; i +=

Hi, I'm trying to plot bar graphs with stacked columns marked with askterisks for certain columns, which follow certain criteria. I've gotten the stacked bar plots with Data set A (please refer to code below). However, I haven't figured how to put asterisks on columns based on some extra information (such as "used/unused") as in Data set B i.e. asterisks at the

Hi I’m working in behavioural ecology I’m interesting in doing a multi-factorial MANOVA with a random effect as of independent variable. Do you know if it’s possible to do that and how I can do it? Thanks Julien Julien Martin Laboratoire de Denis Réale Chaire de Recherche du Canada en écologie comportementale Canadian Research Chair in behavioural ecology Département des sciences

#include <float.h> should define DBL_MIN and FLT_MIN, as well as DBL_MIN_EXP and FLT_MIN_EXP (such that 10eFLT_MIN_EXP is a valid float) This is a standard header file. But gcc now hides it away (mine was in /usr/lib/gcc-lib/i386-linux/<compiler-version>/include/ ) BTW - any corrections for my surf.chimie.uqam.ca/~kruus/vorbis/ stuff? Erik Kruus <kruus@on2.com> --- >8

Martin, You got some fine response already from other r-users. Indeed, centering explanatory variables is quite common and it really can have benefits in numerical accuracy. However one might also use centering for response variable and even standardization i.e. x - mean(x) / sd(x). This might increase the robustness of multivariate analysis, where variables have very different scales. For

Hi. I've put up some docs related to my first encounter with the vorbis source code at http://surf.chimie.uqam.ca/~kruus/vorbis/encoder/index.html It documents a newbie's first look and [serious] attempts to trace through examples/example_encoder.c . Don't laugh too hard! Even though the docs are quite un-doc-like, I've stopped working on them indefinitely, since Monty

Dear useRs: I use pROC package to compute the bootstrap C.I. of AUC. The command was as follows: roc1<-roc(all$D,all$pre,ci=TRUE,boot.n=200) However, the result was: Area under the curve: 0.5903 95% CI: 0.479-0.7016 (DeLong) Why the C.I. was computed by the Delong Method? Yao Zhu Department of Urology Fudan University Shanghai Cancer Center Shanghai, China [[alternative HTML version