File: trout2.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 (10 lines) | stat: -rw-r--r-- 249 bytes parent folder | download | duplicates (4) 
1
2
3
4
5
6
7
8
9
10
 
with(trout, {
conc <- N <- dead <-rep(0,6)
for(i in 1:6){
  conc[i] <- Concentr[i*4]
  for(j in 1:4){
      N[i] <- N[i] + Trouts[(i-1)*4+j]
      dead[i] <- dead[i] + Dead[(i-1)*4+j]
  }}
sm.binomial.bootstrap(log(conc), dead, N, 0.5, nboot=50)
})