|
1 |
| -globalVariables(c("cdslen", "end", "start", 'tga', "tlen", ".")) |
| 1 | +globalVariables(c("cdslen", "end", "start", 'tga', "tlen", ".","gid")) |
2 | 2 |
|
3 | 3 | # ==============================================================================
|
4 | 4 | # setGeneric
|
@@ -236,7 +236,8 @@ setMethod("filterRepTrans",
|
236 | 236 | }
|
237 | 237 |
|
238 | 238 | # return data.frame
|
239 |
| - tranLength <- data.frame(gid = unique(tg$gene_id), |
| 239 | + tranLength <- data.frame(gname = unique(tg$gene_name), |
| 240 | + gid = unique(tg$gene_id), |
240 | 241 | tid = tida[x],
|
241 | 242 | tlen = sum(tmp$width),
|
242 | 243 | cdslen = cdslen)
|
@@ -344,13 +345,16 @@ setMethod("getLongName",
|
344 | 345 | if("5UTR" %in% ginfo.s$type | "five_prime_utr" %in% ginfo.s$type){
|
345 | 346 | cds.sp <- cds.len + cds.st - 1
|
346 | 347 | }else{
|
347 |
| - cds.sp <- cds.len + cds.st |
| 348 | + if("CDS" %in% ginfo.s$type){ |
| 349 | + cds.sp <- cds.len + cds.st - 1 |
| 350 | + }else{ |
| 351 | + cds.sp <- cds.len + cds.st |
| 352 | + } |
348 | 353 | }
|
349 | 354 |
|
350 |
| - |
351 | 355 | # name
|
352 |
| - tname <- paste(gname,tid,cds.st,cds.sp,tlen,sep = sep) |
353 |
| - res <- data.frame(gid = gname,tid = tid,tname = tname) |
| 356 | + tname <- paste(gname,gid,tid,cds.st,cds.sp,tlen,sep = sep) |
| 357 | + res <- data.frame(gname = gname,gid = gid,tid = tid,tname = tname) |
354 | 358 | return(res)
|
355 | 359 | }) -> resName
|
356 | 360 | # }) %>% do.call("rbind",.) %>% data.frame() -> resName
|
|
0 commit comments