我试图计算以下文件中包含的省略号的偏心率:
https://www.dropbox.com/s/cqwghg7uvj1xl4d/ellipses.csv?dl=0
它组织为100个点,定义每组100个椭圆。组是A,B,C ......组中的每个椭圆表示在" rep"列。
我以为我可以使用dplyr做到这一点,但我遇到了麻烦,因为我是dplyr的新手。我在post上找到了如何从椭圆点计算偏心率。步骤是:计算椭圆的质心,找到半长轴和半短轴,然后计算偏心率。
我理解如何group_by,但我不知道如何处理汇总部分以完成多步计算。这是计算质心的原因。
ellipse_df <- read.csv("~/Desktop/ellipses.csv")
ellipse_df %>%
group_by(group, rep) %>%
summarise(
mean(iso1),mean(iso2))
显然这并不多。我希望输出保持数据的分组,因此每组获得100个偏心率值。
感谢您的帮助!