我在R Studio中手动运行时有一个正常运行的脚本,但是当我通过包装器从另一个程序运行它时它不起作用。
我在调试输出中收到此信息:
[912] Error in xj[i] : only 0's may be mixed with negative subscripts
[912] Calls: GetTopN -> cor -> is.data.frame -> [ -> [.data.frame
如果我在收到错误之前保存图像然后在R Studio中加载它,那么当我执行GetTopN(10)
时会出现同样的错误。但是,如果我在R Studio中重新运行语句actionlist<- sqlQuery(channel,al_string)
然后执行GetTopN(10)
,一切正常。
我甚至试图在关键调用之前将图像保存在R Studio中,然后在执行GetTopN(10)之前通过包装器加载它并且我得到了相同的错误。
我检查了所有相关变量(crs,z,x,n)似乎具有正确的值。我不知道这可能是什么原因,我真的很感激一些帮助!
以下是正在执行的内容(按顺序):
#INIT:
library(RODBC)
library(stats)
channel<- odbcConnect("data")
crs<-mat.or.vec(3000,5) #will hold correlations
n1<-seq(-33,0)
#Get whole series
z <- sqlQuery(channel,"SELECT RPos,M1,M2,M3,M4 FROM `data`.`z` ")
al_string <- "SELECT RPos,OpenTime FROM z JOIN actionlist on(OpenTime = pTime)"
trim_string<- "DELETE FROM ActionList WHERE OpenTime NOT IN (SELECT OpenTime FROM ReducedList)"
GetTopN<-function(n)
{
for(i in 1:nrow(actionlist))
{
crs[i,1]<-actionlist$OpenTime[i]
for(j in 2:ncol(z))
{
crs[i,j]<-cor(z[actionlist$RPos[i]+n1,j],x[,j])
}
}
avc <- (cbind(crs[,1],rowSums(crs[,2:5])))
sorted <- crs[order(avc[,2], decreasing=T),1]
topx<- head(sorted,n)
bottomx <- tail(sorted,n)
DF<-as.data.frame(c(topx,bottomx),row.names=NULL)
colnames(DF)[1]<-'OpenTime'
sqlSave(channel,dat=DF,tablename='ReducedList',append=F,rownames=F,safer=F)
sqlQuery(channel,trim_string)
}
curpTime <- 1275266400
actionlist<- sqlQuery(channel,al_string)
x<- sqlQuery(channel,paste('SELECT pTime,M1,M2,M3,M4 FROM z WHERE pTime <= ',curpTime,' AND
pTime > ',curpTime,'-(300*34) ORDER BY pTime ASC'))
GetTopN(10)
如果有帮助的话,我也保存了我的工作区(4.7mb):workspace
如果连接到我的MYSQL数据库会有所帮助,它应该在74.73.17.163:3306
答案 0 :(得分:14)
问题:actionlist$RPos[1000]
的值为21. n1
范围从-31到0.当您添加它们时,您会得到一个混合了正值和负值的向量,这不是允许子集化。
我是如何到达那里的:首先检查traceback()
:
traceback()
5: `[.data.frame`(z, actionlist$RPos[i] + n1, j) at #8
4: z[actionlist$RPos[i] + n1, j] at #8
3: is.data.frame(x) at #8
2: cor(z[actionlist$RPos[i] + n1, j], x[, j]) at #8
1: GetTopN(10)
这告诉我问题最有可能在actionlist$RPos[i] + n1
。然后我只添加了一个简单的print(i)
语句来告诉我哪个迭代是问题。 (或者,您可能刚刚手动检查了actionlist$RPos + n1
的问题点。
答案 1 :(得分:0)
对于tidyverse
用户,如果尝试对已分组的df进行子集设置,dplyr
会出现此错误,您需要先调用ungroup
。