similar to: aggregate by column names

Hi, I'm wondering if you can help me, this is a really simple query but I keep getting confused. I have run a GLM to see how boldness varies over time following a particular treatment. The results are as follows... Call: glm(formula = boldtwentyfour ~ treatment + boldcontrol) Deviance Residuals: Min 1Q Median 3Q Max -1.7577 -0.5469 0.0456 0.5515 1.5327

Hi, I have been trying unsuccessfully to plot data using different colors based on a variable within a subset of an imported file. The file I am reading is about 20000 lines long and has a column (in the example called FILE) that contains approximately 100 unique entries. I would like to plot a subset of the data from the file and key the color from the FILE column, This is what my file looks like

Hello, I'm running mixed models in GAMM4 with 2 (non-nested) random intercepts and I want to include a spline term for one of my exposure variables. However, when I include a spline term, I always get reported degrees of freedom of less than 1, even when I know that my spline is using more than 1 degree of freedom. For example, here is the code for my model: >

Dear R list, I've been attempting to interpret the results from a three-way ANOVA. I think I understand contrasts and the R defaults for these (treatment contrasts). My question is: what is the intercept in this test? As far as I can tell, its NOT the expected value of a point that belongs to the first level of all three explanatory factors (because there is only one point that satisfies these

Hi, i've been trying to apply the lme apprach to growth curves of children, but lme keeps running for ever and ever as soon as I use a reasonable basis. First Example: Data are 39 boys from the Berkeley growth study, each one measured 31 times at the ages of 1.00 1.25 1.50 1.75 2.00 3.00 4.00 5.00 6.00 7.00 8.00 8.50 9.00 9.50 10.00 10.50 11.00 11.50 12.00 12.50 13.00 13.50

Using caret on the Titanic data from Kaggle, I tried various models, including rfRules which produces a model, partly described as such: > caret.rfRules.cv$finalModel $model len freq err [1,] "2" "0.0368" "0" [2,] "2" "0.032" "0.05" [3,] "2"

I had found the author's (Wolfgang Viechtbauer) earlier meta-analytic code in R, MiMa, useful. so I have been exploring metafor using an example dataset from MiMa. metafor provides a lot more. However, MiMa provided parameter estimates, standard errors, z values, etc. for individual moderators in the meta-analysis, but I don't see how to obtain these from metafor. Have you any help

The test is manova. I tried to use manova() function, I used the code below:fit <- manova(Y ~ X)summary(fit, test="Wilks")but I get p values for intercept and regression coefficient as in anova() function, not for the hull model. Date: Mon, 7 Feb 2011 00:57:43 -0800 Subject: Re: [R] FW: multivariate regression From: djmuser@gmail.com To: denizsigirli@hotmail.com CC:

How can I run multivariate linear regression in R (I have got 3 dependent variables and only 1 independent variable)? I tried lm function, but it gave different R2 and p values for every dependent variable. I need one R2 and p value for the model. [[alternative HTML version deleted]]

Hello, I?m using zero-inflated Poisson regression via the zeroinfl() function to analyze data on seed-set of plants, but for some reason, I don?t seem to be getting the output for all three levels of my two categorical predictors. More about my data and model: My response variable is the number of viable seeds (AVInt), and my two categorical predictors are elevation (Elev) and Treatment

Dear Sir or Madam? I am a doctor of urology,and I am engaged in developing a nomogram of bladder cancer. May I ask for your help on below issue? I set up a dataset which include 317 cases. I got the Binary Logistic Regression model by SPSS.And then I try to reconstruct the model ?lrm(RECU~Complication+T.Num+T.Grade+Year+TS)? by R-Project,and try to internal validate the model through

Hi all, I computed a two way anova with two within-subjects fixed effects: diffi (2 levels) and cond (4 levels) and one random effect subject, and one dependent variable accuracy: summary(aov(accuracy ~ diffi*cond + Error(sub/(diffi+cond)), data=phd.df)) Error: sub Df Sum Sq Mean Sq F value Pr(>F) Residuals 4 136.549 34.137 Error: sub:diffi Df Sum Sq Mean Sq F value

All, I'm attempting to run this: table.CAPM(series[,"Strat.Return",drop=FALSE],series[,"spy.Return",drop=FALSE]) and getting this error Error in as.vector(data[, i]) : subscript out of bounds I've searched around and cannot find a solution to the problem.  I've used this in the past without problem and I'm not sure what I did different this time. Other

Hello, I would like to perform a Partial least square discriminate analysis (PLSDA) in R. To do this I use the package mixOmics. I could perform the PLSDA in R. however I would also like to perform a leave-one-out cross validation in order to assess the performance of my model. My supervisor told me that I should focus on the R2/Q2 ratios. However when I read the instruction for running the

Dear ALL, I'm trying to figure out what the percentage effects are in a logistic regression. To be more clear, I'm not interested in the effect on y of a 1-unit increase in x, but on the percentage effect on y of a 1% increase in x (in economics this is also often called an "elasticity"). For example, if my independent variables are in logs, the betas can be directly

I rebuild mesa with debug symbols, and now top functions using CPU looks like this: CPU: AMD64 family15h, speed 3800 MHz (estimated) Counted CPU_CLK_UNHALTED events (CPU Clocks not Halted) with a unit mask of 0x00 (No unit mask) count 100000 samples % image name symbol name ------------------------------------------------------------------------------- 222978 45.1489

G'Day Tanya, Is it too late to bring in the following patches to fix some major brokenness in the AuroraUX tool chain for 2.6? http://llvm.org/viewvc/llvm-project/cfe/trunk/lib/Driver/Tools.cpp?r1=84468&r2=84469&view=diff&pathrev=84469 http://llvm.org/viewvc/llvm-project/cfe/trunk/lib/Driver/Tools.cpp?r1=84265&r2=84266&view=diff&pathrev=84266

Hi Tanya, > 1) Compile llvm from source and untar the llvm-test in the projects > directory (name it llvm-test or test-suite). Choose to use a > pre-compiled llvm-gcc or re-compile it yourself. I compiled llvm and llvm-gcc with separate objects directories. Platform is x86_64-linux-gnu. > 2) Run make check, report any failures (FAIL or unexpected pass). Note > that you need to

On Oct 20, 2009, at 6:02 AM, Duncan Sands wrote: > Hi Tanya, > >> 1) Compile llvm from source and untar the llvm-test in the projects >> directory (name it llvm-test or test-suite). Choose to use a pre- >> compiled llvm-gcc or re-compile it yourself. > > I compiled llvm and llvm-gcc with separate objects directories. > Platform is x86_64-linux-gnu. > Ok.

LLVMers, 2.6 pre-release2 is ready to be tested by the community. http://llvm.org/prereleases/2.6/ If you have time, I'd appreciate anyone who can help test the release. To test llvm-gcc: 1) Compile llvm from source and untar the llvm-test in the projects directory (name it llvm-test or test-suite). Choose to use a pre- compiled llvm-gcc or re-compile it yourself. 2) Run make check,