bioinformatics-ptp/detectRUNS

return mean values and standard error/deviation

Opened this issue · 2 comments

The data obtained from the summaryRuns option is great, however, for the "count and count_chr" tables it would be very useful to obtain the mean ROH values ( plus standard dev/error) according to the number of samples per population. That way we can really compare them. Right now we get a total number which could be biased when the pops size is different. Same for the ROH length cause we get the mean but not the standard deviation of the values.

I think it would be straightforward to achieve for someone with good R skills not like me hah but I think that by adding these options to the ddply function "SE = sd(N) / sqrt((length(N))), SD = sd(N))" ??

Anyway it is just a suggestion but I think it would be very much appreciated.
Cheers!

bunop commented

Dear @scisea ,
Thank you for your interest in detectRUNs!

for the "count and count_chr" tables it would be very useful to obtain the mean ROH values ( plus standard dev/error) according to the number of samples per population.

This seems reasonable to me. @filippob , what do you think? could this feature be added in the next release?