我可以在R中以时间作为因变量(y)进行ANOVA检验吗?

Question

我想测试一下蝌蚪的起源是否会影响它的蜕变时间，也就是说，如果一只蝌蚪的起源在北方的一个池塘里，它的蜕变会比一只蝌蚪的起源在南方的蜕变早吗?为了做到这一点，我希望使用方差分析测试。我的代码是这样的:

DATE_OF_METAMORPHOSIS -> dmy_hms(DATE_OF_METAMORPHOSIS)
one.way <- aov(DATE_OF_METAMORPHOSIS ~ POPULATION, data = info_table)

它确实计算了一些东西，但我不确定R如何处理DATE_OF_METAMORPHOSIS变量中的日期。如果有人知道我的代码是否做了我用文字描述的，或者知道R是否知道如何使用日期作为延续变量以及如何做到这一点，我将很乐意帮助!tnx !

示例:
info_table

ID | POPULATION | DATE_OF_METAMORPHOSIS|
---|------------|----------------------|
1 | 1 | 7/19/2021  12:01:00 AM|
2 | 2 | 7/29/2021  12:01:00 AM|
3 | 3 | 8/1/2021  12:01:00 AM|
4 | 1 | 8/4/2021  12:01:00 AM|
5 | 2 | 5/16/2021  12:01:00 AM|
6 | 3 | 5/14/2021  12:01:00 AM|

Answer 1

这两个错误是:

1. datetime列的格式为mdy hms，而不是ymd hms;
2. 没有POPULATION列，名称为origin。

纠正这两个错误，这是一个简单的问题运行aov。

x <- '
ID | origin | DATE_OF_METAMORPHOSIS|
1 | 1 | 7/19/2021  12:01:00 AM|
2 | 2 | 7/29/2021  12:01:00 AM|
3 | 3 | 8/1/2021  12:01:00 AM|
4 | 1 | 8/4/2021  12:01:00 AM|
5 | 2 | 5/16/2021  12:01:00 AM|
6 | 3 | 5/14/2021  12:01:00 AM|'
info_table <- read.table(textConnection(x), header = TRUE, sep = "|")[-4]
info_table$DATE_OF_METAMORPHOSIS <- lubridate::mdy_hms(info_table$DATE_OF_METAMORPHOSIS)
one.way <- aov(DATE_OF_METAMORPHOSIS ~ origin, data = info_table)
summary(one.way)
#>             Df    Sum Sq   Mean Sq F value Pr(>F)
#> origin       1 8.885e+12 8.885e+12   0.743  0.437
#> Residuals    4 4.782e+13 1.196e+13

^{在2022-05-24由reprex包(v2.0.1)创建}