File: muscle.q

package info (click to toggle)
sm 2.2-6.0-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 1,540 kB
  • sloc: f90: 259; ansic: 21; makefile: 2
file content (7 lines) | stat: -rw-r--r-- 225 bytes parent folder | download | duplicates (4) 
1
2
3
4
5
6
7
 
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)
})