File: trw_lf.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 (13 lines) | stat: -rw-r--r-- 350 bytes parent folder | download | duplicates (4) 
1
2
3
4
5
6
7
8
9
10
11
12
13
 
with(trawl, {
ind       <- (Year == 0 & Zone == 1)
score     <- Score1[ind]
longitude <- Longitude[ind]
par(mfrow=c(1,2))
plot(longitude, score, pch="o")
ls.fit <- lm(score ~ longitude)
abline(ls.fit$coefficients, lty=3)
plot(longitude, residuals(ls.fit), pch="o")
abline(0,0,lty=3)
par(mfrow=c(1,1))
print(summary(lm(score ~ poly(longitude,2))))
})