<div dir="ltr">Olá,<div><br></div><div>Estou trabalhando neste código aqui:</div><div><br></div><div><div>Código:</div><div><br></div><div><div># Para selecionar o vm_data>0</div><div>rownames.positive <- function(data, i) {</div><div> return <- colnames(data)[!<a href="http://is.na">is.na</a>(data[i,])]</div><div>}</div><div><br></div><div># Para selecionar os rets>0</div><div>rownames.positive1 <- function(data, j) {</div><div> return <- colnames(data)[!<a href="http://is.na">is.na</a>(data[j,])]</div><div>}</div><div><br></div></div><div><br></div><div><div>valid.tickers <- function(vm_data,rets, years, years1) { #years = rownames(vm_data), years1 = rownames(rets)</div><div> nomeativos <- list()</div><div> for (i in 1:dim(vm_data)[[1]]){ </div><div> for (j in 1: dim(rets)[[1]]){ # i=2; j=42</div><div> a <- rownames.positive(vm_data, years[i])</div><div> b <- rownames.positive1(rets[,a], years1[j])</div><div> nomeativos[[j]] <- b</div><div> }}</div><div> names(nomeativos)<-years1[1:j]</div><div> return(nomeativos)</div><div>}</div></div><div><br></div><div><div><br></div><div>Estou tentando fazer dois filtros em um banco de dados, no qual primeiro verifico se existe vm_data>0 através da função rownames.positive e, consequentemente, para os ativos que passaram neste filtro, aplico a função rownames.positive1 para identificar rets>0. Com isso, uso a função valid.tickers para aplicar os dois filtros e retornar uma lista com apenas os ativos (AT) válidos.</div><div><br></div><div>Minhas variáveis vm_data e rets estao abaixo:<br></div><div><br></div><div>vm_data</div><div><br></div><div><div> At1 At2 At3 At4 At5 At6</div><div>2000-01-01 NA 7.478553 106.05138 154.9142 246.72923 84.74537</div><div>2000-02-01 56.45321 22.248116 49.97758 124.3428 100.43189 57.89444</div><div>2000-03-01 10.53851 222.545990 15.49008 36.3656 47.19245 51.59810</div></div><div><br></div><div>head(rets,10)</div><div><div> At1 At2 At3 At4 At5 At6</div><div>2000-01-01 -0.84032326 1.18251249 0.2791886 -0.3506028 -1.2780974 0.61644019</div><div>2000-01-03 1.93239850 -0.05272626 0.2532625 1.3811129 1.1809074 -2.00703218</div><div>2000-01-04 0.56596405 0.36471560 1.2898374 1.0211276 -1.0492678 -0.13936848</div><div>2000-01-05 -0.72244969 -2.13653307 0.6436690 -1.4931114 1.4653181 -0.65127247</div><div>2000-01-06 -0.15135535 -1.15261317 0.1577168 0.1466359 1.4614018 1.29129399</div><div>2000-01-07 0.08811143 -0.93320399 -0.6568275 -1.1365693 -0.1147392 1.16723028</div><div>2000-01-08 -1.01018908 -0.79147444 -0.2458063 0.6530664 -0.2351105 -1.56632861</div><div>2000-01-10 -1.41642902 -0.73998066 NA -0.1959526 1.6433762 -1.83665410</div><div>2000-01-11 -1.88839674 -1.07155233 -0.9909747 1.5976778 1.1155120 0.02490766</div><div>2000-01-12 -0.68444922 -0.69488280 -0.3216356 0.0174740 -0.2322832 0.48464970</div></div><div><br></div><div>tail(rets)</div><div><br></div><div><div> At1 At2 At3 At4 At5 At6</div><div>2000-02-29 -0.061260477 0.3558457 1.0351393 0.5690812 -0.5696858 -0.4889139<br></div><div>2000-03-01 1.023446137 0.2470349 0.2481944 0.7226727 0.9731186 -0.6501507</div></div><div><br></div><div>O problema que não estou conseguindo resolver é que meu "for" para o contador j está errado. Quero que para cada i (mês) de vm_data selecionar apenas os rets>0 para os dias desse respectivo mês usando contador j. Da forma como está no meu código valid.tickers está retornando através de nomeativos que na data $'2000 - 01 - 01 os AT de 1 a 6 passaram nos dois filtros, o que não é verdade, pois para todo o mês 01, o AT1 já não deveria aparecer, pois a variável vm_data na linha 1 tem AT1 = NA. Assim, valid.tickers parece armazenar apenas o i=3 para j= 1:dim(rets)[[1]], pois se eu fizer manualmente i=1 e j=1, o resultado está correto, mas quando uso o "for" não.</div><div><br></div><div>Portanto, como posso fazer i=1,2 e 3 sendo que para cada i quero delimitar que o contador j possa varrer apenas aquelas observações da variável rets nas quais o mês e ano de rets são iguais aos de i? Quero um resultado final com uma lista de length=52, que é a dim(rets)[[1]], no qual, por exemplo, nomeativos[[j]] = b retorne que $'2000-01-01' os AT de 2 a 6 sao válidos (pois, vm_data para o mês 01 no AT1 = NA e no rets todos os demais AT são maiores que zero), de igual maneira, $'2000-02-18' retorne apenas os AT 1, 2, 4, 5 e 6 (pois, não há NA no vm_data para o mês 2, mas há um NA na variável rets para o AT3 nessa data).</div><div><br></div><div>Obrigada pela ajuda!</div><div><br></div><div><br></div><div><div><br></div><br><div class="gmail_signature"><div><b><span style="background-color:rgb(255,255,255);color:rgb(0,0,102)"><br style="background-color:rgb(255,255,255);color:rgb(0,0,102)"></span></b><br><br style="color:rgb(0,0,102);background-color:rgb(0,0,102)"><b><span style="color:rgb(0,0,102);background-color:rgb(255,255,255)"></span><br style="color:rgb(0,0,102)"></b><br></div></div>
</div></div></div></div>