Início / coding / Perguntas / 79060100
Accepted
user8229029
Asked: 2024-10-07 05:17:40 +0800 CST

Usando funções apply em R para executar funções elemento por elemento em matrizes

  • 772

Quero calcular graus-dias crescentes com várias bases, usando temperatura diária mínima e temperatura diária máxima. Gostaria de fazer isso sem um loop for, se possível, para minimizar meu código.

bases <- c(40,45,50)
tmin <- runif(10,30,70)
tmax <- runif(10,55,95)

Quero encontrar o número de graus-dia de crescimento para cada base para cada um dos dez dias para os quais tenho dados falsos de temperatura. Deve haver 10 valores de saída para cada uma das 3 bases.
Tentei o mapply da seguinte forma:

gdd_func <- function(tmin,tmax,bases){
   (tmin + tmax)/2 - bases}

test <- mapply(gdd_func,tmin,tmax,bases)

Isso produz uma saída incorreta onde eu honestamente não sei o que está fazendo. Quero que a saída seja igual a executar a função acima 3 vezes diferentes com as diferentes bases. Cada saída teria 10 valores gdd correspondentes às diferentes bases. Como eu faria isso usando uma função apply de algum tipo? Ou preciso de algo mais?

3 respostas

  1. Best Answer

    Acho que você deve encapsular tmine tmaxdentro de um list, por exemplo,

    > mapply(gdd_func, list(tmin), list(tmax), bases)
           [,1]      [,2]      [,3]
 [1,] 19.119648 14.119648  9.119648
 [2,]  6.119249  1.119249 -3.880751
 [3,] 19.523350 14.523350  9.523350
 [4,] 13.023240  8.023240  3.023240
 [5,] 31.083015 26.083015 21.083015
 [6,] 32.281061 27.281061 22.281061
 [7,] 18.062584 13.062584  8.062584
 [8,] 18.047198 13.047198  8.047198
 [9,] 21.518514 16.518514 11.518514
[10,] 28.440563 23.440563 18.440563
    • 3

  2. Só por diversão e mesmo que essa pergunta já tenha uma boa resposta, uma alternativa de loop/*aplicação gratuita, explorando a reciclagem:

    gdd_func(tmin, tmax, rep(bases, each=length(tmin))) |> 
matrix(ncol = length(bases))

    O resultado acima produz a mesma saída que a resposta de @ThomasIsCoding.

    • 3

  3. outeré provavelmente a melhor opção aqui.

    outer((tmin + tmax)/2, bases, "-")
#>           [,1]      [,2]      [,3]
#>  [1,] 13.70856  8.708559  3.708559
#>  [2,] 17.66631 12.666310  7.666310
#>  [3,] 21.00194 16.001937 11.001937
#>  [4,] 14.82080  9.820800  4.820800
#>  [5,] 21.35614 16.356138 11.356138
#>  [6,] 33.96390 28.963897 23.963897
#>  [7,] 28.37265 23.372649 18.372649
#>  [8,] 30.85054 25.850542 20.850542
#>  [9,] 32.30424 27.304244 22.304244
#> [10,] 40.03540 35.035405 30.035405
    • 2

relate perguntas