Stacks Image 20108

For Academics

Merge data

You will often want to combine different datasets with a column in common. This is easy to do in R.

For example the mapping of sequencing reads of different tissues to contigs generated from all the reads.

It could also be different information on the same individuals.

In the sequencing example, the mapping datasets may not fully overlap. For example contig 1 may be represented in tissue 1, but not tissue 2, and hence dataset1 but not dataset2.

We would like to merge the datasets based on their common column (contig). The resulting dataset can be used for statistical analysis, for example for differences in the read number per tissue.

This is straightforward in R using the merge command.

First load the datasets:

tissueData1<-read.table("/Users/username/Desktop/tissueData1", header=T)
str(tissueData1
tissueData2<-read.table("/Users/username/Desktop/tissueData2.txt", header=T)
str(tissueData2)
tissueData3<-read.table("/Users/username/Desktop/tissueData3.txt", header=T)
str(tissueData3)
tissueData4<-read.table("/Users/username/Desktop/tissueData4.txt", header=T)
str(tissueData4)

Now merge by adding one dataset at a time. We merge based on the total_contig and tissue_contig columns, which have overlapping values (i.e. contig1, contig2, etc).

merged1=merge(tissueData1, tissueData2, by.x="tissue1_contig", by.y="tissue2_contig", all=T )
merged2=merge(merged1, tissueData3, by.x="total_contig", by.y="tissue3_contig", all=T )
mergedData=merge(merged3, tissueData4, by.x="total_contig", by.y="tissue4_contig", all=T )

Note: all=T retains the data from rows that did not have a common value in the columns used for merging. In this case NAs are automatically added to the values in the merged row.

The merge command can do more tricks, to learn more about it:

?merge
Previous Post 5 / 10 Post

Tag:

Sex chromosome papers RSS


The genome sequence of the sycamore periphyllus aphid, Periphyllus acericola (Walker, 1848)
Link

The genome sequence of the Straw Grass-veneer moth, Agriphila straminella (Denis & Schiffermüller), 1775
Link

The genome sequence of the Golden-tabbed robberfly, Eutolmus rufibarbis (Meigen, 1820)
Link