Piszę funkcję, w której chcę przekazać niektóre argumenty do funkcji crrstep (pakiet „crrstep”), ale napotkałam problem: jakoś argument „zdarzenie” w mojej funkcji nie jest rozpoznawany, gdy wprowadzam go w crrstep. Domyślam się, że crrstep wygląda w innym środowisku niż to, na które chcę, ale nawet po wielu godzinach poszukiwań rozwiązań w sieci nie wydaje mi się, jak to rozwiązać (jestem dość niedoświadczony w programowaniu ..) . Każda pomoc byłaby bardzo mile widziana!

Oto kilka danych symulacyjnych (skorygowany przykład z dokumentacji crrstep) i przykład mojego kodu:

n <- 500
ftime <- rexp(n)
fstatus <- sample(0:2,n,replace=TRUE)
testdata <- matrix(runif(8*n),nrow=n)
testdata <- cbind(ftime,fstatus,testdata)
dimnames(testdata)[[2]] <- c('ftime','fstatus','x1','x2','x3','x4','x5','x6','x7','x8')
testdata <- as.data.frame(testdata)
formula1 <- ftime ~ 1 + x1 + x2 + x3 + x4 + x5 + x6 + x7 + x8
rm(fstatus,ftime,n)

test.fun <- function(x,data,event){
require(crrstep)
select.mod<- crrstep(formula=x,,etype=event, failcode=1, cencode=0,data=data, 
                   direction = "backward", criterion = "AIC", crr.object = TRUE, 
                   trace = FALSE)
#Rest of function omitted for now
print(select.mod)
}

#Test
test.fun(x=formula1,data=testdata,event=fstatus) 
#I get: Error in eval(expr, envir, enclos) : object 'event' not found"

Wielkie dzięki! Obrabować