File: citrate.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 (12 lines) | stat: -rw-r--r-- 373 bytes parent folder | download | duplicates (4) 
1
2
3
4
5
6
7
8
9
10
11
12
 
with(citrate, {
Citrate<-as.matrix(citrate)
nSubj<-dim(Citrate)[1]
nTime<-dim(Citrate)[2]
Time<-(1:nTime)
plot(c(min(Time),max(Time)), c(min(Citrate),max(Citrate)),
    type="n", xlab="time", ylab="Citrate")
for(i in 1:nSubj) lines(Time,as.vector(Citrate[i,]))
pause()
a <- sm.rm(y=Citrate, display.rice=TRUE)
sm.regression(Time,a$aux$mean,h=1.2,hmult=1,add=TRUE,lty=3)
})