taxa de falsos positivos na diferença entre as médias

votos
0

Eu tenho tentado usar o seguinte código para avaliar a taxa de falsos positivos de uma diferença de meios estimativa.

sm_n <- 10
dat <- data.frame(b=rep(c(a,b),c(4,6)),
            x = c(0,1,0,0,0,1,0,0,1,1),
            nb=rep(c(4,6),c(4,6)))
dat$y0 <- 1000*dat$x+c(0,0,0,0,1,2,3,400,5000,60000)
## summary(lm(y0~x,data=dat))$r.squared
## dat$Z <-  ifelse(dat$b==a,complete_ra(N=4), complete_ra(N=6,m=2))
set.seed(12345)
dat$Z  <- block_ra(blocks = dat$b, block_m = c(2,2))
dat$tau <- c(2500,0,2500,0,25000,25000,0,0,0,50000)
dat$y1 <- dat$y0 + dat$tau
trueATE <- with(dat,mean(y1-y0))
trueRankATE <- with(dat,mean(rank(y1)-rank(y0)))
dat$Y <- with(dat,Z*y1 + (1-Z)*y0)
dat$ZF <- factor(dat$Z)
dat %<>% group_by(b) %>% mutate(pi=mean(Z), #prob treated
                  nbwt=Z/pi + (1-Z)/(1-pi))

A estimativa em si funciona bem, e assim que faz a função que gera uma nova atribuição aleatória com nenhuma relação sistemática entre Z e qualquer outra coisa. No entanto, uma vez que eu tente executar a função para extrair os valores de p, eu recebo o seguinte erro:

'Erro no jogo (x, mesa, NoMatch = 0D): 'match' requer argumentos vetor'

O problema parece estar a ocorrer em 'difference_in_means (Y ~ newexp (b), os blocos = b, dados = DAT)'

est3 <- difference_in_means(Y~Z,blocks=b,data=dat)

summary(est3) 

newexp <- function(b){
    ## A new random assignment with no systematic relationship between Z and anything else
    Z <- block_ra(blocks = dat$b, block_m = c(2,2))
    return(Z)
}

get_p_val_from_difference_in_means<-function(){
    ## First, shuffle to break the relationship, to make the truth zero
    ## Then, test
    theest <- difference_in_means(Y~newexp(b), blocks=b,data=dat)
    thep <-  coef(summary(theest))[1,4]
    thep
}

est3ps <- replicate(1000, get_p_val_from_difference_in_means())

mean(est3ps < .05)
plot(ecdf(est3ps ), ylim = c(0,1), xlim = c(0,1))
abline(0,1)
Publicado 10/10/2019 em 00:50
fonte usuário
Em outras línguas...                            

Cookies help us deliver our services. By using our services, you agree to our use of cookies. Learn more