In writing our my solution to the problem, I realized that I need to change one of the hints. Use this
p.poly <- step(s.1)*step(s.2)*(s.2/(s.1+s.2)) + step(-s.1)*step(s.2)
Instead of what I suggested in the notes.
Posted by Kent at February 24, 2004 9:37 AM | TrackBack