vdjtools icon indicating copy to clipboard operation
vdjtools copied to clipboard

Shared clonotype abundance plot labels up-side-down.

Open suntaosimon opened this issue 7 years ago • 4 comments

as figures below,

test paired strict table collapsed

NonOverlapping should be in grey which is on the top.

looks like it's right in the document sample.

image

suntaosimon avatar Sep 11 '17 20:09 suntaosimon

Hi, can you confirm that your clonotype table have frequencies that 1) sum to 1, 2) are strictly positive?

mikessh avatar Sep 12 '17 06:09 mikessh

Yes. They all sum to 1 and positive. image

I have tried many data and they all have the same problem in Shared clonotype abundance plot.

Modified the following parts in "intersect_pair_area_v2.r" that seems help,


#dd <- dd[with(dd, order(cdr3nt, levels(dd$cdr3aa))), ]
df.0$cum <- rev(cumsum(rev(df.0$expr)))
df.0$cum <- (c(df.0$cum[-1],0) + df.0$cum) / 2
df.1$cum <- rev(cumsum(rev(df.1$expr)))
df.1$cum <- (c(df.1$cum[-1],0) + df.1$cum) / 2

img4654

suntaosimon avatar Sep 12 '17 14:09 suntaosimon

The very same thing happens to me.

tiagobrc avatar Jan 18 '18 22:01 tiagobrc

830new.paired.strict.table.collapsed.pdf The very same thing happens to me, I Modified the following parts in "intersect_pair_area.r"

#dd <- dd[with(dd, order(cdr3nt, levels(dd$cdr3aa))), ] df.0$cum <- rev(cumsum(rev(df.0$expr))) df.0$cum <- (c(df.0$cum[-1],0) + df.0$cum) / 2 df.1$cum <- rev(cumsum(rev(df.1$expr))) df.1$cum <- (c(df.1$cum[-1],0) + df.1$cum) / 2

But the labels don't match perfect, how to fix it.

haojw2815 avatar Dec 25 '23 08:12 haojw2815