Início / coding / Perguntas / 79012855
Accepted
Salvador
Asked: 2024-09-23 08:38:30 +0800 CST

Como escrever uma função para estimar intervalos de confiança

  • 772

Estou tentando escrever uma função para estimar intervalos de confiança em torno da média. Quero estimar intervalos por wkou monthusando pdcomo meu var, no entanto, depois de tentar por um tempo, estou obtendo NAs para todos os meus wke além disso, a função ignora o monthe retorna wkapenas por. Alguém poderia apontar o que estou fazendo errado? O que devo group_by dentro da função para poder alternar entre wke monthao usar a cifunção? Dados e código abaixo:

x <-  structure(list(wk = c(1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 
2, 3, 3, 3, 3, 3, 3), month = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
1, 1, 1, 1, 1, 1, 1, 1, 1, 1), xd = c(0.0370098838217444, 0.0391215970961887, 
0.0409865246890313, NA, NA, NA, NA, 0.015188, 0.007513, 0.012716, 
0.013962, 0.014259, 0.023553, 0.032122, 0.03953, 0.028946, 0.030769, 
0.030815, 0.029187, 0.022604), td = c(0.02701157274146, 0.0284158620689655, 
0.0296560389182058, NA, NA, NA, NA, 0.0125, 0.007396, 0.010856, 
0.011685, 0.011882, 0.018063, 0.023761, 0.028687, 0.021649, 0.022861, 
0.022892, 0.021809, 0.017432), pd = c(317.308439683869, 0, 126.719553152898, 
NA, NA, NA, NA, 2671.6, 3540.6976744186, 1270.35740604274, 1067.69362430466, 
688.099646524154, 317.444499806234, 420.941879550524, 280.475476696762, 
250.681324772507, 159.048160622895, 258.125109208457, 450.868907331836, 
0), year = c(2007, 2007, 2007, 2007, 2007, 2007, 2007, 2007, 
2007, 2007, 2007, 2007, 2007, 2007, 2007, 2007, 2007, 2007, 2007, 
2007)), class = "data.frame", row.names = c(NA, -20L))

x

library(tidyverse)
ci <- function(dat, by, var, conf.int = 0.90,...) {
dat <- na.omit(dat)
 dat %>%
  group_by(wk) %>%   #Should I group by wk here?
  do(data.frame(rbind(smean.cl.normal(.$var, conf.int = conf.int))))
}

ci(dat = x, by = month, var = pd) # Returns NAs

1 respostas

  1. Best Answer

    Eu usaria summariseinstead of doe then unnest_wider. E você precisa colocar o argumento "by" entre chaves. Se você group_by(wk)then sua função irá agrupar pela wkvariável não importa o que você especifique como argumento "by". Você também não precisa omitir o NA porque smean.cl.normalele omite por padrão.

    library(Hmisc)
library(dplyr)

ci <- function(dat, by, var, conf.int = 0.90,...) {
  #dat <- na.omit(dat)
  dat %>%
    group_by({{by}}) %>%   
    dplyr::summarise(mean=list(smean.cl.normal({{var}}, conf.int=conf.int))) %>%
    unnest_wider(mean)
}

    ci(dat = x, by = month, var = pd)
  month  Mean Lower Upper
  <dbl> <dbl> <dbl> <dbl>
1     1  739.  300. 1177.

    ci(dat = x, by = wk, var = pd)
# A tibble: 3 x 4
     wk  Mean Lower Upper
  <dbl> <dbl> <dbl> <dbl>
1     1  148. -121.  417.
2     2 1425.  528. 2323.
3     3  233.  111.  355.
    • 1

relate perguntas