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!
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!