Расчет процентного содержания в двумерном семействе

Question

Я хотел бы получить процентное соотношение 96u да от трех разных условий.

Как вы можете видеть, что я работаю с этим подгонкой (потому что я работаю с биномиальными данными)

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

Но каким код следует использовать, если я хочу, чтобы стать процентом из 96 U да для каждого условия seperatly и/или в сочетании с типом яйца

Answer 1

я не уверен, если я правильно понимаю, что вы имеете в виду думаю, что это то, что вы ищете

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

где df обозначает ваш входной кадр данных. Если вы хотите частот, сгруппированных по типу яиц или как условие и тип яйца, вам необходимо изменить переменные группировки соответственно в group_by()

Answer 2

Я думаю, что вы ищете что-то вроде

Construct " предсказания кадра», то есть соответствующие комбинации предикторов:

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

Предскажите:

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

в качестве альтернативы, так и для прогнозы, основанные только на состоянии, проверьте lsmeans package (не забывайте type="response": example)