R语言 在某些参数为常量而某些参数不为常量时应用函数

63lcw9qa  于 2023-04-03  发布在  其他
关注(0)|答案(1)|浏览(113)

我需要计算CVfromCI。在这个函数中,参数lowerupperpen是不同的;参数designalpharobust是常量,怎么样才能让代码更短,目前我每次都需要从头到尾写一遍。

library(PowerTOST)

CVfromCI(pe = 0.95, lower = 0.86, upper = 1.029, n = 24, design = "2x2", alpha = 0.05, robust = FALSE)
CVfromCI(pe = 0.94, lower = 0.897, upper = 1.027, n = 24, design = "2x2", alpha = 0.05, robust = FALSE)
CVfromCI(pe = 0.99, lower = 0.88, upper = 1.025, n = 24, design = "2x2", alpha = 0.05, robust = FALSE)
7d7tgy0s

7d7tgy0s1#

我们可以使用mapply来应用带有多个参数的函数CVfromCI

library(PowerTOST)

mapply(CVfromCI, 
       pe = c(0.95, 0.94, 0.99),
       lower = c(0.86, 0.897, 0.88),
       upper = c(1.029, 1.027, 1.025),
       n = 24,
       design = "2x2",
       alpha = 0.05,
       robust = FALSE)
# [1] 0.1824596 0.1371548 0.1547650
# Warning messages:
# 1: sigma based on pe & lower CL more than 10% different than
# sigma based on pe & upper CL. Check input. 
# 2: sigma based on pe & lower CL more than 10% different than
# sigma based on pe & upper CL. Check input. 
# 3: sigma based on pe & lower CL more than 10% different than
# sigma based on pe & upper CL. Check input.

我们也可以使用purrr包中的pmap_dbl。请注意,在使用pmap_dbl时,我们首先以列表的形式提供多个参数,然后提供函数。

library(purrr)

pmap_dbl(list(pe = c(0.95, 0.94, 0.99),
              lower = c(0.86, 0.897, 0.88),
              upper = c(1.029, 1.027, 1.025),
              n = 24,
              design = "2x2",
              alpha = 0.05,
              robust = FALSE),
         CVfromCI)
# [1] 0.1824596 0.1371548 0.1547650
# Warning messages:
# 1: sigma based on pe & lower CL more than 10% different than
# sigma based on pe & upper CL. Check input. 
# 2: sigma based on pe & lower CL more than 10% different than
# sigma based on pe & upper CL. Check input. 
# 3: sigma based on pe & lower CL more than 10% different than
# sigma based on pe & upper CL. Check input.

相关问题