Olá, tenho uma data.frame (df) 8 linhas x 8 colunas. Estou chamando as colunas de "ST[i]".
df <- structure(list(ST1 = c(58.69, 58.5, 58.5, 58.69, 58.69,
58.5, 58.69, 58.69), ST2 = c(68.7, 68.42, 68.42, 68.7,
68.7, 68.42, 68.7, 68.7), ST3 = c(69.15, 69.15, 68.83,
69.15, 69.15, 69.15, 69.15, 69.15), ST4 = c(78.99,
80.29, 78.99, 77.7, 78.99, 80.29, 78.99, 77.7), ST5 = c(75.65,
75.65, 75.65, 72.57, 71.07, 68.14, 66.7, 66.7), ST6 = c(71.35,
71.35, 79.83, 86.38, 83.09, 87.49, 87.49, 78.76), ST7 = c(73.61,
72.48, 78.22, 102.06, 82.33, 79.97, 112.48, 91.39), ST8 = c(77.57,
77.57, 77.57, 77.57, 77.57, 77.57, 77.57, 79.17)), .Names = c("ST1",
"ST2", "ST3", "ST4", "ST5", "ST6", "ST7", "ST8"), row.names = 122:129, class = "data.frame")
Eu preciso executar um comando waveslim::mra()
para cada uma dessas colunas.
O resultado deste comando, para cada coluna, será uma list com 12 vetores. Portanto, o produto final será 8 lists, cada list com 12 vetores.
Posso fazer, por exemplo:
dwc_ST1 <- mra(ST1, wf = "la8", method = "dwt", J=3, boundary = "reflection"
dwc_ST2 <- mra(ST2, wf = "la8", method = "dwt", J=3, boundary = "reflection")
...
dwc_ST8 <- mra(ST8, wf = "la8", method = "dwt", J=3, boundary = "reflection")
Mas seria muito bom se simplificasse isso. Pensei no for()
, em algo como:
names_dwc <- as.character(rep(NA, ncol(df)))
for (k in 1:ncol(df)){
names[k] <- paste('dwc_ST', k, sep = "")
dwc_ST[k] <- mra(df[ ,k], wf = "la8", method = "dwt", J=3, boundary = "reflection")
}
Mas tem aparecido:
Error in ST[k] <- mra(df,[ , k], wf = "la8", : object 'dwc_ST' not found
Eu pensei que eu teria solucionado este problema criando essas 8 ST's, com a linha (que utilizei acima):
names[k] <- paste('dwc_ST', k, sep = "")
Parece que não estou sabendo como criar os arquivos para serem gerados os resultados. Alguém sabe me dar uma ajuda?
dput
), para o seu problema ser reproduzível. Por último, ao utilizar funções que não fazem parte doR base
, indique no início do seu código qual o pacote necessáriodput
na sua pergunta.