similar to: Anova interaction not tested

Dear Colleagues, Apologies for a long email to ask what I feel may be a very simple question; I figure it's better to overspecify my situation. I was asked a question, recently, by a colleague in my department about pre-aggregating variables, i.e., computing the mean of defined subsets of a data frame. Naturally, I thought of the 'by' and 'tapply' functions, as

Hey Folks: New to R and learning as I go, but really, I know just enough to get myself into trouble. I've waded through everything up till now, and don't see anything in the search that is directly helpful for the POS that I've created. The GOAL: All I want in the world is a program that performs 1-way ANOVA's on every column in a data set (taking the first column as the

Hi all, I appreciate the help this list has given me before. I have a question which has been perplexing me. I have been working on doing a Bayesian calculating inserting studies sequentially after using a non-informative prior to get a meta-analysis type result. I created a function using three iterations of this, my code is below. I insert prior mean and precision (I add precision manually

Dear list, I would like to export a matrix to a TXT-File by using write.csv (not necessarily). Is there a way to add a header (with additional informations concerning the project) spanning multiple lines to this file before the actual data are listed up? Should look like this: date: filename: number of permutations: ------------ data (as a matrix) Any suggestions? Thnx in advance.

Hello. I'm currently stuck with the same "what does interaction.depth really mean" stuff. Did you find out what the right answer is? Best regards, Boris Yangel. [[alternative HTML version deleted]]

Dear All, I am working with a generalized linear mixed-effects model with poisson error using the lme4 package in R. I created a model with the lmer function including some main effects, three two way interactions and two random effects. The model works well, but I have troubles when removing on of the two-way interactions. The Warning message: "In mer_finalize(ans) : false convergence

Suppose I have a small dataframe > aa Target Eaten ID 50 TPP 0 1 51 TPP 1 2 52 TPP 3 3 53 TPP 1 4 54 TPP 2 5 50.1 GPA 9 1 51.1 GPA 11 2 52.1 GPA 8 3 53.1 GPA 8 4 54.1 GPA 10 5 And I want to reshape it into ID TPP GPA 1 1 0 9 2 2 1 11 3 3 3 8 4 4 1 8 5 5 2 10 I realise that

When I try to install Rmpi, it has a dependency rsprng which, as its description says: Provides interface to SPRNG 2.0 APIs I installed it on Mepis by installing the appropriate debs with minimal drama. In the rpm world of Fedora, no such luxury. Installing rsprng fails because: checking sprng.h usability... no checking sprng.h presence... no checking for sprng.h... no Cannot find sprng 2.0

I got Rmpi to compile with little difficulty, but had a tricky time setting the LD_LIBRARY_PATH to use the OpenMPI libs. I now get a different error when I try to load Rmpi > require(Rmpi) Loading required package: Rmpi librdmacm: couldn't read ABI version. librdmacm: assuming: 4 libibverbs: Fatal: couldn't read uverbs ABI version. CMA: unable to open

Dear list members, i am currently want to install Rpy2 in a linux box which has R 2.4.0 installed RPy requries R 2.7.0 or above but i have no root previlleges so my question is how to install R 2.7.0 on my own directory? and replace the system installed R 2.4.0 when i input R command from the bash Shell? thanks for any reply. Best Regards! -- Daofeng Li,PhD Candidate China Agricultural

Or, to avoid accusing you of lying. what you think is "vanilla" probably isn't. What exactly did you do? On Unix-likes, I would do something like this echo 'options(repos=list(CRAN="cran.r-project.org"));install.packages("Rcpp")' | R --vanilla (or maybe https://cloud.r-project.org is better...) -pd > On 2 Feb 2018, at 08:15 , Jeff Newmiller

When i tried to install the hunspell package, I got this error message: Error: package ?Rcpp? 0.12.3 was found, but >= 0.12.12 is required by ?hunspell? So I set about installing a new version of Rcpp but I get this message: Error in unloadNamespace(pkg_name) : namespace ?Rcpp? is imported by ?dplyr? so cannot be unloaded How does one get around that? I tried installing Rcpp in a

Dear all,   I did ANOVA for many response variables (Var1, Var2, ....Var75000), and i got the result of p-value like below. Now, I want to select those predictors, which have pvalue less than or equal to 0.05 for each response variable. For example, X1, X2, X3, X4, X5 and X6 in case of Var1, and similarly, X1, X2.......X5 in case of Var2, only X1 in case of Var3 and none of the predictors in case

Hi all, I try to make a split-plot with poisson errors using glmmPQL, but I have some doubts about the model simplification. Look my system: Block = 3 blocks Xvar1 = 2 levels Xvar2 = 13 levels Yvar = Count data Response I need know about the behaviour of Var1, Var2 and interaction Var1:Var2. Look the levels: > levels(Xvar1) [1] "A" "B" > levels(Xvar2) [1]

Dear R-users, I’d like to make a nested anova on my data, and since I’m discovering both R and statistic, I’d like to be sure that I’m not doing something stupid. Here is my data: I’ve measured some variable responses (Y, for example leaf size) for different plants grown on three differing conditions (A, B and C). These plants come from two different species (sp1 and sp2) and within each species,

On Thu, 18-May-2017 at 05:46PM +0200, Martin Maechler wrote: |> ..... |> |> Being pretty "stretched" time wise currently, I'm happy for |> timezone-portable propositions to change the test. Meantime, anyone who lives where DST happpens in December who wants to get through the remaining tests can avoid this one by changing the line > stopifnot(length(fd) == 10,

On Wed, 17-May-2017 at 01:21PM +0200, Peter Dalgaard wrote: |> |> Anyways, you might want to |> |> a) move the discussion to R-devel |> b) include your platform (hardware, OS) and time zone info System: Host: MTA-V1-427894 Kernel: 3.19.0-32-generic x86_64 (64 bit gcc: 4.8.2) Desktop: KDE Plasma 4.14.2 (Qt 4.8.6) Distro: Linux Mint 17.3 Rosa Machine: System:

The 3.2.4 release had two annoyances which we would rather not have in an "ultra-stable" release, designed to hang around for the duration of the 3.3 series. One was a relatively minor Makefile issue affecting system using R's bundled lzma library. The other, rather more serious, affected printing and formatting of POSIXlt objects, which would unpredictably get the Daylight Savings

Hi, I am trying to reconcile anova table in R (summary(lm)) with individual t.test. datafilename="http://personality-project.org/R/datasets/R.appendix1.data" data.ex1=read.table(datafilename,header=T) #read the data into a table summary(lm(Alertness~Dosage,data=data.ex1)) gives: Call: lm(formula = Alertness ~ Dosage, data = data.ex1) Residuals:    Min     1Q Median     3Q    Max

Hi, I'm a beginner in both statistics and the usage of R. I have been stuck for days with a problem I hope somebody is willing to help me solving. I'm trying to create some kind of graphical plot etc. of my results of the following glm of interactions. I17<-glm(formula = response ~ var1+ var2 + var3 + var4 + var1:var3 + var1:var4, data = data) Is there any possibility to solve