Re: [R] Converting a categorical variable to multiple dichotemousvariables

From: Jeremy Hetzel <jthetzel_at_gmail.com>
Date: Tue, 12 Apr 2011 08:53:25 -0700 (PDT)


Shane,

Does this work?

# Your simulated data

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)
growth=0
simtest=data.frame (subject=subject, treat=treat, gender=gender, eth=round(eth,digits=0),
cogat=round(cogat,digits=0),map=round(map,digits=0),growth) simtest<-transform(simtest,
growth=rnorm(1000,m=ifelse(treat==0,0.1,0.5),s=0.03)) head(simtest)

# From Pablo above, but with unique() instead of levels(), since your eth is
numeric, not a factor
simtest.dichotomous <- outer(simtest$eth, unique(simtest$eth), "==") + 0

# Combine

simtest.combined <- cbind(simtest, simtest.dichotomous)

head(simtest.combined)

Jeremy

Jeremy Hetzel
Boston University



R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Tue 12 Apr 2011 - 16:02:41 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 12 Apr 2011 - 16:10:30 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive