similar to: partial SS for anova

Hello all, I have data: Genes time rep vart dye y trt 130911 sa1-d07 030min 1 col g 9.636244 o 145771 sa1-d07 030min 1 col r 8.107577 c 93335 sa1-d07 030min 1 ler g 7.409566 o 94821 sa1-d07 030min 1 ler r 5.107160 c 10119101 sa1-d07 030min 2 col g 8.336862 o 11605101 sa1-d07 030min 2 col r 7.824530 c 725313 sa1-d07

Hello all, I have problem with setting up random effects. I have a model: y=x1+x2+x1*x2+z1+z1*x2 where x1, x2, x1*x2 are fixed effects and z1, z1*x2 are random effects (crossed effects) I use library(nlme) 'lme' function. My question is: how I should set up random effects? I did lme(y~x1+x2+x1:x2, data=DATA, random=~z1+z1:x2, na.action='na.omit') but it did not work.

I'm trying to configure Squid ntlm autentication on Samba4 DC. I followed Squid and Samba's documentation and i got success when I login with user natalia.silva, but if I log with natalia.vaz i get the error -- Nat?lia Vaz Silva Administradora de redes

Hello, I was using version 2.13.2 and I have just downloaded the latest version 2.14.1. However, I''m trying to install the packages I was using and when I look for them in the packages list, I can´t find many in the CRAN binaries (e.g. "vegan"). I do find them in the CRAN sources but the installation fails. I tried downloading the version 2.14.0 and I had the same problem. I

I have two questions I would appreciate assistance with: (1) I believe "by" is the command used to split a file. In the following example, "mydat" is the dataframe , "group" the variable I want to split the analysis by & "WK1FREQ,WK2FREQ" the variables attach(mydat) by (GROUP) cor.test (WK1FREQ,WK2FREQ) I have also tried: (by (GROUP) cor.test

Dear all, I'm using R version 3.2.0 in Ubuntu 12.04 LTS. While trying to install any package I get the error: "/usr/bin/ld: cannot find -lgfortran" If I open a Terminal and run "ldconfig -p | grep libgfortran"?, I get:? "libgfortran.so.3 (libc6,x86-64) => /usr/lib/x86_64-linux-gnu/libgfortran.so.3" I've been an R user for two years and this is the

I'm an SPSS user trying to make the transition to R. Can someone help me translate the following SPSS code into R?: GLM Total_tp1 Total_tp2 WITH Age Sex   /WSFACTOR=Time 2 Repeated   /METHOD=SSTYPE(3)   /CRITERIA=ALPHA(.05)   /WSDESIGN= Time   /DESIGN= Age Sex Age*Sex. Also. can anyone recommend any resources to help SPSS users learn to things in R? Thanks, -kristi [[alternative HTML

Hello, R users, I have two factors (treat, section) anova design experiment where there are 3 replicates. The objective of the experiment is to test if there is significant difference of yield between top (section 9 to 11) and bottom (section 9 to 11) of the fruit tree under treatment. I found that there are interaction between two factors. I wonder if I can contrast means from levels of

Dear Thomas Lumley, and R-help list members, I have read your article "Network meta-analysis for indirect treatment comparisons" (Statist Med, 2002) with great interest. I found it very helpful that you included the R code to replicate your analysis; however, I have had a problem replicating your example and wondered if you are able to give me a hint. When I use the code from the

Hello, R experts, I have a list called dataHP which has 30 elements (m1, m2, ..., m30). Each element is a 7x6 matrix holding yield data from two factors experimental design, with treatment in column, position in row. For instance, the element 20 is: dataHP[[20]] col1 col2 col3 trt1 trt2 trt3 [1,] 22.0 20.3 29.7 63.3 78.5 76.4 [2,]

Dear R people: I wish to switch from SPSS to R, but there is one particular type of ANOVA design that cannot be done in R. Or more likely, it can be done, but it is nowhere documented. The problem is typical for psychologists: You have a repeated measures design with different groups of subjects. Now, this can be done with the aov command, but the number of subjects in both groups must be

I am confused whether Student's sleep data "show the effect of two soporific drugs" or Control against Treatment (one drug). The reason is the next: > require(stats) > data(sleep) > attach(sleep) > extra[group==1] numeric(0) > group [1] Ctl Ctl Ctl Ctl Ctl Ctl Ctl Ctl Ctl Ctl Trt Trt Trt Trt Trt Trt Trt Trt Trt [20] Trt Levels: Ctl Trt > sleep$group [1] 1 1 1 1 1

I am trying to define groupings from levels of factor variables and this the warning message that R give "& not meaningful for factors". The nature of my task is this. I have a variable stage which has the levels (1B, 2A, 2B) - these are the AJCC TNM stages of cancer, and another variable diameter with factor levels ("=< 4", "4 - 6.5, > 6.5; limit values are

Dear, I have a dataset with 4 subjects (see ID in example), and 4 treatment (see TRT in example) which are tested on 2 locations and in 3 blocs. By using Lattice dotplot, I made a graph that shows the raw data per location and per bloc. In that graph, I would like to have a reference line per bloc that refers to the first treatment (T1). However, I can not find how to do that. I can make

Hello, I probably have a syntax error in trying to generate an expected survival curve from a weighted cox model, but I can't see it. I used the help sample code to generate a weighted model, with the addition of a "weights=albumin" argument (I only chose albumin because it had no missing values, not because of any real relevance). Below are my code with the resulting error

Hello, It seems that, unlike Splus, the anova.lm function of R does not take the argument ssType=3 into account. How can one get an anova table with the adjusted sum of squares in R? Thanks in advance, Patrick Giraudoux University of Franche-Comté Department of Environmental Biology EA3184 af. INRA F-25030 Besançon Cedex tel.: +33 381 665 745 fax.: +33 381 665 797 http://lbe.univ-fcomte.fr

Hello, I make a samba server + ldap. I create some users and get login in mydomain with success. After login I can't write in paste user.But is denied. #some users created brasil:x:1012:513:Brasil Telhados:/dev/null:/bin/false atleticomg:x:1013:513:Clube Atletico MG:/home/atleticomg:/bin/bash #groups crated informatica:*:1001:atleticomg,natalia,faria

Hello, I'm having trouble in using "assign(paste ..." command . I could create several dataframes following trinomial distribution using it but it could not be used to check their row means of the created dataframe. For example, the following works: probTrt=matrix(0,4,3); probTrt; #malf, death, normal probTrt[1,]=c(0.064,0.119,0.817);#for Trt 1 probTrt[2,]=c(0.053,0.125,0.823);#for

You need to give the model formula that gave your output. There are two sources of variation (at least), within and between locations; though it looks as though your analysis may have tried to account for this (but if so, the terms are not laid out in a way that makes for ready interpretation. The design is such (two locations) that you do not have much of a check that effects are consistent over

Dear R users, I am trying to use xyplot to draw group mean and CI. The following is the sample code. But I want: 1. Use different colors and symbols to draw individual points, CI and the lines connect group means from different time points; 2. Add jitters to x axis to allow CIs not be overlapped Could anyone modify the attached code to achieve this? Thanks library(lattice)