peak.lme<-lme(mu ~ 0 + betaxolol + bimatoprost + brimonidine + brinzolamide + dorzolamide + latanoprost + placebo + travoprost, ~1|trtpair,data=peakdata, weights=varConstPower(form=~sigma, fixed=list(power=1)))
谢谢！