inst/scripts/muscle.q

with(muscle, {
TypeI <- TypeI.P+TypeI.R+TypeI.B
sm.poisson(log(TypeI), TypeII, 0.25, display="se")
pm <- glm(TypeII ~ log(TypeI), family=poisson)
lines(sort(log(TypeI)), fitted(pm)[order(log(TypeI))], 
      col=4, lty=6)
})

Try the sm package in your browser

Any scripts or data that you put into this service are public.

sm documentation built on May 29, 2024, 2:28 a.m.