similar to: Suggestion on Model Def - reg in nlme

Hi R Forum I am a newbie to R and I have been amazed by what I can get my team to accomplish just by implementing Scripting routines of R in all my team's areas of interest.. Recently i have been trying to adopt R scripting routine for some analysis with longitudanal data.. I am presenting my R script below that I have tried to make to automate data analysis for longitudanal data by employing

Dear R Group I have a character vector from which I want to select a few elements and create a new character vector. I need a GUI to do this in R Script. Can someone help? a<-c("A","B","C","D","E") ## I want to have a GUI in R that will display elements in "object a" as a drop down list.. from there I want to be able to select a

Dear R Group I had an observation that in some cases, when I use the randomForest model to create partialPlot in R using the package "randomForest" the y-axis displays values that are more than -1! It is a classification problem that i was trying to address. Any insights as to how the y axis can display value more than -1 for some variables? Am i missing something! Thanks Regards

Dear RGroup I have a requirement for which I am seeking help. I am looking at automating the last column calculation through R when having the data of the other columns as a dataframe, In excel I can do using the formula function as given below, however, hereagain for the number of observations that come under control, I need to write the formula seperately for each cell and then I can block

Dear R Group Can any body point me to a link from where I can get zip archive of the package HH for R 2.12.0? Thanks & Regards Vijayan Padmanabhan "What is expressed without proof can be denied without proof" - Euclide. Please visit us at www.itcportal.com ****************************************************************************** This Communication is for the exclusive use

Dear Group I have to populate a column in a R dataframe containing say 99 rows, where for every 3 row one row (randomly determined) of the column should take the value 1 and the other 2 row of that column should take the value 0. I used the command dataframe$new<-rep(sample(c(0,1,0),3,replace=FALSE),99) This command repeats the same patter of response across the entire column.. I would

Dear Group I was referring to a conjoint analysis scenario using R from the paper referred below: Agricultural Information Research 17(2),2008,86-94 available online at www.jstage.jst.go.jp/ This paper describes the data modelling of a conjoint study design based on conditional logit procedure. I understand that Heirarchical Bayes is asymptotically equivalent to Conditionallogit. However it

Dear R Users I am looking at carrying out Multivariate Supervised Discretization on a set of continuous attributes to determine cut off values on each attribute which would allow best explanation of a binary outcome variable. I understand that there is an algorithm for doing this (The HyperCluster Finder Algorithm from the paper -Muhlenbach F., & Rakotomalala R. (2002). Multivariate

Dear R Group I have a matrix object in R, where i want to retain only those columns which have each row from a different group.. see the example below.. team1<-c("a1","a2","b1","b2","b3","c1","c2") teams<-combn(team1, 3) I want a function which will delete all columns in the teams matrix that have more than 1 row

Dear R Group in package Deducer, how should i use the add.cross.strata.test() function to pass the chi.squared.test function as htest argument? ##########Example########### library(Deducer) dat<-data.frame(a=rnorm(100)>.5,b=rnorm(100)>0,c=rnorm(100)>(-.5)) tables<-contingency.tables( row.vars=a, col.vars=b, stratum.var=c,data=dat)

Hello everybody, problem solved, there was a typo. I wrote Type instead of Material Best ----Messaggio originale---- Da: angelo.arcadi@virgilio.it Data: 22-giu-2012 11.05 A: <r-help@r-project.org> Ogg: Error with glht function: Error in mcp2matrix(model, linfct = linfct) : Variable(s) 'Type' have been specified in 'linfct' but cannot be found in 'model'!

Dear R members, I tried to fit an lme model and to use the glht function of multcomp. However, the glht function gives me some errors when using weights=varPower(). The glht error makes sense as glht needs factor levels and the model works fine without weights=. Does anyone know a solution so I do not have to change the lme model? Thanks Sibylle --> works fine

Hi, I am experimenting with using glht() from multcomp package together with coxph(), and glad to find that glht() can work on coph object, for example: > (fit<-coxph(Surv(stop, status>0)~treatment,bladder1)) coxph(formula = Surv(stop, status > 0) ~ treatment, data = bladder1) coef exp(coef) se(coef) z p treatmentpyridoxine -0.063 0.939 0.161

I've run an ancova, edadysexo is a factor with 3 levels,and log(lcc) is the covariate (continous variable) I get this results > ancova<-aov(log(peso)~edadysexo*log(lcc)) > summary(ancova) Df Sum Sq Mean Sq F value Pr(>F) edadysexo 2 31.859 15.9294 803.9843 <2e-16 *** log(lcc) 1 11.389 11.3887 574.8081 <2e-16 ***

Dear R-user I used lme to fit a linear mixed model inlcuding weights=varPower(). Additionally I wanted to use glht to calculate Tukey-Kramer multiple comparision. error: > glht(modelF, linfct=mcp(Species="Tukey")) Error in glht.matrix(model = list(modelStruct = list(reStruct = list(SubPlot = -0.305856275920955, : ?ncol(linfct)? is not equal to ?length(coef(model))? >

Hi, I have a fairly simple repeated measures-type data set I've been attempting to analyze using the lme function in the nlme package. Repeated searches here and other places lead me to believe I have specified my model correctly. However, I am having trouble with post-hoc tests. From what I gather, other people are successfully using the glht function from the multcomp package to

Hi, i was working with this model > mq<-glm(rojos~edadysexo*zona*estacion,quasipoisson) and i get this minimal adequate model > anova(mq5,test="F") Df Deviance Resid. Df Resid. Dev F Pr(>F) NULL 518 64799 edadysexo 2 1556.5 516 63243 8.9434 0.0001524 *** zona 4

Dear list members, I get the following error when using the glht function to perform a post hoc analysis for an ANOVA with repeated measures: require(nlme) lme_H2H_musicians = lme(H2H ~ Emotion*Material, data=musicians, random = ~1|Subject) require(multcomp) summary(glht(lme_H2H_musicians, linfct=mcp(Type = "Tukey")), test = adjusted(type = "bonferroni")) Error in

Dear all, I was trying to understand how "multcomp" package works by running the examples given in the documentation. However I still don't understand when it comes to multiple comparison set by user (please refer to "Ksub" in the code). Therefore I run 2 other cases along with the original example (case 1), with the expectation I'll get the point from the output. The

I am running a linear model with achiev as the outcome and major as my iv (5 levels). The lm statement runs fine, but for the glht command I get the following error. I noted that someone else asked the same question a while back but received no reply. I am hoping someone might know what is happening. anovaf2<-lm(achiev ~ major, data=data_mcp) > pairwise<- glht(anovaf2,linfct =