We have already read this website. And referred the code to write our code. We guess the attachment code we writed main problem appears in the following part. print Param[c={"R0" "B" "lambda" "simga"} label=""] LogLik; p = param[1,]; print p; f = LogLik(p); print f; print Param[c={"R0" "B" "lambda" "simga"} label=""] LogLik; In addition, we don't use date to test . We just use sas to run and revise our code. We want to make sure our code is correct first. If possible, wish you could give us more direction to revise our code again.
... View more