0

我有以下数据集:

require(xts)
days<-seq(as.Date("2016-05-01", format="%Y-%m-%d"), as.Date("2016-07-31", format="%Y-%m-%d" ), by= 1)
df.days<-xts(x = data.frame(Empty=rep(NA, length(days))), order.by=days)
df.norm <- rnorm(10000, 40, 5)
df.norm <- df.norm[df.norm > 30 & df.norm < 50]
SAMPLE <-sample(df.norm, length(df.days[.indexwday(df.days)==1]), replace=FALSE) 
df.weeks<-xts(SAMPLE, order.by = index(df.days[.indexwday(df.days)==1]))
df.data<-merge.xts(df.weeks,df.days)
df.data<-df.data[,1]
colnames(df.data)[1]<-"rate1"
df.data<-as.data.frame(df.data)
rownames(df.data)<-NULL
df.data<-data.frame(times=rownames(df.data), coredata(df.data))

rate1 是一个与时间相关的函数,在几天内趋向于接近 40 的值。

我想知道 rate1 在什么日期等于 40。

做我需要的一种可能性是密集填充 rate1 并减去我的目标值 (40)。一个点越接近零,它就越接近目标值。并包括一个小的“噪声”缓冲区,接受零范围内的所有点。

绝对值(f[t]-\Theta)\leq \epsilon

如何估计 R 中 rate1 的时间值等于 40?另外请建议R中是否有比我使用的更合适的方法?欢迎所有建议。

4

0 回答 0