Berechnung Prozentsatz in einem Glemer mit einer Binomialfamilie

Question

Ich möchte den Prozentsatz von 96u ja meiner drei verschiedenen Bedingungen erhalten.

Wie Sie sehen, ich bin mit diesem fit arbeitet (weil ich mit binomischen Daten arbeiten bin)

fit96u=glmer(cbind(96u.yes,96u.No)~Condition+Type.eggs+(1|number),family=binomial(link=logit),data=data) 

Aber die Code sollte ich verwenden, wenn ich der Prozentsatz werden wollen von 96 u ja für seperatly jede Bedingung oder/und mit einer Art von Ei kombiniert

Answer 1

ich bin nicht sicher, ob ich das richtig verstehen, was Sie dies bedeuten, denken ist, was Sie für

suchen

library(dplyr) 
percByCondition <- df %>% group_by(Condition) %>% summarize(percYes = mean(96u.yes)) 

Dabei steht df für den Eingabedatenrahmen. Wenn Sie die Frequenzen nach Art der Ei oder beide Zustand und Art von Ei gruppiert werden sollen, müssen Sie die Gruppenvariablen entsprechend in group_by ändern()

Answer 2

Ich glaube, du bist für so etwas wie

Baukonstruktion suchen " Prädiktionseinzelbildspeicher“, dh relevante Kombinationen von Prädiktoren:

pframe <- with(my_data, 
    expand.grid(Condition=levels(Condition), 
       Type.Eggs=levels(Type.Eggs))) 

Predict:

predict(fit96u,newdata=pframe, 
      type="response", ## predict on probability (not logit) scale 
      re.form=~0)  ## ignore grouping when predicting 

Alternativ und für die Vorhersagen basieren nur auf Condition, überprüfen Sie die lsmeans package (nicht vergessen type="response": example)