|
|
@@ -28,5 +28,17 @@ def Q(t):
|
|
|
#screening intervals
|
|
|
t=numpy.linspace(0,9,10)
|
|
|
|
|
|
-three_parameter_screening.probabilty_of_clinical_incidence(beta,w,mu,t,Q)
|
|
|
+#generate incidence (interval cancers I_i)
|
|
|
+I=three_parameter_screening.probabilty_of_clinical_incidence(beta,w,mu,t,Q)
|
|
|
+
|
|
|
+#generate screen detection (D_i)
|
|
|
+D=three_parameter_screening.probabilty_of_preclinical_diagnosis(beta,w,mu,t,Q)
|
|
|
+
|
|
|
+#somehow generate a vector of data
|
|
|
+n,s,r=getData()
|
|
|
+
|
|
|
+#calculate likelihood
|
|
|
+alpha=[beta]
|
|
|
+L=three_parameter_screening.likelihood_function(D,I,n,s,r,alpha,beta)
|
|
|
+
|
|
|
```
|
