[英]Using dates in R for Theil-Sen
I am trying to use dates as my X variable in a Theil-Sen slope estimation and I am having difficulty using the R package zyp 我试图在Theil-Sen斜率估算中将日期用作我的X变量,并且在使用R包zyp时遇到困难
library(zyp)
myDates <- as.Date(c("2009-11-24","2009-12-03","2010-01-19","2010-02-18","2010-03-04"))
myData <- c(10.17,10.91,11.72,12.5,13.43)
a <- c("a","b","c","d","e")
df = data.frame(a,myData,myDates)
zyp.sen(myDates~myData, df)
And then I get the following error: 然后我得到以下错误:
Error in `+.Date`(y - slope * x, intercept) :
binary + is not defined for "Date" objects
I am also trying to use zyp.trend.dataframe(df, 1)
and I get essentially nothing: 我也在尝试使用zyp.trend.dataframe(df, 1)
,但我什么也没得到:
a lbound trend trendp ubound tau sig nruns autocor valid_frac linear intercept
1 a NA NA NA NA NA NA NA NA NA NA NA
2 b NA NA NA NA NA NA NA NA NA NA NA
3 c NA NA NA NA NA NA NA NA NA NA NA
4 d NA NA NA NA NA NA NA NA NA NA NA
5 e NA NA NA NA NA NA NA NA NA NA NA
I am assuming that I am accessing zyp incorrectly. 我假设我不正确地访问zyp。 Does anyone know what I am doing wrong ? 有人知道我在做什么错吗?
TIA TIA
Converting Dates to their underlying integer values should succeed: 将Date转换为它们的基础整数值应该成功:
myDates <- as.integer(
as.Date(c("2009-11-24","2009-12-03","2010-01-19","2010-02-18","2010-03-04")))
myData <- c(10.17,10.91,11.72,12.5,13.43)
a <- c("a","b","c","d","e")
df = data.frame(a,myData,myDates)
zyp.sen(myDates~myData, df)
Call:
NULL
Coefficients:
Intercept myData
14204.66 36.12
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.