[R] add a list name into the list content in a new output
Chel Hee Lee
chl948 at mail.usask.ca
Tue Jan 20 04:44:23 CET 2015
I am just using the first two components of your output.
> op
$bp
chr pos ci.low ci.high lod
c7.loc45 7 47.7 36.71 56.7 6.11
c15.loc8 15 12.0 3.96 22.8 5.29
$hr
chr pos ci.low ci.high lod
c2.loc54 2 59.8 14.8 87.8 4.19
D15MIT184 15 22.8 12.0 36.0 3.15
>
> result <- lapply(names(op), function(x){
+ col1 <- x
+ col2 <- row.names(op[[x]])
+ mat <- op[[x]]
+ row.names(mat) <- NULL
+ rval <- cbind(col1, col2, mat)
+ names(rval) <- c("trait", "locus", names(mat))
+ rval
+ })
>
> result
[[1]]
trait locus chr pos ci.low ci.high lod
1 bp c7.loc45 7 47.7 36.71 56.7 6.11
2 bp c15.loc8 15 12.0 3.96 22.8 5.29
[[2]]
trait locus chr pos ci.low ci.high lod
1 hr c2.loc54 2 59.8 14.8 87.8 4.19
2 hr D15MIT184 15 22.8 12.0 36.0 3.15
>
Is this what you are looking for? I hope this helps.
Chel Hee Lee
On 01/19/2015 02:14 PM, Shuhua Zhan wrote:
> Dear All,
> I'd like to add a list name into the list contents to make a new output. The list is a list of data.frame derived from summary command in Rqtl. I want to add this list name to the data frame with a given column name such as "trait" and output this entire list as a table tab delimited as below.
> Here is the list generated by summary command in Rqtl:
>> summary(out.all, threshold=3, format="tabByCol")
> bp:
> chr pos ci.low ci.high lod
> c7.loc45 7 47.7 36.71 56.7 6.11
> c15.loc8 15 12.0 3.96 22.8 5.29
>
> hr:
> chr pos ci.low ci.high lod
> c2.loc54 2 59.8 14.8 87.8 4.19
> D15MIT184 15 22.8 12.0 36.0 3.15
>
> bw:
> chr pos ci.low ci.high lod
> c15.loc16 15 20 11 30 6.75
>
> heart_wt:
> chr pos ci.low ci.high lod
> c12.loc49 12 51.2 28.2 62.2 3.65
>
>
> The new output I want:
>
> trait locus chr pos ci.low ci.high lod
> bp c7.loc45 7 47.7 36.71 56.7 6.11
> bp c15.loc8 15 12.0 3.96 22.8 5.29
>
> trait locus chr pos ci.low ci.high lod
> hr c2.loc54 2 59.8 14.8 87.8 4.19
> hr D15MIT184 15 22.8 12.0 36.0 3.15
>
> trait locus chr pos ci.low ci.high lod
> bw c15.loc16 15 20 11 30 6.75
>
> heart_wt locus chr pos ci.low ci.high lod
> trait c12.loc49 12 51.2 28.2 62.2 3.65
>
> I appreciate for your help for any suggestions!!
> Joshua
>
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> 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.
>
More information about the R-help
mailing list