主页 / coding / 问题 / 79060100
Accepted
user8229029
Asked: 2024-10-07 05:17:40 +0800 CST

使用 R 中的应用函数对矩阵逐个元素运行函数

  • 772

我想使用多个基准、最低日温度和最高日温度来计算生长度日数。如果可能的话，我想不使用 for 循环来执行此操作以最小化我的代码。

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

我想找出我拥有的十天假温度数据中每个基地的生长度日数。这 3 个基地中的每一个应该有 10 个输出值。
我尝试了 mapply，如下所示：

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

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

这会产生错误的输出，老实说，我不知道它在做什么。我希望输出等于使用不同的基数运行上述函数 3 次。每个输出将有 10 个对应于不同基数的 gdd 值。我如何使用某种应用函数来做到这一点？或者我还需要其他东西吗？

3 个回答

  1. Best Answer

    我认为你应该将tmin和封装tmax在 中list，例如，

    > 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. 只是为了好玩，即使这个问题已经得到了一个很好的答案，一个循环/ *应用免费的替代方案，利用回收：

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

    以上产生与@ThomasIsCoding 的答案相同的输出。

    • 3

  3. outer可能是这里最好的选择。

    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

相关问题