glm с asin (sqrt ...) ссылкой

Question

В моем наборе данных есть доля времени, которое насекомые проводят на цветах возле химической приманки. tfac это время суток. Будет запущен квазибиномиальный glm, но он чувствует себя не так - это не подсчет данных. Как насчет:

timeflower<-c(0.00000000,0.11998400,0.00000000,0.00000000,0.37103531,0.95535096,
0.99931082,0.00000000,0.99778271,0.26382646,0.00000000,0.00000000,
0.00000000,0.99942890,0.00000000,0.00000000,0.00000000,0.99988507,
0.99988414,0.99957465,0.87557604,0.99991936,0.99991313,0.00000000,
0.00000000,0.99961553,0.00000000,0.63080219,0.37604695,0.76897359,
0.90740904,0.82835184,0.99984851,0.99996048,0.89557430,0.03569983,
0.00000000,0.35877969,0.39578965,0.36295767,0.00000000,0.99993175,
0.49829998,0.49998175,0.35941774,0.59701493,0.08104158,0.00000000,
0.00000000,0.74677242,0.57419931,0.82720751,0.00000000,0.12395962,
0.62372233,0.58498388,0.00000000,0.00000000,0.95219155,0.89045866,
0.00000000,0.19296553,0.00000000,0.00000000)
tfac<-as.factor(c(9,11,13,15,9,11,13,15,9,11,13,15,9,11,13,15,9,11,13,15,9,,11,13,15,9,
11,13,15,9,11,13,15,9,11,13,15,9,11,13,15,9,11,13,15,9,11,13,15,9,11,
13,15,9,11,13,15,9,11,13,15,9,11,13,15))
#####angular link
vang <- function() {
    ## link
    linkfun <- function(y) asin(sqrt(y))
    ## inverse link
    linkinv <- function(eta)  sin(eta)^2
    ## derivative of invlink wrt eta
    mu.eta <- function(eta) { 2*sin(eta)*cos(eta) }
    valideta <- function(eta) TRUE
    link <- "asin(sqrt(y))"
    structure(list(linkfun = linkfun, linkinv = linkinv,
                   mu.eta = mu.eta, valideta = valideta, 
                   name = link),
              class = "link-glm")
}
vv<-vang()
m3<-glm(timeflower~tfac,family=gaussian(link=vv))
summary(m3)