data nlmm; input y base trt age Visit subject; cards; 5 1.01160091167848 0 3.43398720448515 -0.3 1 3 1.01160091167848 0 3.43398720448515 -0.1 1 3 1.01160091167848 0 3.43398720448515 0.1 1 3 1.01160091167848 0 3.43398720448515 0.3 1 3 1.01160091167848 0 3.40119738166216 -0.3 2 5 1.01160091167848 0 3.40119738166216 -0.1 2 3 1.01160091167848 0 3.40119738166216 0.1 2 3 1.01160091167848 0 3.40119738166216 0.3 2 2 0.405465108108164 0 3.2188758248682 -0.3 3 4 0.405465108108164 0 3.2188758248682 -0.1 3 0 0.405465108108164 0 3.2188758248682 0.1 3 5 0.405465108108164 0 3.2188758248682 0.3 3 4 0.693147180559945 0 3.58351893845611 -0.3 4 4 0.693147180559945 0 3.58351893845611 -0.1 4 1 0.693147180559945 0 3.58351893845611 0.1 4 4 0.693147180559945 0 3.58351893845611 0.3 4 7 2.80336038090653 0 3.09104245335832 -0.3 5 18 2.80336038090653 0 3.09104245335832 -0.1 5 9 2.80336038090653 0 3.09104245335832 0.1 5 21 2.80336038090653 0 3.09104245335832 0.3 5 5 1.90954250488444 0 3.36729582998647 -0.3 6 2 1.90954250488444 0 3.36729582998647 -0.1 6 8 1.90954250488444 0 3.36729582998647 0.1 6 7 1.90954250488444 0 3.36729582998647 0.3 6 6 1.09861228866811 0 3.43398720448515 -0.3 7 4 1.09861228866811 0 3.43398720448515 -0.1 7 0 1.09861228866811 0 3.43398720448515 0.1 7 2 1.09861228866811 0 3.43398720448515 0.3 7 40 2.56494935746154 0 3.73766961828337 -0.3 8 20 2.56494935746154 0 3.73766961828337 -0.1 8 21 2.56494935746154 0 3.73766961828337 0.1 8 12 2.56494935746154 0 3.73766961828337 0.3 8 5 1.74919985480926 0 3.61091791264422 -0.3 9 6 1.74919985480926 0 3.61091791264422 -0.1 9 6 1.74919985480926 0 3.61091791264422 0.1 9 5 1.74919985480926 0 3.61091791264422 0.3 9 14 0.916290731874155 0 3.3322045101752 -0.3 10 13 0.916290731874155 0 3.3322045101752 -0.1 10 6 0.916290731874155 0 3.3322045101752 0.1 10 0 0.916290731874155 0 3.3322045101752 0.3 10 26 2.56494935746154 0 3.58351893845611 -0.3 11 12 2.56494935746154 0 3.58351893845611 -0.1 11 6 2.56494935746154 0 3.58351893845611 0.1 11 22 2.56494935746154 0 3.58351893845611 0.3 11 12 2.11021320034659 0 3.17805383034795 -0.3 12 6 2.11021320034659 0 3.17805383034795 -0.1 12 8 2.11021320034659 0 3.17805383034795 0.1 12 4 2.11021320034659 0 3.17805383034795 0.3 12 4 1.50407739677627 0 3.13549421592915 -0.3 13 4 1.50407739677627 0 3.13549421592915 -0.1 13 6 1.50407739677627 0 3.13549421592915 0.1 13 2 1.50407739677627 0 3.13549421592915 0.3 13 7 2.35137525716348 0 3.58351893845611 -0.3 14 9 2.35137525716348 0 3.58351893845611 -0.1 14 12 2.35137525716348 0 3.58351893845611 0.1 14 14 2.35137525716348 0 3.58351893845611 0.3 14 16 3.07961375753469 0 3.25809653802148 -0.3 15 24 3.07961375753469 0 3.25809653802148 -0.1 15 10 3.07961375753469 0 3.25809653802148 0.1 15 9 3.07961375753469 0 3.25809653802148 0.3 15 11 2.52572864430826 0 3.25809653802148 -0.3 16 0 2.52572864430826 0 3.25809653802148 -0.1 16 0 2.52572864430826 0 3.25809653802148 0.1 16 5 2.52572864430826 0 3.25809653802148 0.3 16 0 1.50407739677627 0 3.3322045101752 -0.3 17 0 1.50407739677627 0 3.3322045101752 -0.1 17 3 1.50407739677627 0 3.3322045101752 0.1 17 3 1.50407739677627 0 3.3322045101752 0.3 17 37 3.32323584019244 0 3.43398720448515 -0.3 18 29 3.32323584019244 0 3.43398720448515 -0.1 18 28 3.32323584019244 0 3.43398720448515 0.1 18 29 3.32323584019244 0 3.43398720448515 0.3 18 3 1.50407739677627 0 3.46573590279973 -0.3 19 5 1.50407739677627 0 3.46573590279973 -0.1 19 2 1.50407739677627 0 3.46573590279973 0.1 19 5 1.50407739677627 0 3.46573590279973 0.3 19 3 1.6094379124341 0 3.04452243772342 -0.3 20 0 1.6094379124341 0 3.04452243772342 -0.1 20 6 1.6094379124341 0 3.04452243772342 0.1 20 7 1.6094379124341 0 3.04452243772342 0.3 20 3 1.09861228866811 0 3.36729582998647 -0.3 21 4 1.09861228866811 0 3.36729582998647 -0.1 21 3 1.09861228866811 0 3.36729582998647 0.1 21 4 1.09861228866811 0 3.36729582998647 0.3 21 3 0.810930216216329 0 3.04452243772342 -0.3 22 4 0.810930216216329 0 3.04452243772342 -0.1 22 3 0.810930216216329 0 3.04452243772342 0.1 22 4 0.810930216216329 0 3.04452243772342 0.3 22 2 1.44691898293633 0 3.46573590279973 -0.3 23 3 1.44691898293633 0 3.46573590279973 -0.1 23 3 1.44691898293633 0 3.46573590279973 0.1 23 5 1.44691898293633 0 3.46573590279973 0.3 23 8 1.94591014905531 0 3.2188758248682 -0.3 24 12 1.94591014905531 0 3.2188758248682 -0.1 24 2 1.94591014905531 0 3.2188758248682 0.1 24 8 1.94591014905531 0 3.2188758248682 0.3 24 18 2.62103882411258 0 3.40119738166216 -0.3 25 24 2.62103882411258 0 3.40119738166216 -0.1 25 76 2.62103882411258 0 3.40119738166216 0.1 25 25 2.62103882411258 0 3.40119738166216 0.3 25 2 0.810930216216329 0 3.68887945411394 -0.3 26 1 0.810930216216329 0 3.68887945411394 -0.1 26 2 0.810930216216329 0 3.68887945411394 0.1 26 1 0.810930216216329 0 3.68887945411394 0.3 26 3 0.916290731874155 0 2.94443897916644 -0.3 27 1 0.916290731874155 0 2.94443897916644 -0.1 27 4 0.916290731874155 0 2.94443897916644 0.1 27 2 0.916290731874155 0 2.94443897916644 0.3 27 13 2.46385324059017 0 3.09104245335832 -0.3 28 15 2.46385324059017 0 3.09104245335832 -0.1 28 13 2.46385324059017 0 3.09104245335832 0.1 28 12 2.46385324059017 0 3.09104245335832 0.3 28 11 2.94443897916644 1 2.89037175789616 -0.3 29 14 2.94443897916644 1 2.89037175789616 -0.1 29 9 2.94443897916644 1 2.89037175789616 0.1 29 8 2.94443897916644 1 2.89037175789616 0.3 29 8 2.2512917986065 1 3.46573590279973 -0.3 30 7 2.2512917986065 1 3.46573590279973 -0.1 30 9 2.2512917986065 1 3.46573590279973 0.1 30 4 2.2512917986065 1 3.46573590279973 0.3 30 0 1.55814461804655 1 2.99573227355399 -0.3 31 4 1.55814461804655 1 2.99573227355399 -0.1 31 3 1.55814461804655 1 2.99573227355399 0.1 31 0 1.55814461804655 1 2.99573227355399 0.3 31 3 0.916290731874155 1 3.40119738166216 -0.3 32 6 0.916290731874155 1 3.40119738166216 -0.1 32 1 0.916290731874155 1 3.40119738166216 0.1 32 3 0.916290731874155 1 3.40119738166216 0.3 32 2 1.55814461804655 1 2.89037175789616 -0.3 33 6 1.55814461804655 1 2.89037175789616 -0.1 33 7 1.55814461804655 1 2.89037175789616 0.1 33 4 1.55814461804655 1 2.89037175789616 0.3 33 4 1.79175946922806 1 3.17805383034795 -0.3 34 3 1.79175946922806 1 3.17805383034795 -0.1 34 1 1.79175946922806 1 3.17805383034795 0.1 34 3 1.79175946922806 1 3.17805383034795 0.3 34 22 2.04769284336526 1 3.40119738166216 -0.3 35 17 2.04769284336526 1 3.40119738166216 -0.1 35 19 2.04769284336526 1 3.40119738166216 0.1 35 16 2.04769284336526 1 3.40119738166216 0.3 35 5 1.25276296849537 1 3.55534806148941 -0.3 36 4 1.25276296849537 1 3.55534806148941 -0.1 36 7 1.25276296849537 1 3.55534806148941 0.1 36 4 1.25276296849537 1 3.55534806148941 0.3 36 2 1.01160091167848 1 3.29583686600433 -0.3 37 4 1.01160091167848 1 3.29583686600433 -0.1 37 0 1.01160091167848 1 3.29583686600433 0.1 37 4 1.01160091167848 1 3.29583686600433 0.3 37 3 2.81839825827108 1 2.99573227355399 -0.3 38 7 2.81839825827108 1 2.99573227355399 -0.1 38 7 2.81839825827108 1 2.99573227355399 0.1 38 7 2.81839825827108 1 2.99573227355399 0.3 38 4 2.32727770558442 1 3.09104245335832 -0.3 39 18 2.32727770558442 1 3.09104245335832 -0.1 39 2 2.32727770558442 1 3.09104245335832 0.1 39 5 2.32727770558442 1 3.09104245335832 0.3 39 2 0.559615787935423 1 3.3322045101752 -0.3 40 1 0.559615787935423 1 3.3322045101752 -0.1 40 1 0.559615787935423 1 3.3322045101752 0.1 40 0 0.559615787935423 1 3.3322045101752 0.3 40 0 1.70474809223843 1 3.13549421592915 -0.3 41 2 1.70474809223843 1 3.13549421592915 -0.1 41 4 1.70474809223843 1 3.13549421592915 0.1 41 0 1.70474809223843 1 3.13549421592915 0.3 41 5 1.17865499634165 1 3.68887945411394 -0.3 42 4 1.17865499634165 1 3.68887945411394 -0.1 42 0 1.17865499634165 1 3.68887945411394 0.1 42 3 1.17865499634165 1 3.68887945411394 0.3 42 11 2.4423470353692 1 3.49650756146648 -0.3 43 14 2.4423470353692 1 3.49650756146648 -0.1 43 25 2.4423470353692 1 3.49650756146648 0.1 43 15 2.4423470353692 1 3.49650756146648 0.3 43 10 2.19722457733622 1 3.04452243772342 -0.3 44 5 2.19722457733622 1 3.04452243772342 -0.1 44 3 2.19722457733622 1 3.04452243772342 0.1 44 8 2.19722457733622 1 3.04452243772342 0.3 44 19 2.2512917986065 1 3.55534806148941 -0.3 45 7 2.2512917986065 1 3.55534806148941 -0.1 45 6 2.2512917986065 1 3.55534806148941 0.1 45 7 2.2512917986065 1 3.55534806148941 0.3 45 1 0.559615787935423 1 3.2188758248682 -0.3 46 1 0.559615787935423 1 3.2188758248682 -0.1 46 2 0.559615787935423 1 3.2188758248682 0.1 46 3 0.559615787935423 1 3.2188758248682 0.3 46 6 2.19722457733622 1 3.25809653802148 -0.3 47 10 2.19722457733622 1 3.25809653802148 -0.1 47 8 2.19722457733622 1 3.25809653802148 0.1 47 8 2.19722457733622 1 3.25809653802148 0.3 47 2 1.01160091167848 1 3.2188758248682 -0.3 48 1 1.01160091167848 1 3.2188758248682 -0.1 48 0 1.01160091167848 1 3.2188758248682 0.1 48 0 1.01160091167848 1 3.2188758248682 0.3 48 102 3.63098547569503 1 3.09104245335832 -0.3 49 65 3.63098547569503 1 3.09104245335832 -0.1 49 72 3.63098547569503 1 3.09104245335832 0.1 49 63 3.63098547569503 1 3.09104245335832 0.3 49 4 1.70474809223843 1 3.46573590279973 -0.3 50 3 1.70474809223843 1 3.46573590279973 -0.1 50 2 1.70474809223843 1 3.46573590279973 0.1 50 4 1.70474809223843 1 3.46573590279973 0.3 50 8 2.32727770558442 1 3.2188758248682 -0.3 51 6 2.32727770558442 1 3.2188758248682 -0.1 51 5 2.32727770558442 1 3.2188758248682 0.1 51 7 2.32727770558442 1 3.2188758248682 0.3 51 1 2.07944154167984 1 3.55534806148941 -0.3 52 3 2.07944154167984 1 3.55534806148941 -0.1 52 1 2.07944154167984 1 3.55534806148941 0.1 52 5 2.07944154167984 1 3.55534806148941 0.3 52 18 2.63905732961526 1 3.04452243772342 -0.3 53 11 2.63905732961526 1 3.04452243772342 -0.1 53 28 2.63905732961526 1 3.04452243772342 0.1 53 13 2.63905732961526 1 3.04452243772342 0.3 53 6 1.79175946922806 1 3.71357206670431 -0.3 54 3 1.79175946922806 1 3.71357206670431 -0.1 54 4 1.79175946922806 1 3.71357206670431 0.1 54 0 1.79175946922806 1 3.71357206670431 0.3 54 3 1.38629436111989 1 3.46573590279973 -0.3 55 5 1.38629436111989 1 3.46573590279973 -0.1 55 4 1.38629436111989 1 3.46573590279973 0.1 55 3 1.38629436111989 1 3.46573590279973 0.3 55 1 1.70474809223843 1 3.25809653802148 -0.3 56 23 1.70474809223843 1 3.25809653802148 -0.1 56 19 1.70474809223843 1 3.25809653802148 0.1 56 8 1.70474809223843 1 3.25809653802148 0.3 56 2 1.83258146374831 1 3.04452243772342 -0.3 57 3 1.83258146374831 1 3.04452243772342 -0.1 57 0 1.83258146374831 1 3.04452243772342 0.1 57 1 1.83258146374831 1 3.04452243772342 0.3 57 0 1.17865499634165 1 3.58351893845611 -0.3 58 0 1.17865499634165 1 3.58351893845611 -0.1 58 0 1.17865499634165 1 3.58351893845611 0.1 58 0 1.17865499634165 1 3.58351893845611 0.3 58 1 1.09861228866811 1 3.61091791264422 -0.3 59 4 1.09861228866811 1 3.61091791264422 -0.1 59 3 1.09861228866811 1 3.61091791264422 0.1 59 2 1.09861228866811 1 3.61091791264422 0.3 59 ; title Fil fra Dale der han viser meg qpoint. Han sier dette fjerner problemet med sensitivitet. title "Fit same model but fix the number of quadrature points to 9"; title "Initial parameter set 1"; proc nlmixed data=nlmm qpoints=9; parms b0=0 b_base=0 b_trt=0 b_base_trt=0 b_age=0 b_visit=0 k=0.36 log_s2u1=0 log_s2u2=-1; eta_nb = b0 + b_base*base + b_trt*trt + b_base_trt*base*trt + b_age*age + b_visit*visit + u1 + u2*base; mean = exp(eta_nb); like = exp(lgamma(y+(1/k)) - lgamma(y+1) - lgamma(1/k) + y*log(k*mean) - (y+(1/k))*log(1+k*mean)); loglike = log(like); model y ~ general(loglike); random u1 u2 ~ normal([0,0], [exp(2*log_s2u1), 0, exp(2*log_s2u2)]) subject=subject; estimate "s2u1" exp(2*log_s2u1); estimate "s2u2" exp(2*log_s2u2); estimate "alpha" 1/k; run;