我正在处理一个不平衡、不规则间隔的横截面时间序列。我的目标是获得"数量"向量的滞后移动平均向量，由"主题"分割。
换句话说，假设已对 Subject_1 观察到以下数量：
[1,2,3,4,5]。我首先需要将它滞后 1，得到 [NA,1,2,3,4]。
然后我需要取 3 阶的移动平均线，得到 [NA,NA,NA,(3 2 1)/3,(4 3 2)/3]
以上所有科目都需要完成。
2 3 4 5 6 7 8 9 10 11 | panel <- data.frame( as.factor(sort(rep(1:6,5))), rep(1:5,6), rnorm(30) ) colnames(panel)<- c("Subject","Day","Quantity") #Make panel DF unbalanced panelUNB <- subset(panel,as.numeric(Subject)!= Day) panelUNB <- panelUNB[-c(15,16),] |
如果面板是平衡的，我将首先使用包
然后我会像这样使用包
中的函数
2 | x,3,align="right",fill=NA,na.rm=TRUE)) |
这将在应用于平衡的"面板"DF 时产生正确的结果。
问题在于
StackExchange 上有一个带有 data.table 的解决方案，它暗示了我的问题的解决方案：Producing a rolling average of an unbalanced panel data set
也许可以修改此解决方案以生成固定长度的移动平均线，而不是"滚动累积平均线"。
这会给你想要的结果吗？
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 | library(zoo) # create time series where each subject have an observation at each time step d1 <- data.frame(subject = rep(letters[1:4], each = 5), day = rep(1:5, 4), quantity = sample(x = 1:4, size = 20, replace = TRUE)) d1 # select some random observations d2 <- d1[sample(x = seq_len(nrow(d1)), size = 15), ] d2 # reshape to wide format with dcast # -> 'automatic' extension from irregular to regular series for each subject, # _given_ that all time steps are represented. # Alternative method below more explicit # fill for structural missings defaults to NA d3 <- dcast(d2, day ~ subject, value.var ="quantity") d3 # convert to zoo time series z1 <- zoo(x = d3[ , -1], order.by = d3$day) ################################ # alternative method to extend time series # time steps to include are given explicitly # create a zero-dimensional zoo series z0 <- zoo(, min(d1$day):max(d1$day)) # extend z1 to contain the same time indices as z0 z1 <- merge(z1, z0) ################################ # lag, defaults to one unit z2 <- lag(x = z1) z2 # calculate rolling mean with window width 3 rollmeanr(x = z2, k = 3) # Handling of NAs: # from ?rollmean: #"The default method of rollmean does not handle inputs that contain NAs. # In such cases, use rollapply instead.": rollapplyr(data = z2, width = 3, FUN = mean, na.rm = TRUE) |
所以，要回答我自己的问题，一种方法是通过 split-lapply(rollingaverage)-unlist:
2 3 4 | Temp <- lapply(Temp, FUN=function (x) rollapplyr( x,2,align="right",fill=NA,na.rm=TRUE, FUN=mean)) QuantityMA <-unlist(Temp) |
然后必须将"QuantityMA"向量添加回主"panelUNB"框架。似乎正在工作。可以使用 ddply 在不平衡面板上完成滞后。
如果有人有其他可能更优雅的解决方案，欢迎分享。
如何匹配未被反斜杠转义的平衡定界符对(其本身未被反斜杠转义)(无需考虑嵌套)？例如对于反引号，我试过了，但是转义的反引号没有像转义那样工作。regex=/(?!$1:"how\\"#expected"how\\`are"上面的正则表达式不考虑由反斜杠转义并位于反引号前面的反斜杠，但我愿意考虑。StackOverflow如何做到这一点？这样做的目的并不复杂。我有文档文本，其中包括内联代码的反引号，就像StackOverflow一样，我想在HTML文件中显示它，内联代码用一些spanMaterial装饰。不会有嵌套，但转义反引号或转义反斜杠可能出现在任何地方。
我想设置一个默认日期，例如实际日期，我该如何设置？还有如何在组合框中设置默认值顺便问一下，date_field_tag和date_field之间有什么区别？ 最佳答案 试试这个:将默认日期作为第二个参数传递。youcorrectlysetthedefaultvalueofcomboboxasshowninyourquestion. 关于ruby-on-rails-date_field_tag，如何设置默认日期？[rails上的ruby]，我们在StackOverflow上找到一个类似的问
我将我的Rails应用程序部署到OpenShift，它运行良好，但我无法在生产服务器上运行“Rails控制台”。它给了我这个错误。我该如何解决这个问题？我尝试更新rubygems，但它也给出了权限被拒绝的错误，我也无法做到。railsc错误:Warning:You'reusingRubygems1.8.24withSpring.UpgradetoatleastRubygems2.1.0andrun`gempristine--all`forbetterstartupperformance./opt/rh/ruby193/root/usr/share/rubygems/rubygems
我正在尝试从Postgresql表(table1)中获取数据，该表由另一个相关表(property)的字段(table2)过滤。在纯SQL中，我会这样编写查询:SELECT*FROMtable1JOINtable2USING(table2_id)WHEREtable2.propertyLIKE'query%'这工作正常:scope:my_scope,->(query){includes(:table2).where("table2.property":query)}但我真正需要的是使用LIKE运算符进行过滤，而不是严格相等。然而，这是行不通的:scope:my_scope,->(que
有没有办法在这个简单的get方法中添加超时选项？我正在使用法拉第3.3。Faraday.get(url)四处寻找，我只能先发起连接后应用超时选项，然后应用超时选项。或者有什么简单的方法？这就是我现在正在做的:conn=Faraday.newresponse=conn.getdo|req|req.urlurlreq.options.timeout=2#2secondsend 最佳答案 试试这个:conn=Faraday.newdo|conn|conn.options.timeout=20endresponse=conn.get(url
“输出”是一个序列化的OpenStruct。定义标题try(:output).try(:data).try(:title)结束什么会更好？:) 最佳答案 或者只是这样:deftitleoutput.data.titlerescuenilend 关于ruby-on-rails-更好的替代方法try(:output).try(:data).try(:name)?，我们在StackOverflow上找到一个类似的问题： https://stackoverflow.c
我需要一个表，其中行实际上是2行表，一个嵌套表是..我怎样才能在Prawn中做到这一点？也许我需要延期..但哪一个？ 最佳答案 现在支持子表:Prawn::Document.generate("subtable.pdf")do|pdf|subtable=pdf.make_table([["sub"],["table"]])pdf.table([[subtable,"original"]])end 关于ruby-on-rails-PrawnPDF:Ineedtogeneratenested
我想在Ruby中创建一个用于开发目的的极其简单的Web服务器(不，不想使用现成的解决方案)。代码如下:#!/usr/bin/rubyrequire'socket'server=TCPServer.new('127.0.0.1',8080)whileconnection=server.acceptheaders=[]length=0whileline=connection.getsheaders想法是从命令行运行这个脚本，提供另一个脚本，它将在其标准输入上获取请求，并在其标准输出上返回完整的响应。到目前为止一切顺利，但事实证明这真的很脆弱，因为它在第二个请求上中断并出现错误:/usr/b
我意识到这可能是一个非常基本的问题，但我现在已经花了几天时间回过头来解决这个问题，但出于某种原因，Google就是没有帮助我。(我认为部分问题在于我是一个初学者，我不知道该问什么......)我也看过O'Reilly的RubyCookbook和RailsAPI，但我仍然停留在这个问题上.我找到了一些关于多态关系的信息，但它似乎不是我需要的(尽管如果我错了请告诉我)。我正在尝试调整MichaelHartl'stutorial创建一个包含用户、文章和评论的博客应用程序(不使用脚手架)。我希望评论既属于用户又属于文章。我的主要问题是:我不知道如何将当前文章的ID放入评论Controller。
我的工作要求我为某些测试自动生成电子邮件。我一直在四处寻找，但未能找到可以快速实现的合理解决方案。它需要在outlook而不是其他邮件服务器中，因为我们有一些奇怪的身份验证规则，我们需要保存草稿而不是仅仅发送邮件的选项。显然win32ole可以做到这一点，但我找不到任何相当简单的例子。 最佳答案 假设存储了Outlook凭据并且您设置为自动登录到Outlook，WIN32OLE可以很好地完成此操作:require'win32ole'outlook=WIN32OLE.new('Outlook.Application')message=
