similar to: lm on group

Dear Dr Skaug and R users, just discovered glmm.admb in R, and it seems a very useful tool. However, I ran into a problem when I compare two models: m1<-glmm.admb(survival~light*species*damage, random=~1, group="table", data=bm, family="binomial", link="logit") m1.1<-glmm.admb(survival~(light+species+damage)^2, random=~1, group="table", data=bm,

Hi, Using package gclus in R, I have created some graphs that show the trends within subgroups of data and correlations among 9 variables (v1-v9). Being interested for more details on these data I have produced also the var-covar matrices. Question: From a pair of two subsets of data (with 9 variables each, I have two var-covar matrices for each subgroup, that differ for a treatment on one

#subject: type III sum of squares - anova() Anova() AnovaM() #R-version: 2.12.2 #Hello everyone, #I am currently evaluating experimental data of a two factor experiment. to illustrate de my problem I will use following #dummy dataset: Factor "T1" has 3 levels ("A","B","C") and factor "T2" has 2 levels "E" and "F". The

Hi all, I have a very easy questions (I hope). I had measure a property of plants, growing in three different substrates (A, B and C). The rest of the conditions remained constant. There was very high variation on the results. I want to do address, whether there is any difference in the response (my measurement) from substrate to substrate?

Dear all, I often import data sets from other programs that contain user defined missing values. As an example consider a variable MXPLOCO with a user defined missing value of as -9. > summary(MXPLOCO) Min. 1st Qu. Median Mean 3rd Qu. Max. NA's -9.000 4.750 5.000 5.349 6.000 8.500 2.000 The following command successfully converts the value -9 to NA. However,

Hi, I'm wondering how I can make my patched up mplayer work with other MPlayer frontends, such as mplayer-plugin, GMPlayer, SMPlayer and KMPlayer. When I play a video using mplayer normally with: mplayer video.ext, It detects the compiz video output and works fine, however when it is called from other players (AFAIK, it just creates a window then asks mplayer to display in it) with mplayer

Dear colleagues, How can I do to "relevel" the intercept? I need that the treatment "Db" be the intercept, and have p-values for the comparisons with the others treatments. I used the function "relevel" but it did not work out to have what I want. Thanks for your help, Xavier T1 <- read.table(file.choose(), h=T) > head(T1) treatment replicate Time

Hi, (R version 2.15.0) I am running a pgm with 1 response (earlier standardized Y) and 44 independent vars (Xi) from the same data =a2: When I run the 'lm' function on single Xi at a time, the beta coefficient for let's say X1 is = -0.08 (se=0.03256) But when I run the same Y with 44 Xi-s with the 'step' function (because I left direction parameter empty, I assume a backward

Hi I am working on corpora of automatically recognized utterances, looking for features that predict error in the hypothesis the recognizer is proposing. I am using the glm functions to do logistic regression. I do this type of thing: * logistic.model = glm(formula = similarity ~., family = binomial, data = data) and end up with a model: > summary(logistic.model) Call:

Dear R developers, I was looking for a function analogous to base::diag() for getting and setting bands of a matrix. The closest I could find was Matrix::band(), but this was not exactly what I wanted for two reasons. Firstly, Matrix::band() returns a matrix rather than just the specified band. Secondly, Matrix::band() cannot be used for setting the values for a matrix band. Setting or

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

Hello, I need some help on a regression I am running please. I am running a multiple regression in R and I am getting weird outputs and would like your help in resolving it. This is the code I run: Reg = lm (Final$Y.t.-Final$Y.t.1. ~ Final$ Y.t.1. + Final$Cor + Final$Gov+ Final$Inv+ Final$TrOp + Final$Pop+ Final$Sch , data = Final) summary(Reg) And the output I get is pasted below. As

I have 11 possible predictor variables and use them to model quite a few target variables. In search for a consistent manner and possibly non-manual manner to identify the significant predictor vars out of the eleven I thought the option "select=T" might do. Example: (here only 4 pedictors) first is vanilla with "select=F" >

Shalini Raghavan 3M Pharmaceuticals Research Building 270-03-A-10, 3M Center St. Paul, MN 55144 E-mail: sraghavan at mmm.com Tel: 651-736-2575 Fax: 651-733-5096 ----- Forwarded by Shalini Raghavan/US-Corporate/3M/US on 08/16/2004 11:25 AM ----- Shalini

Let's say i have a = c(1, 4, 5) b = c(2, 6, 7) and i have matrix m, what's an efficient way of access m[1, 2], m[4, 6], m[5, 7] like of course m[a, b] = is not going to do, but what's an expression that will allow me to have that list? Thanks! -- View this message in context: http://www.nabble.com/Simplest-question-ever...-tf3329894.html#a9258932 Sent from the R help mailing list

Dear all, I am unable to understand how Hline() works in tabular(). I've read the vignette and the help page, and here this example compiles perfectly fine: latex( tabular( Species + Hline() + 1 ~ Heading()*mean*All(iris), data=iris) ) However, if I try it on my own data it fails. Consider this: set.seed(1) Xa <- data.frame(p=rep(c("First group","Second

Hello together, this is ma first post, so please aplogize me if post this in the wrong section. I have problem concerning ma two matrizes. After a regressione and so on, I got two matrizes Matrixres contains the results of ma calculation. Matrixr contains my detiene, which where Aldo used for the regression. Please ser the following code: #Datei einlesen residual =

There is a problem when importing an spss-file containing explicitly declared missing values in R using the read.spss function from the foreign package. I'm not sure these problems are the same in every version of spss, I am using the latest version 16.0.2. I included http://www.nabble.com/file/p18776776/missingdata.sav missingdata.sav and

I have a data set for different time intervals. The data has three comment lines before data for each time interval. For each time interval there are 500 data points. I want to change the dataset such that I have the following format: t1 t2 t3 ................ 0.00208 0.00417 0.00625 ................. a1 a2 a3 ...................

Hi all, I have this list of strings [1] "1 ,2 ,3" "4 ,5 ,6" Is there an efficient way to convert it to data.frame: V1 V2 V3 1 1 2 3 2 4 5 6 Like I can use strsplit to get to a list of split strings.. and then use say a = strsplit(mylist, ",") data.frame(V1 = lapply(a, function(x){x[1]}), V2 = lapply(a, function(x){x[2]}),.....) but i'm