(我的第一个数据示例示例,我希望这足够了)
我有这个xts对象weights
> dput(head(weights,4))
structure(c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.333333333333333, -0.333333333333333,
-0.333333333333333, -0.333333333333333, 0, 0, 0, 0, 0, 0.333333333333333,
0.333333333333333, 0.333333333333333, 0.333333333333333, 0.333333333333333,
0.333333333333333, 0.333333333333333, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0.333333333333333, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, -0.333333333333333, -0.333333333333333, -0.333333333333333,
-0.333333333333333, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0.333333333333333, 0.333333333333333, 0.333333333333333, 0.333333333333333,
0, 0, 0, 0, -0.333333333333333, -0.333333333333333, -0.333333333333333,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.333333333333333), .indexCLASS = "Date", .indexTZ = "UTC", tclass = "Date", tzone = "UTC", class = c("xts",
"zoo"), index = structure(c(946512000, 949276800, 951782400,
954460800), tzone = "UTC", tclass = "Date"), .Dim = c(4L, 39L
), .Dimnames = list(NULL, c("ABB.N", "ACTELION.N", "ADECCO.N",
"ALUSUISSE.N", "BALOISE.N", "CIBA.SC.N", "CLARIANT.N", "CS.GROUP.N",
"ELEKTROWATT.I", "EMS.CHEMIE", "GEBERIT.N", "GIVAUDAN.N", "JULIUS.BAER.N",
"KUDELSKI.I", "LAFARGEHOLCIM.N", "LONZA.N", "NESTLE.N", "NOBEL.BIOCARE.N",
"NOVARTIS.I", "NOVARTIS.N", "RICHEMONT.N", "ROCHE.GS", "SAIRGROUP.N",
"SBV.N", "SERONO..B..I", "SGS.N", "SULZER.N", "SWATCH.I", "SWATCH.N",
"SWISS.LIFE.HOLDING.AG.N", "SWISS.RE.N", "SWISSCOM.N", "SYNGENTA.N",
"SYNTHES.N", "TRANSOCEAN.N", "UBS.GROUP.N", "UBS.I", "UNAXIS.N20",
"ZURICH.INSURANCE.N")))
和两个数字向量Pos
和Neg
> dput(head(Pos,4))
c(`1999-12-30` = 0.34191022613305, `2000-01-31` = 0.411653443805752,
`2000-02-29` = 0.393078629707615, `2000-03-31` = 0.376686025847744
)
> dput(head(Neg,4))
c(`1999-12-30` = -1.32059915686837, `2000-01-31` = -1.27821843244263,
`2000-02-29` = -1.3142472840903, `2000-03-31` = -1.32647296498383)
我现在想将weights
中的每个负/正值除以Neg
或Pos
中的相应行(如果Neg
和Pos
等于转换为矩阵/数据框)
我尝试了for-loops,apply-functions和带条件的简单除法,但是我却没有得到想要的结果,对于如此简单的任务来说,这并不会那么困难...
From the last problem I had从xts-objects进行子集设置似乎也存在问题...
这些代码对我不起作用:(主要是尺寸或下标错误...)
For循环:
weights_adj <- apply(weights,1,function(x) if (coredata(x)[coredata(x)>0]){
coredata(x)[coredata(x)]/Pos
} else if (coredata(x)[coredata(x)<0]){
coredata(x)[coredata(x)]/Neg
})
应用功能:
weights_adj <- apply(weights,1,function(x) (coredata(x)[coredata(x)>0])/Pos)
其中+适用:
which_negative_weight <- t(apply(weights,1,function (x) which(x<0)))
result_negative <- apply(weights,1,function(x) x[,which_negative_weight[1,]])
其中+用于:
for (i in 1:nrow(weights)){
(weights[i,which_negative_weight[i,]]/Neg[i,])
}
@ K.Hua的解决方案行之有效,他很快为我提供了帮助,但是@ G.Grothendieck的解决方案更简单,更快捷!
答案 0 :(得分:0)
您尝试将大小为N的向量的样本除以大小为N的整个向量。此外,我认为在列上循环更容易:
weights_adj = weights
neg <- apply(weights,2,function (x) x<0)
for (i in 1:ncol(weights)){
weights_adj[neg[,i],i] = weights[neg[,i],i]/Neg[neg[,i]]
weights_adj[!neg[,i],i] = weights[!neg[,i],i]]/Pos[!neg[,i]]
}
答案 1 :(得分:0)
尝试一下:
weights / ((weights > 0) * Pos + (weights < 0) * Neg + (weights == 0))