Próbuję przyzwyczaić się do omawiania problemów w R. Chciałbym wywołać tę funkcjęglm() wewnątrz funkcji, ale nie działa, najwyraźniej z powodów ograniczających, których nie udało mi się naprawić za pomocą funkcjiassign() lubeval().

Oto uproszczona wersja:

ao <- function (y, x, phi = seq (0,1,0.1), dataset, weights) {
    logLikvector <- rep(0,length(phi))  # vector of zeros to be replaced thereafter
    for (i in 1:length(phi)) {          # loop to use glm()   
        fit <- glm (y ~ x, data = dataset, family = binomial, weights = weights)         
        logLikvector[i] <- logLik(fit)      # get log likelihood
    }
    logLikvector
}

Teraz chcę użyć funkcji ao () w moim zbiorze danych

    ao (y = Prop, x = Age, dataset = mydata, weights = Total) 

To nie działa, ale działa następująco:

ao (y = mydata$Prop, x = mydata$Age, dataset = mydata, weights = mydata$Total)

Czy ktoś wie, co robić?

Każda pomoc byłaby bardzo mile widziana !!!

Tak, oto jak powtórzyć mój problem z zestawem danych, którego używam

library("MASS")
data(menarche)
mydata <- menarche
mydata$Prop <- mydata$Menarche / mydata$Total