similar to: Simulation Questions

I'm new to R, but I'm trying to write a program for a dissertation that generates a dataset as follows... subject=1:1000 treat=rbinom(1*1000,1,.13) gender=rbinom(1*1000,1,.5) eth=runif(1*1000, min=1, max=4) cogat=rnorm(1*1000, 100, 16) map=rnorm(1*1000, 200, 9) simtest=data.frame (subject=subject, treat=treat, gender=gender, eth=round(eth,digits=0),

# R for Windows will not send your bug report automatically. # Please copy the bug report (after finishing it) to # your favorite email program and send it to # # r-bugs at r-project.org # ###################################################### This report is joint from Richard Heiberger <rmh at temple.edu> and Burt Holland <bholland at temple.edu>. Burt Holland is the coauthor

Dear R-users I analysed the same data set by two different ways; analysis of covariance by using lm and anova functions and multiple comparison by using simtest function in the multcomp library. The output from the analysis of covariance is; > y<-lm(D~Cond+Q1,data=x) > anova(y) Analysis of Variance Table Response: D Df Sum Sq Mean Sq F value Pr(>F) Cond 2

Hi, I've got a problem building a contrast matrix for the Dunnet contrast in package multcopm. The following works fine: > summary(simtest(adiff ~ trial)) Simultaneous tests: Dunnett contrasts Data: adiff by trial Contrast matrix: trial1 trial2 trial3 trial4 trial5 trial2-trial1 -1 1 0 0 0 trial3-trial1 -1 0 1 0 0

Hi, all: Two examples are shown below. I want to use the multiple comparison of Dunnett. It succeeded in upper case "example 1". However, the lower case "example 2" went wrong. In "example 2", the function pmvt return NaN, so I cannot show this simtest result. Is there any solution? (I changed the variable "maxpts" to a large number in front of the

Hi! I use simtest fonction of multcomp package to compile a Dunnett's test. I have 10 treatments and one control group, so i create a matrix with: m<-matrix(0,10,11) m[1,1]<--1 m[1,2]<-1 m[2,1]<--1 m[2,3]<-1 m[3,1]<--1 m[3,4]<-1 m[4,1]<--1 m[4,5]<-1 m[5,1]<--1 m[5,6]<-1 m[6,1]<--1 m[6,7]<-1 m[7,1]<--1 m[7,8]<-1 m[8,1]<--1 m[8,9]<-1

Hello, I have problems with installation of my second DomU. I cannot ping it from the outside world and I cannot access outside world from this domain. I see it from Dom0 and my first DomU, same as my second DomU can ping Dom0 and my first DomU. I use xen 3.0.3 on ubuntu 6 LTS. DomU is running on debian. I also network-bridge script in the xend config.

Dear WizaRds, I am sorry to bother you with a newbie question, but although I tried to solve my problem using the various .pdf files (Introduction, help pages etc.), I have come to a complete stop. Please be so kind as to guide me a little bit along my way of exploring multivariate analysis in R. I want to test wether the means-vector mu1 of X, consisting of the means per column of that matrix

Hello all I am working with the package multcomp but I have problems with the function simtest; the program say that can not find this function, nevertheless I doesn't have any problem with the function glht that it is in the same package. Someone knows what could be the problem? Thank you [[alternative HTML version deleted]]

Dear R-listers, I am trying to compute simultaneous confidence intervals with simint from the package multcomp. 230 measures (abundance) have been taken in 23 sites (factor) of a data.frame (donnees: a file can be sent on request, saved with save(donnees,file="donnees")). I would like to get all pairwise comparisons with : mc<- simint(ren~ID,type="Tukey",data=donnees) I

hello, after an anova I use pairwise.t.test(), it gives only p.value and I want the t.stat. I try to get these by computing the Welch approximation of the degree of freedom and using the qt(p.value,df) function but when I test this method with t.test results (the function gives p.value and t.test), I doesn't find the same t.stat. I also use the simtest(x~y,type="Tukey") function

Hi,. We have 2 values (first formant F1, second formant F2) for a given phoneme for six languages. We want to see whether the languages are significantly different one from another for this given phoneme. We have done a manova on our data and it works well, but we doesn't allow us to see which pair of languages are different. If we have only one formant for the phoneme, we would use

hi list, i have to ask you again, having tried and searched for several days... i want to do a TukeyHSD after an Anova, and want to get the adjusted p-values after the Tukey Correction. i found the p.adjust function, but it can only correct for "holm", "hochberg", bonferroni", but not "Tukey". Is it not possbile to get adjusted p-values after

Dear R-users, We expect to develop statistic procedures and environnement for the computational analysis of our experimental datas. To provide a proof of concept, we plan to implement a test for a given experiment. Its design split data into 10 groups (including a control one) with 2 mesures for each (ref at t0 and response at t1). We aim to compare each group response with control response

Hello R users, I am having difficulting getting multcomp to run. I have a dataframe attached with a numeric variable q12a and a numeric variable quota (which is really a classification variable). quota has 10 levels and unequal sample sizes. a12a has some missing data. I am interested in doing pairwise testing across the 10 quota groups on q12a. Using the ctest package the following code

I've never seen anything written about multiple comparisons, as in the multcomp package or with TukeyHSD, but using a glm. Do such procedures exist? Are they sensible? Are there any packages in R that implement such comparisons? Thank you. -- Ken Knoblauch Inserm U371 Cerveau et Vision 18 avenue du Doyen Lepine 69675 Bron cedex France Tel: +33 (0)4 72 91 34 77 Fax: +33 (0)4 72 91 34 61

Hi, I use glht to make multcomp, using Tukey, from a glm model. It is possible to get a more simplified output of result? Somethink like ordering by letters. Thanks Ronaldo -- Human kind cannot bear very much reality. -- T. S. Eliot, "Four Quartets: Burnt Norton" -- > Prof. Ronaldo Reis J?nior | .''`. UNIMONTES/Depto. Biologia Geral/Lab. de Ecologia | : :' :

Hi group, I'm trying to do a Tukey test to compare the means of a factor ("treatment") with three levels in an lme model that also contains the factors "site" and "time": model = response ~ treatment * (site + time) When I enter this model in csimtest, it takes all but the main factor "treatment" as covariables, not as factors (see below). Is it

Dear list, I have been studying R and I would like the aid of more experienced to solve the problems of the analysis below: r = gl(3, 8, label = c('r1', 'r2', 'r3')) e = rep(gl(2, 4, label = c('e1', 'e2')), 3) y = c(26.2, 26.0, 25.0, 25.4, 24.8, 24.6, 26.7, 25.2, 25.7, 26.3, 25.1, 26.4, 19.6, 21.1, 19.0, 18.6, 22.8, 19.4, 18.8, 19.2, 19.8, 21.4,

I would like estimate a number of contrasts from a one-way ANOVA model. I see that the lm command has a contrasts option, but I can't figure out how to use it! Any help that can be offered would be greatly apreciated. Here is my model statement: Model<-lm(log2PM~P+T+P*T) where P has 16 levels, T(treatment) has 12 levels and I am interested in looking at different treatment