[R-es] ceiling con max duplica tamaño matriz

[Marcelino de la Cruz Rot]

Hola Sebastián:

No nos has enseñado cómo es idx_epsilon.

Haz

head(idx_epsilon)

y comprueba que sus valores vayan de 1 a 469. Si no incluyen estos 
valores es normal que después del bucle las posiciones 1:469 de v1 se 
queden como estaban.

Saludos,
Marcelino


El 22/08/2022 a las 22:06, Sebastian Kruk escribió:
> Estimados,
>
> Tengo:
> v1
> v_DEM_i_d_m
>
> Donde
>
>> IN_ES_NUM = 469
>> v1 = 1:IN_ES_NUM
>> head(v_DEM_i_d_m)
>              Mes 1        Mes 2        Mes 3        Mes 4        Mes 5
>        Mes 6        Mes 7        Mes 8        Mes 9       Mes 10
> Mes 11       Mes 12
> ES 1  7.811032258  8.982678571  9.563935484  8.411500000  7.788354839
> 8.976933333  9.517419355  8.902258065  9.328600000  9.391096774
> 9.866800000 10.317838710
> ES 2  0.641419355  1.174714286  1.353838710  1.299900000  1.387096774
> 0.905033333  1.140935484  0.880387097  1.237766667  1.292677419
> 1.327166667  1.387838710
> ES 3  0.003225806  0.003571429  0.003225806  0.003333333  0.003225806
> 0.003333333  0.003225806  0.003225806  0.003333333  0.003225806
> 0.003333333  0.003483871
> ES 4 13.475258065 13.240285714 15.168838710 11.492233333 12.792774194
> 12.989566667 14.996064516 12.911774194 13.835400000 14.635612903
> 14.541266667 15.206000000
> ES 5 18.711258065 23.263535714 24.875967742 23.478000000 27.372548387
> 27.509800000 28.649096774 26.986935484 33.831833333 20.810419355
> 24.669400000 25.797096774
> ES 6 16.112419355 17.903678571 18.250548387 16.712066667 17.353741935
> 17.137233333 18.501000000 17.613258065 19.115766667 20.773354839
> 20.990566667 21.950096774
>
>> class((v_DEM_i_d_m))
> [1] "matrix" "array"
>
>> str(v_DEM_i_d_m)
>   num [1:469, 1:12] 7.81103 0.64142 0.00323 13.47526 18.71126 ...
>   - attr(*, "dimnames")=List of 2
>    ..$ : chr [1:469] "ES 1" "ES 2" "ES 3" "ES 4" ...
>    ..$ : chr [1:12] "Mes 1" "Mes 2" "Mes 3" "Mes 4" ...
>
> Si hago:
>
> for (i in idx_epsilon) {
> v1[i] = ceiling(max(v_DEM_i_d_m[i,]))
> }
>
> v1 queda así:
>
>> head(v1)
> 1 2 3 4 5 6
>
>> class(v1)
> [1] "numeric"
>
>> str(v1)
>   Named num [1:938] 1 2 3 4 5 6 7 8 9 10 ...
>   - attr(*, "names")= chr [1:938] "" "" "" "" ...
>
> No entiendo porque los primeros 469 son los números de 1 a 469.
>
> ¿Que estoy haciendo mal?
>
> Saludos,
>
> Sebastián.
>
> _______________________________________________
> R-help-es mailing list
> R-help-es using r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es


-- 
Marcelino de la Cruz Rot
Depto. de Biología y Geología
Física y Química Inorgánica
Universidad Rey Juan Carlos
Móstoles España