Dear List, I have a table i have read into R: Name Yes/No John 0 Frank 1 Ann 0 James 1 Alex 1 etc - 800 different times. What i want to do is shuffle yes/no and randomly re-assign them to the name. I have used sample() and permute(), however there is no way to do this 1000 times. Furthermore, i want to copy the data into a excel spreadsheet in the same order as the data was input so i can build up a distribution of the statistic for each name. When i use shuffle the date gets returned like this - [1] 1 0 0 1 0 1 0 0 0 1 1 1 1 1 1 1 1 1 0 1 0 0 0 0 0 0 1 0 0 0 1 0 1 [34] 0 1 0 0 0 0 0 0 0 1 0 1 1 0 1 0 0 0 0 0 0 1 1 0 0 1 0 1 1 1 0 0 0 [67] 0 0 0 0 0 0 0 1 1 1 0 1 0 0 1 1 0 1 0 1 0 1 1 1 0 0 1 0 0 1 1 1 1 [100] 1 1 1 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 0 0 1 0 0 1 1 0 0 0 1 0 1 0 0 [133] 0 0 0 0 0 0 1 0 1 1 0 1 1 1 0 0 0 0 0 0 0 1 0 0 1 0 0 0 1 0 0 1 0 [166] 0 0 0 1 1 1 0 1 0 1 0 1 0 0 0 1 0 0 0 0 0 0 1 0 1 1 1 0 0 1 1 0 1 [199] 0 0 0 0 1 0 1 0 0 0 0 0 0 1 1 1 1 0 0 0 0 0 1 0 1 1 0 0 0 1 0 0 1 [232] 0 0 0 1 1 0 1 0 0 1 0 0 0 0 1 0 1 0 1 0 1 0 1 0 0 1 0 0 0 0 0 0 1 [265] 0 1 0 0 0 1 0 0 0 1 0 1 1 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 1 1 0 1 [298] 0 1 1 0 0 0 0 1 0 1 0 1 0 0 0 0 0 0 0 1 0 0 0 1 1 0 0 0 1 0 0 1 0 [331] 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 1 0 0 1 0 0 1 0 1 0 0 0 0 1 0 1 1 etc Rather than like this John 0 Frank 1 Ann 0 James 1 Alex 1 Can anyone suggest a script that would achieve this? Thanks Peter
Hi, On Mon, Oct 18, 2010 at 7:37 AM, Peter Francis <peterfrancis at me.com> wrote:> Dear List, > > I have a table i have read into R: > > Name ? ?Yes/No > > John ? ?0 > Frank ? 1 > Ann ? ? ? ? ? ? 0 > James ? 1 > Alex ? ?1 > > etc ?- 800 different times. > > What i want to do is shuffle yes/no and randomly re-assign them to the name.I guess you mean that you have a data.frame object. Let's say this is called "response". You can get 800 permutations of response like so: R> responses <- lapply(1:800, function(x) { x <- response x[,2] <- sample(x[,2]) x }) You can also look at "replicate" and use it in a similar fashion.> I have used sample() and permute(), however there is no way to do this 1000 times. Furthermore, i want to copy the data into a excel spreadsheet in the same order as the data was input so i can build up a distribution of the statistic for each name.It seems like you can do that in R, but if you want to use excel, you can include a call to `write.table` before on x before you return it from the lapply/replicate loop, then have excel import a tab/comma delimited file. There are some R packages that I believe write to excel format directly, but I can't recall what they were ... this might be one: http://www.omegahat.org/RExcelXML/ .. I'm pretty sure there is another, but it escapes me. -steve -- Steve Lianoglou Graduate Student: Computational Systems Biology ?| Memorial Sloan-Kettering Cancer Center ?| Weill Medical College of Cornell University Contact Info: http://cbio.mskcc.org/~lianos/contact
Hi:
One way to permute your sample 1000 times is to use the r*ply() function
from the plyr package. Let s denote your original vector, randomly generated
as follows:
s <- rbinom(800, 0.6) # simulates your Yes/No vector
library(plyr)
u <- raply(1000, sample(s)) # generates a 1000 x 800 matrix
As Steve mentioned, you can also use replicate():
v <- replicate(1000, sample(r)) # generates an 800 x 1000 matrix
# With mydata as your original data frame, the 1000 reshuffles can be
attached with
df <- data.frame(mydata, t(u)) # or df <-
data.frame(mydata, v)
names(df)[3:1002] <- paste('S', 1:1000, sep = '')
write.csv(df, file = 'myExcelFile.csv', quote = FALSE, row.names =
FALSE)
Substitute your input vector for s in the raply() call. [Note that
replicate() is about 10 times faster than raply() for this example.]
The result should be an 800 x 1002 (two original variables + 1000
reshuffles) data frame. Writing to a csv file is a convenient intermediary
between R and Excel, but there are several ways to write data from R to
Excel. See
http://rwiki.sciviews.org/doku.php?id=tips:data-io:ms_windows
HTH,
Dennis
On Mon, Oct 18, 2010 at 4:37 AM, Peter Francis <peterfrancis@me.com>
wrote:
> Dear List,
>
> I have a table i have read into R:
>
> Name Yes/No
>
> John 0
> Frank 1
> Ann 0
> James 1
> Alex 1
>
> etc - 800 different times.
>
> What i want to do is shuffle yes/no and randomly re-assign them to the
> name.
>
> I have used sample() and permute(), however there is no way to do this 1000
> times. Furthermore, i want to copy the data into a excel spreadsheet in the
> same order as the data was input so i can build up a distribution of the
> statistic for each name. When i use shuffle the date gets returned like
this
> -
>
> [1] 1 0 0 1 0 1 0 0 0 1 1 1 1 1 1 1 1 1 0 1 0 0 0 0 0 0 1 0 0 0 1 0 1
> [34] 0 1 0 0 0 0 0 0 0 1 0 1 1 0 1 0 0 0 0 0 0 1 1 0 0 1 0 1 1 1 0 0 0
> [67] 0 0 0 0 0 0 0 1 1 1 0 1 0 0 1 1 0 1 0 1 0 1 1 1 0 0 1 0 0 1 1 1 1
> [100] 1 1 1 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 0 0 1 0 0 1 1 0 0 0 1 0 1 0 0
> [133] 0 0 0 0 0 0 1 0 1 1 0 1 1 1 0 0 0 0 0 0 0 1 0 0 1 0 0 0 1 0 0 1 0
> [166] 0 0 0 1 1 1 0 1 0 1 0 1 0 0 0 1 0 0 0 0 0 0 1 0 1 1 1 0 0 1 1 0 1
> [199] 0 0 0 0 1 0 1 0 0 0 0 0 0 1 1 1 1 0 0 0 0 0 1 0 1 1 0 0 0 1 0 0 1
> [232] 0 0 0 1 1 0 1 0 0 1 0 0 0 0 1 0 1 0 1 0 1 0 1 0 0 1 0 0 0 0 0 0 1
> [265] 0 1 0 0 0 1 0 0 0 1 0 1 1 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 1 1 0 1
> [298] 0 1 1 0 0 0 0 1 0 1 0 1 0 0 0 0 0 0 0 1 0 0 0 1 1 0 0 0 1 0 0 1 0
> [331] 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 1 0 0 1 0 0 1 0 1 0 0 0 0 1 0 1 1
>
> etc
>
> Rather than like this
>
> John 0
> Frank 1
> Ann 0
> James 1
> Alex 1
>
> Can anyone suggest a script that would achieve this?
>
> Thanks
>
> Peter
>
> ______________________________________________
> R-help@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.
>
[[alternative HTML version deleted]]