首页 文章
  • 0 votes
     answers
     views

    R data.table合并/完全外连接与na.fill / nomatch基于公式

    我需要的是以有效的方式执行完全外部联接与某种智能na.fill / nomatch . 我已经使用循环完成了它但我想使用矩阵代数或data.table操作来加速这个过程 . 以下数据是股票未结订单信息的样本,在要求未结订单和出价未结订单的数据集之间执行完全外部联接 . 要求数据集,B是出价 . 两个数据集都存储原子订单及其累积总和 . 任务是按累计值匹配所有询价单和出价单,反之亦然 . 填充示例数...
  • -3 votes
     answers
     views

    不同动物园对象之间的相关矩阵

    我有几个动物园对象的文件看起来像这样(每个文件在不同的日期开始和结束): code pp 1942-06-01 4016 0 1942-06-02 4016 NA 1942-06-03 4016 0 1942-06-04 4016 0 1942-06-05 NA 0 1942-06-06 NA 0 我想在9月,10月和11月的所有文件的pp之间做一个相...
  • 3 votes
     answers
     views

    XTS to.weekly返回不同的每周 endpoints

    我有一个问题, xts 中的 endpoints() 函数(以及使用 endpoints 的 to.weekly 函数)有时会在星期五作为周末返回,有时会返回星期一 . 我的数据集名为 sp2 . > head(sp2) [,1] 2012-01-09 1.78 2012-01-10 1.78 2012-01-11 1.77 2012-01-12 1.80 2012-...
  • 2 votes
     answers
     views

    R:按组插入NA

    我想在数据帧的变量中执行线性插值,其考虑:1)两点之间的时间差,2)获取数据的时刻和3)用于测量变量的个体 . 例如,在下一个数据帧中: df <- data.frame(time=c(1,2,3,4,5,6,7,1,2,3), Individuals=c(1,1,1,1,1,1,1,2,2,2), Value=c(1, 2, 3, NA, 5...
  • 0 votes
     answers
     views

    预测xts对象不生成未来预测点的日期

    我有一个数据框,其中一列代表我的每日观察值,一个日期列使用as.Dates函数格式化为“2014-10-01” . 我正在创建一个xts对象,执行以下操作: xtsObject <- as.xts(vector$values, order.by = vector$dates) 然后我将xtsObject传递给ets预测函数,并在结果对象上调用预测来预测接下来的三个时间步骤,如下所示: fi...
  • 0 votes
     answers
     views

    按组计算滚动加权总和

    我有一个在过去12小时内加权交易总和的工作示例 . 现在,我添加了一个帐户列,并希望按组分别计算此加权金额 . 代码将按如下所示运行 . 取消注释以 # account 开头的行,将 account 列添加到 df . How can I modify the second to last line of code such that it computes rollapplyr separ...
  • 0 votes
     answers
     views

    使用as.yearmon(动物园)将字符值转换为月度数据时出错

    当月份的值以“Jan”或“January”等文本形式出现时会出现问题 . 例如,以下行 as.yearmon("Jan-2000", "%b-%Y") 返回错误: Error in charToDate(x) : > character string is not in a standard unambiguous format 这里出了什么问题?...
  • 1 votes
     answers
     views

    在R中填充缺失值的时间序列数据

    我正在尝试将我的面板数据中的年度值扩展到年度季度值 . 这是每个季度重复的年度值 . 例如,我希望获得2000Q1,2000Q2,2000Q3,2000Q4,2001Q1,......,2001Q4的年度收入的重复值 . 因此数据框将是id,年度季度,收入 . 我使用两步法,但有一些问题需要处理 . 如果缺少季度起始值,那么我将需要季度缺失(NA) . 情况1: annual_data <-...
  • 6 votes
     answers
     views

    用于绘制时间序列数据的x轴刻度的日期格式

    数据文件的日期是格式,即1975M1,1975M2,... 2011M12,用于时间序列数据 . 当使用R绘制此数据时,我希望x轴显示刻度轴上的月份 . 对于要正确读取的日期,我已经尝试将M替换为 - 以获得%Y-%m格式,但这对于来自hydroTSM包的drawTimeAxis似乎不太好,这可能需要%Y-%M-%d格式 . 它给出了刻度尺维度的维数不正确的错误 . 另一种解析和格式化数据的方法,...
  • 0 votes
     answers
     views

    R read.zoo()可以将's accept object'的值输入作为参数 . [R

    我创建了一个名为data1的data.frame,并通过read.zoo函数读取它,如下所示 . data2<-read.zoo(data1) 现在,我使用以下代码为现有data.frames的整个列表执行此操作 . names_of_dataframes <- ls.str(mode = "list") data2<-read.zoo(name_of_da...
  • 2 votes
     answers
     views

    从每日时间序列到R xts对象中的每周时间序列

    我正在使用zoo和xts包来分析财务数据 . ts套餐不是很合适,因为金融系列有周末没有数据 . 我在xts包中读到了apply function availbale apply.daily(x, FUN, ...) apply.weekly(x, FUN, ...) apply.monthly(x, FUN, ...) apply.quarterly(x, FUN, ...) apply.ye...
  • 2 votes
     answers
     views

    使用R将连续时间序列数据转换为每日每小时表示

    我把时间序列数据用xts表示为 library(xts) xtime <-timeBasedSeq('2015-01-01/2015-01-30 23') df <- xts(rnorm(length(xtime),30,4),xtime) 现在我想计算不同日期之间的共同化,因此我希望以矩阵形式表示 df : 为此,我用了 p_mat= split(df,f="day...
  • 0 votes
     answers
     views

    每周在R中分割时间序列

    我想每周在 R 分割 xts/zoo 时间序列 . 时区设置为"Asia/Kolkata" Sys.setenv(TZ="Asia/Kolkata") library(xts) seqs<- seq(as.POSIXct("2016-01-01"),as.POSIXct("2016-01-30"), by = &...
  • 0 votes
     answers
     views

    如何将每日 Value 汇总限制为月度和年度总计?

    我有一个动物园对象(stn)的每日数据的长记录(日期Vs下雨) . 我用过hydoTSM包 (m < - daily2monthly(stn,FUN = sum,na.rm = TRUE)) 这工作正常,但我想只允许聚合几个月,至少20个每日记录 . 对于年度,应提供至少10个月的记录 . 任何人都可以建议一个更好的方法吗?
  • 2 votes
     answers
     views

    缺少 Value - 有马模型

    我有关于产品销售的每日时间序列,我的系列从2016年1月1日至2017年8月31日开始 . 考虑到它是一个为期六天的星期(我的星期一星期一开始,星期六结束)星期日没有数据,我知道在运行Arima模型之前我需要先填写缺失值 . 这是我需要帮助的地方:我已经读过我可以用 na.approx 或 NA 填充缺失值,但我不知道该怎么做 . 你可以在这里看到我的系列: https://drive.googl...
  • 0 votes
     answers
     views

    获取具有每日频率的时间序列的每周季节图表

    如何在R中分解每日频率的时间序列?我在一年中每天都有一个网站的访问者数量 . 我想展示每周一季的图表 . my_data = read.csv("time series test.csv", header = TRUE) my_zoo_ts = read.zoo("time series test.csv", sep = ",", fo...
  • 16 votes
     answers
     views

    R:合并两个不规则的时间序列

    我有两个多变量时间序列x和y,两者都覆盖了大致相同的时间范围(一个在另一个之前两年开始,但它们在同一天结束) . 两个系列都以日期列旁边的空列形式缺少观察结果,并且在某种意义上,其中一个系列具有在另一个系列中找不到的几个日期,反之亦然 . 我想创建一个数据框(或类似),其中列列出了在x或y中找到的所有日期,没有重复日期 . 对于每个日期(行),我想在y的观察值旁边水平叠加来自x的观测值,其中NA填...
  • 3 votes
     answers
     views

    每小时Holt-Winter时间序列预测(预测)

    我对R包很新,我正在处理时间序列 . 我必须 Build 一个预测模型来预测未来的点击 . 预测的时间间隔需要每小时一次 . 我的示例时间序列: DateTime Clicks (06/23/13 00:00:00) 757 (06/23/13 01:00:00) 714 (06/23/13 02:00:00) 776 (06/23/13 03:00:00) 87...
  • 0 votes
     answers
     views

    使用R在不同的面板中绘制多个时间序列

    我试图绘制9个不同面板的9个站(A到I)的时间序列月平均温度数据 . 但我使用动物园来保持月份格式 . 我的简单数据集如下所示 . +--------+------+------+------+------+ | Time | A | B | C | D | +--------+------+------+------+------+ | Jan-84 | 28.2 |...
  • 0 votes
     answers
     views

    用dplyr和lsfit滚动回归

    我正试图用 dplyr 进行滚动回归 . 我正在使用包 zoo 和 lsfit 中的 rollapplyr 作为我尝试过的'm only interested in the beta of the regression. Here': library(dplyr); library(zoo) df1 = expand.grid(site = seq(10), ...
  • 1 votes
     answers
     views

    每天给定间隔的最大斜率

    我有一组时间序列数据,其中地面温度每隔10分钟从三个不同的位置测量多天(实际上是2年的数据) . 我感兴趣的是计算每个站点每天60分钟间隔的最大斜率(温度升高速率) . 所以基本上我想每天工作10分钟,每次60分钟,并计算每个窗口的坡度,然后确定最大坡度以及白天发生的时间 . 然后,我想将此函数应用于数据集中的每一天 . 日期/时间采用以下格式(%m /%d /%y%H:%M) . 我正在想象使用...
  • 3 votes
     answers
     views

    用dplyr滚动回归

    我有一个“日期”,“公司”和“返回”的数据框,可通过以下代码重现: library(dplyr) n.dates <- 60 n.stocks <- 2 date <- seq(as.Date("2011-07-01"), by=1, len=n.dates) symbol <- replicate(n.stocks, paste0(sample(LET...
  • 1 votes
     answers
     views

    R:rollapplyr和lm因子错误:rollapplyr是否更改了变量类?

    这个问题 Build 在前一个问题的基础之上,这个问题对我很有帮助 . R: Grouped rolling window linear regression with rollapply and ddply 难道你不知道代码在扩展到真实数据而不是示例数据时不能正常工作吗? 我有一个有点大的数据集,具有以下特征 . str(T0_satData_reduced) 'data.frame': 4...
  • 0 votes
     answers
     views

    R中多个股票的滚动Beta回归

    我正在尝试计算宽度为过去12个月的多个股票的滚动Beta回归 . 我有following dataset 好像: 我正在搜索很多帖子,但不知怎的,我没有让它为我的数据框架工作 . func1 <- . %>% { roll_regres.fit(x = cbind(1, .$MKT_ex), y = .$r_rf, width = 12L)$coefs }...
  • 4 votes
     answers
     views

    在过去365天窗口中执行运行总计的有效方法

    这就是我的数据框架: 库(data.table) df <- fread(' Name EventType Date SalesAmount RunningTotal Runningtotal(prior365Days) John Email 1/1/2014 0 0 ...
  • 0 votes
     answers
     views

    R中时间序列相关系数的运行窗口

    我想用R中的运行窗口计算12个变量的相关系数 . 我的数据存储在动物园对象中,%m . %d . %Y%H:%M:%S索引,12个变量中的每一个都有1343个观察值 . 我不知道我将要使用的窗口大小,但如果需要我可以更改它 . @Joshua Ulrich发布了here如何使用rollapplyr计算滚动相关性,但这个例子只有两个变量 . 由于我的R经验有限,我不确定如何使用其中一个应用族函数来运...
  • 1 votes
     answers
     views

    R仅删除具有NA的组

    我有一个类似于以下结构生成的数据帧: library(dplyr) df1 <- expand.grid(region = c("USA", "EUR", "World"), time = c(2000, 2005, 2010, 2015, 2020), ...
  • 1 votes
     answers
     views

    合并两个不规则的动物园时间序列正在弄乱结构

    我正在使用包含数千行的贸易数据集 . 每条记录都有一个基于符号和日期的唯一键 . 给定符号的交易记录是不规则的,因此使用动物园将是自然的选择 . 我需要使用lag和merge来创建一个新的数据集 . 但是,我不知道如何在动物园中设置多列索引以便使用滞后函数 . 下面是一个示例数据集和预期输出 . df = data.frame( dt = as.Date(c("2015-01-0...
  • 1 votes
     answers
     views

    使用LCOF和NOCB方法填写时间序列缺口,但确认时间序列中断

    There are edits to this post at the end. 我有一个关于一群人的每日膳食记录的大型数据集 . 每个人都有随机丢失的数据 . 这是一个人的例子(我最终将这个解决方案推广到人口): > str(final_daily) 'data.frame': 387 obs. of 10 variables: $ Date : chr ...
  • 2 votes
     answers
     views

    R根据时间间隔通过线性增加来计算NA

    PROBLEM 我需要从最后一次观察到的值开始,每52周间隔用最后一次观察到的非NA值1来估算NA 's in my data frame that comes from a repeated measures study. On this particular outcome, I need to impute the NA' . EXAMPLE 包含目标插补目标的示例数据框 . df <...

热门问题