# EXAMPLE FOR USERS:

x.eg<-c(
0.529,0.645,0.517,0.228,0.554,0.338,1.005,0.358,0.416,0.227,0.918,0.439,
1.446,0.309,0.578,0.234,0.358,0.853,0.528,0.206,1.009,0.249,0.921,0.284,
2.240,0.300,0.330,0.114,0.759,0.215,0.049,0.383,0.653,0.242,1.149,0.827,
0.497,0.232,0.447,0.530,0.823,0.362,0.695,0.218,0.379,0.281,0.713,0.454,
0.378,0.666)

y.eg<-c(
10.530,14.684,13.826,,9.780,11.438,8.622,11.480,12.603,12.710,10.319,
12.793,11.301,13.283,12.865,9.780,9.780,,8.622,13.641,12.281,11.243,
9.780,11.920,12.193,,9.780,10.478,11.264,,9.780,8.622,,8.622,,8.622,
15.555,11.317,11.484,11.366,12.952,11.173,,8.622,8.622,,8.622,10.763,
9.780,10.758,13.618,,9.780,,9.780,,8.622,,8.622,12.291,12.573,,9.780)

ND.limiteg<-as.vector(rep(9.315,50))

QL.limiteg<-as.vector(rep(10.245,50))

flag.eg<-c(
"[no","[no","[no","<QL","[no","ND","[no","[no","[no","[no","[no","[no",
"[no","[no","<QL","<QL","ND","[no","[no","[no","<QL","[no","[no","<QL",
"[no","[no","<QL","ND","ND","ND","[no","[no","[no","[no","[no","[no",
"ND","ND","ND","[no","<QL","[no","[no","<QL","<QL","ND","ND","[no",
"[no","<QL")

names.vectoreg	<-("x.eg")
y.nameeg		<-("y.eg")

cbind(x.eg,y.eg,flag.eg,ND.limiteg,QL.limiteg)

x.eg  y.eg flag.eg ND.limiteg QL.limiteg 
0.529 10.53  [no   9.315    10.245  
0.645 14.684 [no   9.315    10.245  
0.517 13.826 [no   9.315    10.245  
0.228 9.78   <QL   9.315    10.245  
0.554 11.438 [no   9.315    10.245  
0.338 8.622  ND    9.315    10.245  
1.005 11.48  [no   9.315    10.245  
0.358 12.603 [no   9.315    10.245  
0.416 12.71  [no   9.315    10.245  
0.227 10.319 [no   9.315    10.245  
0.918 12.793 [no   9.315    10.245  
0.439 11.301 [no   9.315    10.245  
1.446 13.283 [no   9.315    10.245  
0.309 12.865 [no   9.315    10.245  
0.578 9.78   <QL   9.315    10.245  
0.234 9.78   <QL   9.315    10.245  
0.358 8.622  ND    9.315    10.245  
0.853 13.641 [no   9.315    10.245  
0.528 12.281 [no   9.315    10.245  
0.206 11.243 [no   9.315    10.245  
1.009 9.78   <QL   9.315    10.245  
0.249 11.92  [no   9.315    10.245  
0.921 12.193 [no   9.315    10.245  
0.284 9.78   <QL   9.315    10.245  
2.24  10.478 [no   9.315    10.245  
0.3   11.264 [no   9.315    10.245  
0.33  9.78   <QL   9.315    10.245  
0.114 8.622  ND    9.315    10.245  
0.759 8.622  ND    9.315    10.245  
0.215 8.622  ND    9.315    10.245  
0.049 15.555 [no   9.315    10.245  
0.383 11.317 [no   9.315    10.245  
0.653 11.484 [no   9.315    10.245  
0.242 11.366 [no   9.315    10.245  
1.149 12.952 [no   9.315    10.245  
0.827 11.173 [no   9.315    10.245  
0.497 8.622  ND    9.315    10.245  
0.232 8.622  ND    9.315    10.245  
0.447 8.622  ND    9.315    10.245  
0.53  10.763 [no   9.315    10.245  
0.823 9.78   <QL   9.315    10.245  
0.362 10.758 [no   9.315    10.245  
0.695 13.618 [no   9.315    10.245  
0.218 9.78   <QL   9.315    10.245  
0.379 9.78   <QL   9.315    10.245  
0.281 8.622  ND    9.315    10.245  
0.713 8.622  ND    9.315    10.245  
0.454 12.291 [no   9.315    10.245  
0.378 12.573 [no   9.315    10.245  
0.666 9.78   <QL   9.315    10.245


web.eg<-maxlikelihd.function(x.eg, y.eg, flag.eg, ND.limiteg, QL.limiteg, difference=0.01,names.vectoreg,y.nameeg)

[1] "CENSORING DETAILS:"
[1] "------------------"

Percentage ND censored:         20 % 
Percentage <QL censored:        20 % 
Total Percentage censored:      40 % 

[1] "TABLE OF RESULTS - Maximum Likelihood"

Outcome variable:  y.eg 

          Coefficient  Std error  z-value  P-value 
          -----------  ---------  -------  -------- 
(Initial)
Intercept    10.5343     0.445   23.6737         0 
x.eg          0.8126     0.675    1.2038    0.2346 

(Final)
Intercept    10.4056    0.4885   21.3003         0 
x.eg          0.9083    0.7305    1.2434    0.2137 

Initial Residual Variance:       3.2156
Final Residual Variance:         3.6776