how to create geno.rna object in https://rpubs.com/akdes/673120

[BioAmelie]

trafficstars

@akdess Hi ,

Can you tell me how to create geno.rna object in https://rpubs.com/akdes/673120. Best, minfang

[Junedays]

I have exactly the same question. The protocol did not show how to create the geno.rna

[Junedays]

I guess maybe there is a typo in the previous line of the codes.

rna <- GRanges(seqnames = Rle(gsub("q", "", gsub("p", "", all.summary$Chromosome))), IRanges(all.summary$Start, all.summary$End))

should be

geno.rna <- GRanges(seqnames = Rle(gsub("q", "", gsub("p", "", all.summary$Chromosome))), IRanges(all.summary$Start, all.summary$End))

But I am not quite sure.