Mercurial > repos > davidvanzessen > clonal_sequences_in_paired_samples
comparison RScript.r @ 55:43ce3ebfc9b5 draft
Uploaded
| author | davidvanzessen |
|---|---|
| date | Tue, 13 Oct 2015 09:20:16 -0400 |
| parents | ef5ffc8a6e37 |
| children | c89d9e89423b |
comparison
equal
deleted
inserted
replaced
| 54:ef5ffc8a6e37 | 55:43ce3ebfc9b5 |
|---|---|
| 231 second.rows.v = table(second.rows$V_Segment_Major_Gene) | 231 second.rows.v = table(second.rows$V_Segment_Major_Gene) |
| 232 second.rows.v = names(second.rows.v[which.max(second.rows.v)]) | 232 second.rows.v = names(second.rows.v[which.max(second.rows.v)]) |
| 233 second.rows.j = table(second.rows$J_Segment_Major_Gene) | 233 second.rows.j = table(second.rows$J_Segment_Major_Gene) |
| 234 second.rows.j = names(second.rows.j[which.max(second.rows.j)]) | 234 second.rows.j = names(second.rows.j[which.max(second.rows.j)]) |
| 235 | 235 |
| 236 second.sum = data.frame(merge = paste(first.clone.sequence, "*", sep=""), | 236 second.sum = data.frame(merge = first.clone.sequence, |
| 237 Patient = patient, | 237 Patient = patient, |
| 238 Receptor = second.rows[1,"Receptor"], | 238 Receptor = second.rows[1,"Receptor"], |
| 239 Sample = second.rows[1,"Sample"], | 239 Sample = second.rows[1,"Sample"], |
| 240 Cell_Count = second.rows[1,"Cell_Count"], | 240 Cell_Count = second.rows[1,"Cell_Count"], |
| 241 Clone_Molecule_Count_From_Spikes = sum(second.rows$Clone_Molecule_Count_From_Spikes), | 241 Clone_Molecule_Count_From_Spikes = sum(second.rows$Clone_Molecule_Count_From_Spikes), |