aboutsummaryrefslogblamecommitdiff
path: root/vignettes/mkin.html
blob: ea057d195a0faefe01b7acecbb956667b4e744fa (plain) (tree)
1
2
3
4
5
6
7
8
9

               
      


      
                        
                                          
                                                       
 
                                                                      
 
                                               
 
 
                                         
 









                                                                                      
   
   
         

                                                                      
 

                                                            
 

                                                                               

        
 






                                                                     
            
 




                       
   













                                                                          
 




















































                                                                                          
 




























                                                                                    
 

 





















                                      
 









                                                                         
 






































                          
 




                    
 


                           
 

                            
 

                
 


                          
 




                           
 

               
 


                    
 

                  
 

              
 

                           
 

                 
 




                          
 



                         
 

                          
 


                                                        
 

                     
 

                   
 



























































                                                              

        

 
 




       

 

                                                            
                                                                  
 
 
















                                                                              


      


                                                                      

                                          








                                                                        





























                                                                                                                                                                                                                                                                                                                                                                                                           
                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               


                                            









































                                                                                                                                             

                                                        
























                                                                                                                                                                                    
      
      

                                                 

                                                                     
    







                                                                                                                                         
     










                                                                                                                                           


                                                                    























                                                                         

                                                                                      



















                                                                        


                                                                                    




















                                                                                                                       




                                            



                                                               
      




                                                                                                                                                                    
      



                                                                                                                                                                 
      




                                                                          
      


                                                                                                                                                                 
      




                                                                                                                                                                                                                                                                  
      



                                                                                                                                 
      



                                                                                                                                    
      




                                                                                                      
      






                                                                         
      




                                                                                                       

      
      


 
                     
 
 




                                                                       
                                                                                               





                                                                 
<!DOCTYPE html>

<html>

<head>

<meta charset="utf-8" />
<meta name="generator" content="pandoc" />
<meta http-equiv="X-UA-Compatible" content="IE=EDGE" />

<meta name="viewport" content="width=device-width, initial-scale=1" />

<meta name="author" content="Johannes Ranke" />


<title>Short introduction to mkin</title>

<script>// Pandoc 2.9 adds attributes on both header and div. We remove the former (to
// be compatible with the behavior of Pandoc < 2.8).
document.addEventListener('DOMContentLoaded', function(e) {
  var hs = document.querySelectorAll("div.section[class*='level'] > :first-child");
  var i, h, a;
  for (i = 0; i < hs.length; i++) {
    h = hs[i];
    if (!/^h[1-6]$/i.test(h.tagName)) continue;  // it should be a header h1-h6
    a = h.attributes;
    while (a.length > 0) h.removeAttribute(a[0].name);
  }
});
</script>
<style type="text/css">detaiks.chunk-details > summary.chunk-summary {
text-align: right;
}
details.chunk-details[open] > summary.chunk-summary::after {
content: "Hide";
}
details.chunk-details[open] > summary.chunk-summary > span.chunk-summary-text {
display: none;
}
</style>

<style type="text/css">
  code{white-space: pre-wrap;}
  span.smallcaps{font-variant: small-caps;}
  span.underline{text-decoration: underline;}
  div.column{display: inline-block; vertical-align: top; width: 50%;}
  div.hanging-indent{margin-left: 1.5em; text-indent: -1.5em;}
  ul.task-list{list-style: none;}
    </style>



<style type="text/css">
  code {
    white-space: pre;
  }
  .sourceCode {
    overflow: visible;
  }
</style>
<style type="text/css" data-origin="pandoc">
pre > code.sourceCode { white-space: pre; position: relative; }
pre > code.sourceCode > span { display: inline-block; line-height: 1.25; }
pre > code.sourceCode > span:empty { height: 1.2em; }
.sourceCode { overflow: visible; }
code.sourceCode > span { color: inherit; text-decoration: inherit; }
div.sourceCode { margin: 1em 0; }
pre.sourceCode { margin: 0; }
@media screen {
div.sourceCode { overflow: auto; }
}
@media print {
pre > code.sourceCode { white-space: pre-wrap; }
pre > code.sourceCode > span { text-indent: -5em; padding-left: 5em; }
}
pre.numberSource code
  { counter-reset: source-line 0; }
pre.numberSource code > span
  { position: relative; left: -4em; counter-increment: source-line; }
pre.numberSource code > span > a:first-child::before
  { content: counter(source-line);
    position: relative; left: -1em; text-align: right; vertical-align: baseline;
    border: none; display: inline-block;
    -webkit-touch-callout: none; -webkit-user-select: none;
    -khtml-user-select: none; -moz-user-select: none;
    -ms-user-select: none; user-select: none;
    padding: 0 4px; width: 4em;
    color: #aaaaaa;
  }
pre.numberSource { margin-left: 3em; border-left: 1px solid #aaaaaa;  padding-left: 4px; }
div.sourceCode
  {   }
@media screen {
pre > code.sourceCode > span > a:first-child::before { text-decoration: underline; }
}
code span.al { color: #ff0000; font-weight: bold; } /* Alert */
code span.an { color: #60a0b0; font-weight: bold; font-style: italic; } /* Annotation */
code span.at { color: #7d9029; } /* Attribute */
code span.bn { color: #40a070; } /* BaseN */
code span.bu { color: #008000; } /* BuiltIn */
code span.cf { color: #007020; font-weight: bold; } /* ControlFlow */
code span.ch { color: #4070a0; } /* Char */
code span.cn { color: #880000; } /* Constant */
code span.co { color: #60a0b0; font-style: italic; } /* Comment */
code span.cv { color: #60a0b0; font-weight: bold; font-style: italic; } /* CommentVar */
code span.do { color: #ba2121; font-style: italic; } /* Documentation */
code span.dt { color: #902000; } /* DataType */
code span.dv { color: #40a070; } /* DecVal */
code span.er { color: #ff0000; font-weight: bold; } /* Error */
code span.ex { } /* Extension */
code span.fl { color: #40a070; } /* Float */
code span.fu { color: #06287e; } /* Function */
code span.im { color: #008000; font-weight: bold; } /* Import */
code span.in { color: #60a0b0; font-weight: bold; font-style: italic; } /* Information */
code span.kw { color: #007020; font-weight: bold; } /* Keyword */
code span.op { color: #666666; } /* Operator */
code span.ot { color: #007020; } /* Other */
code span.pp { color: #bc7a00; } /* Preprocessor */
code span.sc { color: #4070a0; } /* SpecialChar */
code span.ss { color: #bb6688; } /* SpecialString */
code span.st { color: #4070a0; } /* String */
code span.va { color: #19177c; } /* Variable */
code span.vs { color: #4070a0; } /* VerbatimString */
code span.wa { color: #60a0b0; font-weight: bold; font-style: italic; } /* Warning */

</style>
<script>
// apply pandoc div.sourceCode style to pre.sourceCode instead
(function() {
  var sheets = document.styleSheets;
  for (var i = 0; i < sheets.length; i++) {
    if (sheets[i].ownerNode.dataset["origin"] !== "pandoc") continue;
    try { var rules = sheets[i].cssRules; } catch (e) { continue; }
    var j = 0;
    while (j < rules.length) {
      var rule = rules[j];
      // check if there is a div.sourceCode rule
      if (rule.type !== rule.STYLE_RULE || rule.selectorText !== "div.sourceCode") {
        j++;
        continue;
      }
      var style = rule.style.cssText;
      // check if color or background-color is set
      if (rule.style.color === '' && rule.style.backgroundColor === '') {
        j++;
        continue;
      }
      // replace div.sourceCode by a pre.sourceCode rule
      sheets[i].deleteRule(j);
      sheets[i].insertRule('pre.sourceCode{' + style + '}', j);
    }
  }
})();
</script>



<style type="text/css">
/* for pandoc --citeproc since 2.11 */
div.csl-bib-body { }
div.csl-entry {
  clear: both;
}
.hanging div.csl-entry {
  margin-left:2em;
  text-indent:-2em;
}
div.csl-left-margin {
  min-width:2em;
  float:left;
}
div.csl-right-inline {
  margin-left:2em;
  padding-left:1em;
}
div.csl-indent {
  margin-left: 2em;
}
</style>

<style type="text/css">body {
background-color: #fff;
margin: 1em auto;
max-width: 700px;
overflow: visible;
padding-left: 2em;
padding-right: 2em;
font-family: "Open Sans", "Helvetica Neue", Helvetica, Arial, sans-serif;
font-size: 14px;
line-height: 1.35;
}
#TOC {
clear: both;
margin: 0 0 10px 10px;
padding: 4px;
width: 400px;
border: 1px solid #CCCCCC;
border-radius: 5px;
background-color: #f6f6f6;
font-size: 13px;
line-height: 1.3;
}
#TOC .toctitle {
font-weight: bold;
font-size: 15px;
margin-left: 5px;
}
#TOC ul {
padding-left: 40px;
margin-left: -1.5em;
margin-top: 5px;
margin-bottom: 5px;
}
#TOC ul ul {
margin-left: -2em;
}
#TOC li {
line-height: 16px;
}
table {
margin: 1em auto;
border-width: 1px;
border-color: #DDDDDD;
border-style: outset;
border-collapse: collapse;
}
table th {
border-width: 2px;
padding: 5px;
border-style: inset;
}
table td {
border-width: 1px;
border-style: inset;
line-height: 18px;
padding: 5px 5px;
}
table, table th, table td {
border-left-style: none;
border-right-style: none;
}
table thead, table tr.even {
background-color: #f7f7f7;
}
p {
margin: 0.5em 0;
}
blockquote {
background-color: #f6f6f6;
padding: 0.25em 0.75em;
}
hr {
border-style: solid;
border: none;
border-top: 1px solid #777;
margin: 28px 0;
}
dl {
margin-left: 0;
}
dl dd {
margin-bottom: 13px;
margin-left: 13px;
}
dl dt {
font-weight: bold;
}
ul {
margin-top: 0;
}
ul li {
list-style: circle outside;
}
ul ul {
margin-bottom: 0;
}
pre, code {
background-color: #f7f7f7;
border-radius: 3px;
color: #333;
white-space: pre-wrap; 
}
pre {
border-radius: 3px;
margin: 5px 0px 10px 0px;
padding: 10px;
}
pre:not([class]) {
background-color: #f7f7f7;
}
code {
font-family: Consolas, Monaco, 'Courier New', monospace;
font-size: 85%;
}
p > code, li > code {
padding: 2px 0px;
}
div.figure {
text-align: center;
}
img {
background-color: #FFFFFF;
padding: 2px;
border: 1px solid #DDDDDD;
border-radius: 3px;
border: 1px solid #CCCCCC;
margin: 0 5px;
}
h1 {
margin-top: 0;
font-size: 35px;
line-height: 40px;
}
h2 {
border-bottom: 4px solid #f7f7f7;
padding-top: 10px;
padding-bottom: 2px;
font-size: 145%;
}
h3 {
border-bottom: 2px solid #f7f7f7;
padding-top: 10px;
font-size: 120%;
}
h4 {
border-bottom: 1px solid #f7f7f7;
margin-left: 8px;
font-size: 105%;
}
h5, h6 {
border-bottom: 1px solid #ccc;
font-size: 105%;
}
a {
color: #0033dd;
text-decoration: none;
}
a:hover {
color: #6666ff; }
a:visited {
color: #800080; }
a:visited:hover {
color: #BB00BB; }
a[href^="http:"] {
text-decoration: underline; }
a[href^="https:"] {
text-decoration: underline; }

code > span.kw { color: #555; font-weight: bold; } 
code > span.dt { color: #902000; } 
code > span.dv { color: #40a070; } 
code > span.bn { color: #d14; } 
code > span.fl { color: #d14; } 
code > span.ch { color: #d14; } 
code > span.st { color: #d14; } 
code > span.co { color: #888888; font-style: italic; } 
code > span.ot { color: #007020; } 
code > span.al { color: #ff0000; font-weight: bold; } 
code > span.fu { color: #900; font-weight: bold; } 
code > span.er { color: #a61717; background-color: #e3d2d2; } 
</style>




</head>

<body>




<h1 class="title toc-ignore">Short introduction to mkin</h1>
<h4 class="author">Johannes Ranke</h4>
<h4 class="date">Last change 18 May 2023 (rebuilt 2024-04-27)</h4>


<div id="TOC">
<ul>
<li><a href="#abstract">Abstract</a></li>
<li><a href="#background">Background</a>
<ul>
<li><a href="#derived-software-tools">Derived software tools</a></li>
</ul></li>
<li><a href="#unique-features">Unique features</a></li>
<li><a href="#internal-parameter-transformations">Internal parameter
transformations</a>
<ul>
<li><a href="#confidence-intervals-based-on-transformed-parameters">Confidence
intervals based on transformed parameters</a></li>
<li><a href="#parameter-t-test-based-on-untransformed-parameters">Parameter
t-test based on untransformed parameters</a></li>
</ul></li>
<li><a href="#references">References</a></li>
</ul>
</div>

<p><a href="https://www.jrwb.de">Wissenschaftlicher Berater, Kronacher
Str. 12, 79639 Grenzach-Wyhlen, Germany</a><br /> Privatdozent at the
University of Freiburg</p>
<div id="abstract" class="section level1">
<h1>Abstract</h1>
<p>In the regulatory evaluation of chemical substances like plant
protection products (pesticides), biocides and other chemicals,
degradation data play an important role. For the evaluation of pesticide
degradation experiments, detailed guidance has been developed, based on
nonlinear optimisation. The <code>R</code> add-on package
<code>mkin</code> implements fitting some of the models recommended in
this guidance from within R and calculates some statistical measures for
data series within one or more compartments, for parent and
metabolites.</p>
<details class="chunk-details"><summary class="chunk-summary"><span class="chunk-summary-text">Code</span></summary>
<div class="sourceCode" id="cb1"><pre class="sourceCode r"><code class="sourceCode r"><span id="cb1-1"><a href="#cb1-1" aria-hidden="true" tabindex="-1"></a><span class="fu">library</span>(<span class="st">&quot;mkin&quot;</span>, <span class="at">quietly =</span> <span class="cn">TRUE</span>)</span>
<span id="cb1-2"><a href="#cb1-2" aria-hidden="true" tabindex="-1"></a><span class="co"># Define the kinetic model</span></span>
<span id="cb1-3"><a href="#cb1-3" aria-hidden="true" tabindex="-1"></a>m_SFO_SFO_SFO <span class="ot">&lt;-</span> <span class="fu">mkinmod</span>(<span class="at">parent =</span> <span class="fu">mkinsub</span>(<span class="st">&quot;SFO&quot;</span>, <span class="st">&quot;M1&quot;</span>),</span>
<span id="cb1-4"><a href="#cb1-4" aria-hidden="true" tabindex="-1"></a>                         <span class="at">M1 =</span> <span class="fu">mkinsub</span>(<span class="st">&quot;SFO&quot;</span>, <span class="st">&quot;M2&quot;</span>),</span>
<span id="cb1-5"><a href="#cb1-5" aria-hidden="true" tabindex="-1"></a>                         <span class="at">M2 =</span> <span class="fu">mkinsub</span>(<span class="st">&quot;SFO&quot;</span>),</span>
<span id="cb1-6"><a href="#cb1-6" aria-hidden="true" tabindex="-1"></a>                         <span class="at">use_of_ff =</span> <span class="st">&quot;max&quot;</span>, <span class="at">quiet =</span> <span class="cn">TRUE</span>)</span>
<span id="cb1-7"><a href="#cb1-7" aria-hidden="true" tabindex="-1"></a></span>
<span id="cb1-8"><a href="#cb1-8" aria-hidden="true" tabindex="-1"></a></span>
<span id="cb1-9"><a href="#cb1-9" aria-hidden="true" tabindex="-1"></a><span class="co"># Produce model predictions using some arbitrary parameters</span></span>
<span id="cb1-10"><a href="#cb1-10" aria-hidden="true" tabindex="-1"></a>sampling_times <span class="ot">=</span> <span class="fu">c</span>(<span class="dv">0</span>, <span class="dv">1</span>, <span class="dv">3</span>, <span class="dv">7</span>, <span class="dv">14</span>, <span class="dv">28</span>, <span class="dv">60</span>, <span class="dv">90</span>, <span class="dv">120</span>)</span>
<span id="cb1-11"><a href="#cb1-11" aria-hidden="true" tabindex="-1"></a>d_SFO_SFO_SFO <span class="ot">&lt;-</span> <span class="fu">mkinpredict</span>(m_SFO_SFO_SFO,</span>
<span id="cb1-12"><a href="#cb1-12" aria-hidden="true" tabindex="-1"></a>  <span class="fu">c</span>(<span class="at">k_parent =</span> <span class="fl">0.03</span>,</span>
<span id="cb1-13"><a href="#cb1-13" aria-hidden="true" tabindex="-1"></a>    <span class="at">f_parent_to_M1 =</span> <span class="fl">0.5</span>, <span class="at">k_M1 =</span> <span class="fu">log</span>(<span class="dv">2</span>)<span class="sc">/</span><span class="dv">100</span>,</span>
<span id="cb1-14"><a href="#cb1-14" aria-hidden="true" tabindex="-1"></a>    <span class="at">f_M1_to_M2 =</span> <span class="fl">0.9</span>, <span class="at">k_M2 =</span> <span class="fu">log</span>(<span class="dv">2</span>)<span class="sc">/</span><span class="dv">50</span>),</span>
<span id="cb1-15"><a href="#cb1-15" aria-hidden="true" tabindex="-1"></a>  <span class="fu">c</span>(<span class="at">parent =</span> <span class="dv">100</span>, <span class="at">M1 =</span> <span class="dv">0</span>, <span class="at">M2 =</span> <span class="dv">0</span>),</span>
<span id="cb1-16"><a href="#cb1-16" aria-hidden="true" tabindex="-1"></a>  sampling_times)</span>
<span id="cb1-17"><a href="#cb1-17" aria-hidden="true" tabindex="-1"></a></span>
<span id="cb1-18"><a href="#cb1-18" aria-hidden="true" tabindex="-1"></a><span class="co"># Generate a dataset by adding normally distributed errors with</span></span>
<span id="cb1-19"><a href="#cb1-19" aria-hidden="true" tabindex="-1"></a><span class="co"># standard deviation 3, for two replicates at each sampling time</span></span>
<span id="cb1-20"><a href="#cb1-20" aria-hidden="true" tabindex="-1"></a>d_SFO_SFO_SFO_err <span class="ot">&lt;-</span> <span class="fu">add_err</span>(d_SFO_SFO_SFO, <span class="at">reps =</span> <span class="dv">2</span>,</span>
<span id="cb1-21"><a href="#cb1-21" aria-hidden="true" tabindex="-1"></a>                             <span class="at">sdfunc =</span> <span class="cf">function</span>(x) <span class="dv">3</span>,</span>
<span id="cb1-22"><a href="#cb1-22" aria-hidden="true" tabindex="-1"></a>                             <span class="at">n =</span> <span class="dv">1</span>, <span class="at">seed =</span> <span class="dv">123456789</span> )</span>
<span id="cb1-23"><a href="#cb1-23" aria-hidden="true" tabindex="-1"></a></span>
<span id="cb1-24"><a href="#cb1-24" aria-hidden="true" tabindex="-1"></a><span class="co"># Fit the model to the dataset</span></span>
<span id="cb1-25"><a href="#cb1-25" aria-hidden="true" tabindex="-1"></a>f_SFO_SFO_SFO <span class="ot">&lt;-</span> <span class="fu">mkinfit</span>(m_SFO_SFO_SFO, d_SFO_SFO_SFO_err[[<span class="dv">1</span>]], <span class="at">quiet =</span> <span class="cn">TRUE</span>)</span>
<span id="cb1-26"><a href="#cb1-26" aria-hidden="true" tabindex="-1"></a></span>
<span id="cb1-27"><a href="#cb1-27" aria-hidden="true" tabindex="-1"></a><span class="co"># Plot the results separately for parent and metabolites</span></span>
<span id="cb1-28"><a href="#cb1-28" aria-hidden="true" tabindex="-1"></a><span class="fu">plot_sep</span>(f_SFO_SFO_SFO, <span class="at">lpos =</span> <span class="fu">c</span>(<span class="st">&quot;topright&quot;</span>, <span class="st">&quot;bottomright&quot;</span>, <span class="st">&quot;bottomright&quot;</span>))</span></code></pre></div>
</details>
<p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAwAAAAKgCAIAAACjiI2DAAAACXBIWXMAAA7DAAAOwwHHb6hkAAAgAElEQVR4nOzdeVxM+/8H8DNtWqZp1x7tC2nSTTshaxKVlChCSFFcu9zsrvVLt2vNnksk0eJeZMmeVEIlRJMRKa3am98f8/12+4WQmc4083r+8X2c+fTpc16Tuef7PmfO+XwoLBaLAAAAABAkQmQHAAAAAOhqKIAAAABA4KAAAgAAAIGDAggAAAAEDgogAAAAEDgogAAAAEDgoAACAAAAgYMCCAAAAAQOCiAAAAAQOCiAAAAAQOCgAAIAAACBgwIIAAAABA4KIAAAABA4KIAAAABA4KAAAgAAAIGDAggAAAAEDgogAAAAEDgogAAAAEDgoAACAAAAgYMCCAAAAAQOCiAAAAAQOCiAAAAAQOCgAAIAAACBgwIIAAAABA4KIAAAABA4KIAAAABA4KAAAgAAAIGDAggAAAAEDgogAAAAEDgogPjZy5cvx40b5+fnt3fvXrKzAAAA8BAKi8UiOwNwy+LFi729vc3NzQcOHHjjxg2y4wAAAPAKXAHiZ6tXr6bT6UVFRYqKitzby/Xr18ePH9+Vv/iTPn78+OrVq67fLwAA8A4UQPxMQkLizJkzQUFBf/75J9lZfsr+/fuNjIxUVFSmTJlSW1vb2r527VrR/xEREbGwsCAIYtu2bVpaWhoaGqGhoSwW6/Xr1zY2NkZGRrt27WL/1sqVK+vr68l5JwAAwBtQAPGzxMTEu3fvxsXFqaiofE//xsbGL2533JPb0tPTN2zYcO3atWfPnpWVlW3evLn1R2FhYY3/s2DBgvnz52dlZW3bti0jIyMvL+/mzZtxcXEbN24MDw9/8uTJkSNHKisrnz17RhCEoaFhl+UHAAAehAKIfyxZssTCwqKuro4giIiIiICAgLi4uLy8PE9PzylTprTtGRERoaOjo6+vv3DhwpaWllu3bvn4+Hh7e69Zs6btNkEQ69ev19XV1dHRWbx4cbuenwdoN+zIkSPPnTvH/pGlpeWNGzfadfjO91VYWDh37lwVFRUajebu7v78+fPP+zx48CA3N9fX11dISEhUVFRCQkJMTExcXFxISEhGRubp06cMBqOurk5UVHT9+vUrVqz4oT8sAADwIRbwET8/vwMHDtTW1urq6hYVFX2xz7Vr18zMzN6/f19dXe3l5bVu3bqbN2/KyMjk5uayWKy22xcvXjQ0NKyoqKirq7O3tz969Gjbn7YdcNy4cZ8Pe+jQoSlTprBYrBcvXujo6Fy9erVdB/Yvto6TkJBg1UZAQEC75KWlpZaWlnFxcZ+/KQcHhydPnrC3f/31VwkJCRkZGfbgFRUVs2bNcnJy+ueff9LS0pYtW9bpPy8AAPANEbILMOCk2bNnL1y4sLKy0sPDQ11d/Yt9rly5UlFR4enpSRBESUlJVVWVo6MjnU5v/VaodTslJcXDw4NGoxEE4e3tnZKSoqOj07Znx8NGR0cvXbq0sbHx1KlTfn5+KSkp7TrY29u3HcHZ2dnZ2flrby06Ojo8PDw8PHzcuHHtfpScnNyzZ08TExOCIG7dupWQkHD9+nVpaekZM2ZER0f7+Pjs2bOH3XPChAn79++/cOFCZGSkqanpunXrevTo8T1/WAAA4DMogPiKtbV1aWlpREREenr61/pQqdRZs2YtXbqUIIi6urqGhobs7GxpaenWDq3bLBaLQqGwt4WFhZubm9v+9JvD0mg0a2vrq1evxsTEnD179vTp0+06ZGRktB0hMTFx7dq1rS/NzMzY0xe1tLRMnTq1oqIiNTX1izczHT9+3N3dnb2dkJDg5uZmaWlJEERAQMC5c+d8fHzYP/rnn38GDBjQ3NwcHh4eFxd34MCByMjIBQsWdPD3BAAAfoV7gPhN7969bW1t5eTkvtZh+PDhR44cKSkpqa+vd3d3v3Xr1td6DhkyJDY2trq6uqGhITo6esiQIR3s94vDenl5bdy4UU5OTltb+5v7dXZ2vttG6+SN8fHxTCbz3LlzrdUPi8XKy8tramoiCKK+vv7ixYsjR45k/8jU1DQpKYl9hSk2NrZfv36tvxIZGRkcHFxTUyMnJ6epqWlqalpaWtrxHxMAAPgVCiC+cuvWrdLS0pSUlA4e86bT6SEhITY2Nrq6urq6uqNGjfpaz5EjR06aNMnc3NzY2HjAgAG+vr4d7PqLw7q4uNy/f9/Pz++H9ttOampqamoqjUaTlpaWlpaeOnVqQ0ODkZHRmzdvCIL4559/zMzMWgs+b2/vMWPGWFhYGBkZKSsrL1q0iN0eExMzduxYcXFxLS0te3v7Pn36REREzJs37zszAAAAn8FM0HzFzs5u06ZNO3funDhx4oQJE8iOAwAAwKNwBYh/nDt3jkajOTg4TJ8+vfW2XwAAAPgcCiA+0dzcvHz58nXr1hEEMXLkSCUlJTqdXlJSQnYuAAAAXoSvwAAAAEDg4AoQAAAACBwUQAAAACBwUAABAACAwEEBBAAAAAIHBRAAAAAIHC6uBdbS0pKVlVVcXNzU1KSurk6n04WEUG8BQFeoqKig0WgUCiU9PT0nJ8fCwsLY2JjsUADAQ7j1GPzVq1fnzJmjqqrKXpOcyWQyGIyoqKiBAwd+5wi3bt26cOECN7IBQOdoaWkFBgaSneLb9u3bt3379uzs7C1bthw4cMDBweHatWsrV66cOXPmd46A4w8Ar+H48YdbBZCpqem5c+d0dXVbWxgMhoeHx717975zhMWLFz958uT7CyYA4KqKiopjx44xGAyyg3xbr169srOzaTSapqZmVlaWvLx8SUmJra1tfn7+F/uPGjWq3aGpqqpKXV19zpw5XZIXAL6BG8cfbn0F1tzcrKqq2rZFWVn5R4stR0fH1sUsAYBcRUVFx44dIzvFd1FSUmpoaCAIQlZWlv3Nu7CwcAf94+Liamtr27Y4OTmpqqouWbKEqzkB4Dtx4/jDrQIoKCjIwsLC1dVVQ0ODQqEwmczz588HBQVxanwGg7F27dqioiIrK6tevXppa2sPGjSIU4MDQLe2fv16W1vbsWPHGhgYODo6Dh8+PDExMTg4+Gv9xcXFxcXF27aIiHDx/kgA4AXc+o88MDDQxcUlOTmZyWQSBKGlpZWUlKSpqcmRwRkMhqGhoaioqJiY2MWLF7W0tGRlZXV0dM6cOYP7rAFgxIgR9+/fj4+Pl5WVNTIyUlJSio2NNTIyIjsXAPAQLp7lqKurW1patj4Fxr4bmiN+/fVXRUXF9PR0IyOjEydOTJo0KScnZ8iQIcnJyc7OzpzaCwB0X7Kysn5+fmSnAADexa0C6OefAuvAy5cvzc3NCwsLtbW1vby8fHx8Xr16ZWdnl5ubiwIIAAAAvolbBdC8efMSExO//ymwGzduvHv3rm1LRkZGnz59vtjZ2tr66NGj8+fPf/z4cUhICIVC6d2795UrV7Zs2cLBtwDQgSdPniQkJJCdggMGDRpkbW1NdgoAgK7GK0+BXblyJScnp23Lw4cPP3369MXOv//++7lz50aMGNHc3Lxr1y5DQ0NjY+MRI0Y4OTlxJDzAN505c+bSpUv29vZkB/kpjx49ys/PRwEEAAKIV54CW716dbsWKysrOTm5L3aWlJR8+fLl+fPnb9y4cfLkyXXr1hkbG5uYmHDyDQB0iMViOTk5hYeHkx3kp0RFRd25c4fsFAAAJOiWT4ERBCEqKuru7u7u7l5cXPz48WN3d3dOjQwAAAB8j4tPgWlqagYEBHBvfLZt27aZmZn5+Pjo6elxe18AAADAH7r3rDk3btywt7dvbGw0MTH5888/yY4DAAAA3QO3rgDt2bPn+vXrn7f/9ddfnNpFfX29t7f30aNHHRwcTExMVq9ebW1t3b9/f06NDwAAAPyKWwXQlClT7t27V1lZyb21o3NycpSVlYcOHUoQxJ49e9zd3S9fvowCCADge+Tk5PznP/95//790KFDZ8+ejdU/QNBw6xMvJSUVGhp68uRJdoHCDcrKykwms66uTlxc3MnJSVpa+u7du1zaFwBAV6qpqTlx4sT79+8HDx5sa2vL8fFfvHgxZMiQX3/9ddSoUZGRkU+fPsVdBCBouFjy9+vXr1+/ftwbX1VV1dnZediwYRMmTEhPT5eWlk5NTc3LyzM0NOTeTgEAuK2iomLAgAHm5uZ6enpTpkwJDg4OCQnh7C6io6OnTZu2cOFCgiCcnJzU1NR27NjRo0cPzu4FgJd175ug9+/fHxgY+OLFC3Nz87S0tFWrVs2aNauD6RYBulJdXV1hYWFLSwvZQQiCICorK/fv3092CvguBw8etLe3Z09yduPGjfDw8ObmZs7uoqqqSkFBgb0tJSUlIiJSV1fH2V0A8LjuXQAJCQl5e3vv3LkzJCSESqXOnTu3rq5u27ZtqIGAdFu3blVRUXFwcNDW1r5x4wZX99XY2PjNPlVVVUePHuVqDOCUwsJCU1NT9ra6urqYmFhZWRlnd+Hs7Lxnz57Hjx/X1dWtWbOmX79+MjIynN0FAI/r3gVQO1u3bn3y5MmSJUt69+59+/ZtsuOA4Lpz587evXtzc3Nfv3596NAhb2/v76lR2tq5c+fEiRMNDAz09PQ2bdpEEASLxZoxY4a+vn6vXr3Yy979/fffs2fPHjZsWEJCQnh4uJ6eHp1OP3z4MEEQV65ccXNzs7OzMzExWb58OUEQixcvfvz48fr16zn/boHTbGxsTp8+XVNTQxBEYmKijIyMkpISZ3fh6Oi4dOlSZ2dneXn59PT048ePc3Z8AN7HP7f9p6amRkVFPX/+fOfOnbdu3fLy8iooKBAWFiY7F/C5nTt3Hjt2rF1jcXFxc3PzmDFj2C8/fvxIp9MlJCTa9hEVFd21a5elpeXXRr5//352djZBEObm5sOGDaNSqc3Nzc+ePausrDQ2Nl60aBFBEBcuXMjMzExLS7t58+ajR4/q6+utra3pdDpBEJcuXSosLKRSqfr6+qGhoZs3by4sLFyxYgVn3z5ww4QJE27cuNGrVy8VFZXq6upTp05xYy/Tp0+fPn06N0YG6Bb4qgCaMGGCsrJyWFiYmZkZi8V6/vw5bogGbvP29v58SdSkpKSUlJStW7cSBPHp0ycXF5eIiIjPv2Lo27dvByOPHj2aSqUSBDF27NjU1NSQkJDly5cfPHgwIyPj48eP7D4ODg5KSkqpqanFxcUTJ04kCKK2tjYjI0NLS2vw4MHs1fT09fUrKyvFxcU59I7hv1gs1vnz5+/fv9+rV68pU6a0K3B/BoVC+eOPP1atWlVSUmJgYCAqKsqpkQGgFf8UQCoqKklJSYcOHSoqKpo5c+bSpUvxRAN0gZ49e/bs2bNdo7Gx8enTp7dv306n00+fPu3r6ztkyJAfHZlCobA3hIWFm5ubU1JSQkND58+fP3v27ISEBPaPpKWlCYKQkJBYuHDhtGnTCIKorq7u0aPHjRs3Wm9xhX379nFjWZ6goKAHDx64urpeunRp7969t27d4myV+cWPFgBwCv8UQC4uLnPnzk1LSzM0NLx586aSktLq1asPHTpEdi4QRJKSknfu3Dl06FBhYeHy5ctdXV07MUhycjL7LpD4+PijR49ev3594sSJ/v7+d+7cef/+fduHywYPHrxixQofH5/a2lozM7PU1NQvDsjxJ4l4VlJSUklJSevLNWvWsE+H/Pz8OLWLd+/enT59+tWrV5KSkgRBjB49+vz5856enpwaH0BwvHz58vjx4/X19e7u7l05mzH/FEDJyclOTk5OTk4MBsPX13fu3LnXr19PTk4eNWoU2dFAEElJSQUFBf3MCHQ63cnJqbS01M/Pz8rKSkNDY+rUqfHx8f369Rs8ePDKlSsHDRrE7ung4DBs2LB+/fo1NjauWLFCU1Pz2bNn7UZTVFSsrq5etWrVmjVrfiZVt/DkyZOVK1d6eXnJy8sTBFFXV5eZmUl8vQDau3dvQUFB25aioqKOL+cwmUxNTU129UMQhJGR0Zs3bziTHkCQZGZmjhgxYtq0aVJSUi4uLn/++Wfnzhg7gX8KoMLCQjqdPn/+fPbLDRs2LF68ePbs2Y8ePcLjndAd2drasuepY1NXV7906VK7PiNGjGBvhIWFhYWFtbYPHTq0dRL21t969OgRF+PykkWLFtna2i5btszX13fo0KGXL1/esWNHB/2lpaXZ90u1EhISav0K8ouMjY2ZTObNmzft7e2Li4vj4+NPnjzJmfQAgmTr1q1r165lf0ltb2+/ZMkSFEA/zMbGZv78+QsXLpSVlb1582Z5efnkyZPv378fEhKCL8IABI2dnV1CQkJwcHBCQkJ9fX3HnSdNmtSu5ezZs+w70L9GXFz82LFjXl5e4uLiZWVlK1eu7OCBPgD4mrdv3+rr67O3DQwMmExml+2afwqgoUOHjh8/XldXV01Nrays7NixYyIiIlu3bqXT6bGxse7u7mQHBPgBc+bMwXyeP4lGox05cuTkyZO5ubncGN/JyenVq1evX79WVVVt/S4MAH7IwIED9+3bZ2trKyoq+scff7R+s98FuFgAtbS0ZGVlFRcXNzU1qaur0+l0ISHuzru4du3a0NBQJpOpr6/PvudRUlLyyJEjbm5u9vb2ysrKXN07AAeJiYmRHYFPeHl5eXl5cWlwERERXV1dLg0OIAiWLFni5+enoqIiKipqYmISExPTZbvmVgF09erVOXPmqKqqqqurEwTBZDIZDEZUVNTAgQO5tEc2eXl59m2PrWxsbPz8/ObMmePm5nbp0iVFRcW5c+fq6OhwNQYAAAB8k7i4+KlTpz5+/NjU1MTxGc87xq0CaN68eYmJiW3PjRgMhoeHx71797i0xw6sXr1aS0vrzp07JiYmr169sra2vn//fu/evbs+CQAAALTT7imErsGt76Sam5tVVVXbtigrK5N1T4OIiEh5efmHDx8GDhxoa2tbW1u7a9cuUpIAAAAAL+DWFaCgoCALCwtXV1cNDQ0KhcJkMs+fP/+T06J0WnV1dWNj44YNG2JjY2/fvp2WloalUuHnffz48eXLl2Sn+CkfPnwgOwIAADm4VQAFBga6uLgkJyezH2nT0tJKSkrS1NT8Wv/o6OiioqK2LUwmk1Nr6wgJCbHvv+7Zs+fChQufPn1qYWHBkZFBYPXu3fvYsWOtS1J0X1gOEwAEExefApOXl/fx8ZGSkmIymampqWVlZR0UQLW1ta3rO7K1tLRw6iszaWnpfv36HTx4sKio6OLFi3JycjNmzODIyCCwpk6dOnXqVLJTAABAJ3GrADp27FhwcLCIiEh4ePiOHTtsbGwWL168aNGir30L9nlFcvXqVfZCjxxx+vTpKVOmvH37VlxcvKGhwcrKilMjAwAAQLfDrQIoPDw8NzdXXFxcU1Pz4sWLdnZ2Hz58GDBgAFm3Aenq6t6+fbu6ulpCQiI8PHzSpEmXL1/m9rxEAAAAwJu4VQG0tLTIy8vTaDQpKSlFRUWCIKhUasdr63QBKpUqLCz822+/NTU1bd68mdwwAAAAQBZuXQHy9PS0srISExMbNmzY1KlTfXx8kpOTR44cyaXd/RAREZETJ04MGDDAxsamK2fdBgAAAB7BrQLo999/d3Z2FhYWtrOzu3LlSnJy8rhx4/z9/bm0ux+loaERHR09adKktLQ0NTU1suMAAABAl+LiU2Ctq14MHTp06NCh3NtR5wwePHj27NkTJky4du2aqKgo2XEAAACg6wj0XcArVqyQkZFZsWIF2UEAAACgSwl0ASQkJHT8+PHY2NhTp06RnQUAAAC6Dhe/AusW5OXl4+Lihg4dqq+v379/f7LjAAAAQFcQ6CtAbP369du7d6+HhwfWRQIAABAQKIAIgiDc3Nw8PDy8vb2bmppaG5OSkuzt7Q0NDWfPnl1WVkZiPAD4Uenp6fHx8RUVFa0tfLBwGwBwEAqg/9q0aZOEhERoaCj7ZXp6ekBAwMqVKy9cuCAkJOTr60tuPAD4fhs3bnRzc4uJiaHT6RkZGezGmTNnkpsKAHiKoN8D1EpISCg6OtrW1nb37t1z5sw5e/bsnDlz2DM3RkREKCgoVFVVcXBtMgDgnsjIyKysLAUFhezsbHd39/T09I7/4/306VN9fX3blqamJk4txgwAvAkF0L+kpaXj4+Pt7OwMDAyEhYVbvw5raWkhCAILhwF0F9LS0lJSUgRBmJqazpkzJyQkJCoqqoP+np6et2/fbttSVVVF+tI9AMBV+D/1/0dHR+fUqVOTJk2ysLDYu3fvX3/9lZaWNm3atBEjRrCPpwDA+7y9ve3s7Pbt20cQREhISFlZmZeX16dPn77WPyEhoez/69+/v4qKShdGBoCuhitA7Q0cOHDXrl2hoaF79uzZt2/f+/fvhw4diskSAbqRVatW2dnZlZaWEgRBoVDOnDkTExMjJydHdi4A4CEogL5g4sSJeXl569evv3r1qqSkJNlxAPjfnDlzvti+e/fuzg3YdvkdYWFhb29vb2/vzg0FAHzpCwVQZmZmZmamj49PcXGxpqZm12fiBWFhYQUFBV5eXnFxccLCwmTHAeBzY8aMITsCQLfR0tKC21J/Xvu/4OHDhydOnLhq1ar6+noHB4cDBw6QEot0FArlwIEDQkJCs2fPJjsLAP9z/syoUaNu3bpFdi4A3nL69GltbW1xcXFra+usrCyy43Rv7QugHTt23L1718zMjEqlZmZmbtq0qdNDt7S0ZGRkJCcnX7hw4eHDh+xnqboRYWHh6OjorKys8PBwsrMACITNmzerqalJSEj07t2bSqXi+A7Q1tOnT0NDQ2NjY2tra4OCgsaPH9928l74Ue0LoJqaGnFxcfa2pKSkqKho58a9evWqiYnJggULoqOjT5069euvvxoaGt64ceOnwnatysrKQ4cODRo0KCoqKjIykuw4APzvjz/+yMnJ8ff3T0xMvHLlirq6OtmJvqyqqmr37t3h4eHd65gG3d2VK1fc3d379+8vLCw8efJkaWnpnJwcskN1Y+0LIB8fn7FjxxYUFBw8eNDJycnHx6dz486bNy8xMfHq1avHjx8/fvx4SkpKSkrKokWLfjpwF/n48WP//v3T0tKkpaXFxMTCwsIOHTpEdigAPtejRw8ajWZubn7v3j0bGxvevAJUXl7ev3//27dvUyiUGTNmbNmyhexEIChkZGTev3/P3m5qaiotLZWVlSU3UrfW/ibo1atXJyYm3rx5882bNxs2bLC3t+/cuM3Nzaqqqm1blJWVu9HMqgcPHhw6dOjevXsJgpg5c6aBgcHKlStlZGTc3NzIjgbAtwYOHOjv779s2bIJEya8e/dORkaG7ERfcOjQIUdHx/379xMEERAQYGRkFBoaKiKCJ2qB68aOHRseHr5kyRILC4uTJ0/+8ssvAvugEke0/4+2tLTU2tra2tq69aWCgkInxg0KCrKwsHB1ddXQ0KBQKEwm8/z580FBQT+bt6u8efPG2NiYva2qqiohIXH48OHJkyfLysoOGTKE3GwA/OrAgQPPnj0zMDDYvn17eno6b152LSoqMjExIQiCxWKxDw5lZWU9e/bs+LdaWlpev34tJyeHU3boNFlZ2du3b+/YsSM2NtbBwSEwMJDsRN1b+wKIvfoVQRAVFRX5+flmZmaZmZmdGDcwMNDFxSU5OZnJZBIEoaWllZSU1EGtunHjxlevXrVtefXqFYknVXZ2dlu2bPH396fRaHFxcfLy8sOGDYuLixs/fnxMTMygQYPICgbAx8rKyhQVFUtLS+l0Op1OJzvOl9na2i5atGjz5s1lZWUmJiYSEhLfrH6ysrI8PT0bGxs/fvw4ZcqUnTt3Yp0N6BwVFZXff/+d7BR8on2FkZaW1rqdnZ39M6ug02i0mTNnUiiU9PT0nJyc6urqDjpbWFi0u9R06dIlCQmJTu/9J7m7u9++fbtXr16KioosFuvUqVMEQdja2p44ccLT0/PChQsDBgwgKxsAv+LUCRhXycvLf/z4saWlpXfv3oWFhWpqat/8FV9f3/DwcG9v76qqqlGjRp08eRKzMgKQj9WhPn36dNzha/bu3WtoaNjQ0LB+/XptbW1fX18tLa19+/Z9/wgDBgxwdnbu3N45paysLD8/n70udKvExERlZeX79++TlQqAFAwGQ0NDo8t29+jRIzqd3mW7a6eD48/ixYs3b95cWlqam5vb2NiooKBQWlrawVAVFRUyMjKtLyMjIwMDAzkcF4DfceP40/4KUGhoaOt2bm6uoaFh5+qq9evXZ2dni4qK7t69OysrS15evqSkxNbWdubMmZ0v1rqcnJzc5+sHjR49+vDhwy4uLvHx8VZWVqQEA+B7pqamjY2NZKf4AvaTOPLy8vLy8tXV1Q0NDR0vmEOlUoWEhIqKijQ0NAiCyMnJYW8AALnaF0BtH/tydnYePHhw58ZVUlJqaGggCEJWVpY9Yzc/LSgxcuTIgwcPurq6ogYC4CBOnYBx1eTJk62trWVlZfX09Pbu3evr69s6d9oXCQkJrVq1avDgwb6+vq9evUpJSWl7pwEAkOXfAujevXsEQbQ7NXnw4EHn/g9+/fr1tra2Y8eONTAwcHR0HD58eGJiYnBw8E/G5R2jR48+dOjQ2LFjz5w54+DgQHYcAH7AqRMwrtLS0kpNTd25c2d2draXl9f06dO/+SshISH9+/e/dOlSv379tm3bhgfBAHgBhfW/uXnYN+UxmczU1NR+/fpRqdR79+6xH3rq3NDl5eXx8fEMBqO2tlZJSWnkyJFGRkbf/+tWVlZKSkoJCQmd23vXuH79+sSJE48ePTp8+HCyswBwV1FRkY2NDYPB4Mbg7BOwz5F1hbVbHH8ABAc3jj//XgH666+/CIIYNmxYTk4O+8pzQUFBQEBAp4eWlZX18/P7+Yi8bNCgQbGxsW5ubvv37x87dizZcQC6q//85z8ER0/AAAA61v4eoNevXxsYGLC3e/fuzZ7FBzpgZ2eXlJTk4uJSXl7+M7MGAAgyjp+AAQB0rH0BZGNj4+Hh4ePjw2Kx/vrrL0tLS7QSkdEAACAASURBVFJidS8WFhYpKSkjR4589+5dN1rvDIDX8OsJWF1d3datWy9duqSoqLhw4UJbW1uyEwHAZwXQvn379uzZc+LECSEhoUGDBuEM7DsZGRmlpqaOGDGitLR048aNmOYVoBM4eALW0tKSlZVVXFzc1NSkrq5Op9PZj6OSIjAw8OPHj+vWrXv58qWHh8fFixf79etHVhgAYPu3AAoODl67dm1YWBhBEOx1TJ89e/brr79GRESQlq5b0dTUvHXrlqurq6en57Fjxzp+MhYAPsepE7CrV6/OmTNHVVVVXV2dIAgmk8lgMKKiogYOHMjRvN+loaHhzJkz7969k5CQcHBwKCoqOnnyJAogANL9WwA5OTmJi4u3TkUPnSAnJ/f3339PmjTJxcUlNjaWRqORnQige+DsCdi8efMSExN1dXVbWxgMhoeHx9eeNVuxYkV+fn7blpycnBcvXnh6erbGY892kZiYeOTIkbaxv9l+8ODB2tpaX19fCoUSHBzMXjy1E+OgHe2C3L5nz56KigqCo/4tgFxdXQmC0NXVvXXr1vTp01euXHnt2rVNmzZxdn98T0JC4syZMyEhIfb29gkJCVpaWmQnAugGOHsC1tzczK6iWikrK7dO+fG5UaNGtVt7NS0tTVpaesKECeyX+vr67I2+ffu2Nn5nu5eX1+vXr6uqqsaNG1dWVrZr167jx49ramr+6DhoR7sgtzs7O3/tBKbz2i2NMXz48JMnT967d2/QoEGpqakjR47k7NIb348X1gL7GXv37tXQ0MCSYcA3umAtsJycnAMHDrBYrBUrVtjZ2aWmpnZunMjISCMjoyVLlkRERPzxxx/Lly/v27fvnj17vn8Ezh5/Pnz4MG3aNA0NDXNz85iYGE4NCyA4uHH8aX9XYHFx8cSJE8+ePTt58mR7e/uXL19yuOASGAEBAbt37x4zZgx7Jfmvyc3NTU1Nraqq6rJgADxr/vz5VCr1/v37N2/e3LRp0/r16zs3TmBg4D///KOjo/Phw4eSkhItLa2kpKRZs2ZxNu33U1BQOHjwIIPBePjwYdsTXAAgUfunwAwNDRcuXBgTE5OWlrZjxw7cxfIzxowZc/ny5XHjxmVkZGzYsKHdQyjNzc3e3t5paWmampr5+fnR0dFDhgwhKyoAL2CfgC1duvTnT8A0NTXxECsAdKD9FaDIyEh1dfVz586pqKiUlpayZyeDTjM1Nb179+6dO3fc3d3bXeY5duxYRUVFfn7+jRs3YmJipk2bRlZIAB7BPgGLjo4eM2YMTsAAgKvaF0BKSkpDhgzJzs5ubGycNWuWnp4eKbH4iZKS0qVLl1RUVKysrPLy8lrbMzIyxowZIyIiQhCEg4NDbW1tSUkJeTEByIcTMADoMu0LoMOHD0+cOHHVqlX19fUODg4HDhwgJRafERMT271796+//jpw4MBz586xG3V1dVvvac/Pz29paVFUVCQvIwD5cAIGAF2mfQG0Y8eOu3fvmpmZUanUzMxMPAbPQf7+/gkJCSEhIYsXL25qapoxY8bTp0+HDRs2Z84cR0fH7du3Y/5oEHA4AQOALtO+AKqpqWmdwlhSUlJUVLTTQ7e0tGRkZCQnJ1+4cOHhw4ctLS2dj8kvLC0t09PTHz9+PGTIkPLy8nv37s2ZM8fMzOzKlStYSBUAJ2AA0GXaPwXm4+MzduzYt2/fHjx48PDhwz4+Pp0bl6emoucpCgoKCQkJGzdu/OWXX6Kiotzc3MhOBMArOHgCBgDQsfYF0OrVqxMTE2/evPnmzZsNGzbY29t3btwfnYpeoAgJCa1YscLJyWnSpEmjRo1avHgxJowGIDh3AgYA8E3tCyCCIJydnZ2dndnbVVVV0tLSnRj3R6eiF0BWVlbjxo3btWvXvn37evXqFRcX17dvX7JDAZCJUydgAADf9G8B9PDhw2XLlpWVlY0ePTo4OHjZsmWFhYVZWVnFxcWdGDcoKMjCwsLV1VVDQ4NCoTCZzPPnzwcFBXEuebd37ty569evv3v3LiEhISgoaMSIEUVFRbgPGgQcR07AAAC+6d8CaNq0aePHj3dycjpy5Ei/fv0CAgJ8fHx69+7duXEDAwNdXFySk5OZTCZBEOyp6DU1Nb/Wf+7cuc+ePWvbkpub27lddxfXr1+fMmWKvLy8r6+vlZVV3759R4wYceTIkXZXzgAEAWdPwAAAvunfAujjx4/h4eEEQRgbG8fHx7O3f0bbqej37dvXQfVDEMTcuXPZpVKroKAgOTm5n8zAy5SVlZ8/f87elpOTk5SUtLa2Njc337Fjh7e3N7nZALoYZ0/AAFpVVVVRqVRcXIfP/VsACQsLszcUFBSkpKR+ctykpKS28xqvWbOmR48eBEH4+fl9sb+JiYmJiUnbFhkZGfYsyfzK39/f0tKypaVFT0/v4MGDv/76a1hYmKurq5+f35kzZyIjI1VUVMjOCNBFOH4CBpCRkeHn5/fq1SthYeH169cHBgaSnQh4S/t5gDjlyZMnAQEBKSkpmZmZmZmZdXV17A0u7a476tmzZ3p6uqamJoPB+P3338PCwgiCsLCwSE9PNzExMTMzO3r0KNkZAboIZ0/A4HucPXvWwsJCU1Nz6tSpHz58IDsOh7W0tHh4eKxcubKysjIzM3Pbtm137twhOxTwln8LoNevX0v/T9vtzo27aNGilJSUgoIC9qKGqqqqO3bs2LFjB4di8wlFRcWlS5du37599OjRrY09evRYu3bt33//vXPnzuHDh7948YLEhADAl+7evRsaGsouC2RlZflvItaCggIhISFPT0+CIHr16uXt7Z2SkkJ2KOAt/xZAtbW1H/6n7Xanh7azs0tISDh69GhoaGh9fT0n0goQOp1+7969kSNH2tjYbNq0qbGxkexEAFzE2RMwgiDS09Pj4+MrKipaWxISEjiR9KsKCwvnz5/v7u6+fft23j/inTt3LjAw0NHRUUNDY/v27bdv366qqiI7FCcpKCiUlpZ++vSJ/bKwsFBJSYncSMBr/i2AenzFz4xOo9GOHDliZWXVdkZE+E4iIiILFixIS0u7ffu2mZnZ1atXyU4EwC2cPQHbuHGjm5tbTEwMnU7PyMhgN86cOZNzedt7//69nZ2dvLz8pEmTrl+/zvsXVMTExOrq6tjbjY2NLBaLz+65lJWVnTBhgqOjo7+///jx42/evMm+GgTQqis+8V5eXl5eXl2wI77Uq1ev8+fPnz9/3t/f38bGZvPmzRoaGmSHAuCwnzzXaicyMjIrK0tBQSE7O9vd3T09PZ3b8wmdPXt2xIgRv/32G0EQLi4uGhoaJSUlvHzJwcvLa+jQob169dLT0/vPf/7j6uoqISFBdigOo9PpMTExFRUVNTU1BgYGVCqV7ETAW7h1EzRw1tixY588eaKrq0un0zds2NB66gYAn5OWlmbfSW1qajpnzpyQkBBu77GiokJRUZG9LSYmJiMjU1lZye2d/gwTE5Nz585duHBh2bJlpqame/fuJTsRh9XV1S1fvjw7OzsvL4/BYAgJCcXGxpIdCngLX13z5G+SkpJr16719/dfuHChiYnJpk2bJkyYgMktAD7n7e1tZ2c3a9asgICAkJAQNzc3Ly+v1ttBPvf06dN285BVVFTIysp+/x5HjBjh4uLi6elpamp64MABMTExHR2dzr+BLmFlZcXHNcHr169VVFTY18spFIq9vT3fT64LPwoFUDejra199uzZa9euLVy4cOfOnVu3brWxsSE7FABvWbVqlZ2dXWlpKUEQFArlzJkzMTExHUys+ueff+bl5bVtefv27Q8VQHQ6fevWrRMmTGAwGPb29mfPnsXJCbm0tbVLSkpycnKMjY0bGxv/+eefxYsXkx0KeAuFZxcotbKyUlJS4vaDG91XS0vLsWPHwsLCLC0tN2zYYGhoSHYi4HNFRUU2NjYMBoPsIF3he44/LBbr9OnTd+/e7d279/Tp0zF9Ea+JiYkJDg62sLDIy8uztbU9evQoqtLuixvHH9wD1F0JCQn5+fnl5eVZWVk5ODjMmjWrqKiI7FAAvMjR0ZEbwwYEBOzcuVNTU/P+/fs2Nja1tbXc2At0mqen56NHj4KCguLj448dO4bqB9pBAdS9SUhILF68OC8vT0FBgU6nL1iw4P3792SHAuAtzc3NHB+zuLj4/Pnzly9fDg0NPX78eK9eveLj4zm+F/hJysrKo0eP7tu3L9lBgBehAOIHcnJyGzZsePLkSXNzs4mJyaJFi1AGAbTixqw8b9++VVdXb310XF9fHwvXA3QvKID4h7Ky8s6dOx89etTQ0GBiYrJgwYJ2D7YACCZuTIFoYmLy7t079uoKRUVFcXFxDg4OHN8LAHAPCiB+o6amxi6DKBQKexIUrCYGwHE9evQ4ceLEtGnTNDU1TU1NFy5caGFhQXYoAPgBKID4k5qa2rZt2/Ly8hQVFW1sbLy8vB4+fEh2KAC+MmjQoIKCgjt37rx79y4oKIjsOADwY1AA8TNFRcW1a9e+fPnSyspq3LhxQ4YMSUhIaGlpITsXAJ8QEhLS0NAQExMjOwgA/DAUQPyPSqWGhoa+ePFi5syZq1evNjY2joiI4LOVnwEAAH4ICiBBISoq6u3tnZaWFhUVlZqa2rt37/nz52NueACA7uXp06dJSUl4xuXncbEAamlpycjISE5OvnDhwsOHD/HNC4+wt7ePiYnJysqi0WhDhgxxcnI6c+ZMY2Mj2bkAAOAbZs2aNXr06MjISHNz80OHDpEdp3vj1lpgV69enTNnjqqqqrq6OkEQTCaTwWBERUUNHDiQS3uEH6KhobF27dqwsLCzZ89GRkYGBwdPnTp1+vTpenp6ZEcDAIAvuHTp0v379/Py8nr06FFQUPDLL794eHhIS0uTnau74lYBNG/evMTERF1d3dYWBoPh4eFx7949Lu0ROkFMTMzLy8vLy+vZs2f79+93cHAwMDDw9/d3d3enUqlkpwMAgH89evRo6NChPXr0IAhCW1tbW1v72bNnmH+h07j1FVhzc7OqqmrbFmVlZZ5deBUMDAy2bNlSWFi4YMGCs2fPamlp+fn5XblyBV9cApDi7t277u7u9vb2v/3226dPn8iOAzzB0NDw9u3b7MPyu3fvCgoK2l5lgB/FrStAQUFBFhYWrq6uGhoaFAqFyWSeP38eU2XwOFFRUVdXV1dX15KSkhMnTixZsuTt27fe3t4+Pj7m5uZkpwMQFE+fPh0/fvymTZv09PR27do1Y8aMEydOkB0KyDd69OiDBw/279/f3Nz88uXLK1eulJWVJTtUN0bh3lUZBoORnJzMvlNdVVV19OjRmpqaX+vs7Ox8586dti1VVVVmZmYPHjzgUjz4Hjk5OSdOnDhx4oSoqOjEiRM9PT379OlDdiggR1FRkY2NDYPBIDtIV7CyslJSUkpISCBl7+Hh4SwWa/Xq1QRBNDQ0KCkpvX37VlJSkpQwwGuuX79eUFAwYMAAExMTsrN0HW4cf7h1BYggiPfv3ysrK0+cOFFGRobdkpCQMGbMmC92Pn36dH19fdsWJycnZWVl7sWD72FsbLx27dq1a9empaWdPHly1KhRVCrVw8PD3d3dzMyM7HQA/KmhoaH1zlYREREhIaGmpiZyIwHvGDRo0KBBg8hOwQ+4dQ/Qxo0b3dzcYmJi6HR6RkYGu7GDJQklJSXl/j8REREKhcKlePCjLC0tt23b9vr164MHD9bW1rq5uenq6i5atOjmzZvNzc1kpwPgK+PGjdu9e3dKSgqDwViwYIG1tTWNRiM7FAC/4VYBFBkZ+fDhw+jo6PPnz0+cOBHzDvMHCoVibW29ZcuWFy9enD17VkpKat68eaqqqv7+/rGxsZWVlWQHBPiCffv2kR3hxwwYMCAiImLJkiV2dnaVlZXHjh0jOxEAH+LWV2DS0tJSUlIEQbAXJA8JCYmKiuLSvoAUZmZmZmZm4eHhhYWFFy5c2L9/v7+/v6Wl5ahRo0aOHIlbhYBESUlJJSUlrS/XrFnDfnLYz8+PvFA/hv04AtkpAPgZt64AeXt729nZsU+8QkJCysrKvLy88DAnX9LS0po7d+7FixeZTOb8+fNfvHjh4uKipaU1Y8aMmJiY0tLS4uJizDQNXenJkycBAQEpKSmZmZmZmZl1dXXsDbJzAQAP4dYVoFWrVtnZ2ZWWlhIEQaFQzpw5ExMTIycnx6XdAS+QkpJycXFxcXEhCOLZs2f//PPPH3/84eXlRaFQhISEJk2aFBkZifkVoQssWrTI1tZ22bJlvr6+Q4cOvXz58o4dO8gOBQC8hYuPwf8kch9DhZ/X0tIiKyurpqY2evTolJSUR48eiYuLm5ubDxo0aODAgXZ2dpjBvXvpdo/BV1ZWBgcHy8vLJyYmPnv2rIOefn5+T548aduSm5trbGyclpbG5YwA8F262WPwIOCysrJqamrS09Nb7wbT1tYODQ29fv365s2bHzx4oK+vb29vb2dnZ2dnx14zDoCDaDTakSNHTp48mZub23HP1atXsy9Xt5o6daqSkhI30wEAyVAAAXe1vcQoLCw8ePDgwYMHEwTR0NDw4MGDW7duRUdHBwUFSUhI2NraWltbW1lZmZubs29ZBfh57NXuOu7Tu3fv3r17t22RlJQUEuLWLZIAwAtQAAG3mJmZUalUCwuLoKCgO3fu5OXlbdu2rfWnYmJitra2tra2ixYtIggiPz//zp079+7dO3r0aG5ubp8+fQYMGGBpaWlhYWFsbCwsLEze+4Buz9HR8dq1a2SnAADeggIIuEVISOjy5cvjxo1bsmQJQRBr164dPnz41zrr6+vr6+v7+voSBFFbW5uRkZGWlnbp0qVNmzYVFRWZmpqa/0/fvn1xfQh+CObqBIDPoQACLrK0tHzz5s379+/l5eVFRL73w8b+OszW1pb9sqqqKiMjIyMjIzU1NSIiIj8/X09Pr9//9O3bV0NDg2vvAPgBu7DmtitXrvz2228MBsPOzm7Lli24rQ2Ax6EAAq7r2bPnz/y6tLT0wIEDBw4cyH7Z0NDw+PHjR48eZWdnb9u27fHjx3V1dX379u3Tp0+fPn1MTEyMjY3V1NQ4ERz4RAeL8HBKXl7e5MmTDxw40Ldv36ioKDc3t3v37nF7pwDwM1AAQTcjJibWv3///v37t7aUlpZmZ2fn5OQ8fvw4Li7u6dOndXV1hoaGxsbGhoaGBgYG7O/XJCQkOBIgIyNj9erVr169srKyWrNmDZbsBYIgEhMTfXx8nJ2dCYJYs2bN0aNHX79+3atXL7JzAcBXoQCCbk9BQcHR0dHR0bG15ePHj7m5ubm5uXl5ecePH3/27NnLly979uypp6enp6enq6vL/l8dHZ0fnYvozZs3zs7O69ev79+///Hjx11dXe/cucMHq/ZWVVUlJSU1NDQMHz4cJV0niIiI1NfXs7dbWloaGxtFRUXJjQQAHUMBBHxITk7OxsbGxsamtaW5ubmwsPD58+fPnz9/8eLFnTt3nj9//vLlSykpKR0dHW1tbfaD0L169WL/79cuFyUlJTk7O0+bNo0gCDMzM319/efPn+vr63fN++KSwsJCe3t7CwsLSUnJRYsWxcfHW1lZkR2qmxk3bpyVlZWpqWm/fv2ioqL69u2L72EBeBwKIBAIwsLC2tra2traw4YNa9teXFxcUFBQUFDw6tWr9PT02NjY169fFxYW0mg0TU1NTU1NLS0tTU1NDQ0NDQ0NTU3NhoaGts/kCwkJtbS0dPm74bB169YFBgYuXbqUIIjY2NglS5bgofEfpaWldfHixQ0bNhw6dMje3n7Lli1kJwKAb0ABBAJNRUVFRUWl7bUituLiYsb/FBYWPnjwoKioiMFgvH37trm5+eLFi9ra2hUVFVVVVSkpKXl5ecrKympqaj179uyOj+jn5+f7+Piwt21tbYODg8nN002ZmZmdOnWK7BQA8L1QAAF8AbswsrS0bNfOYrEuXry4bt26/Px8LS2tSZMmZWVlXbx48d27d+wH/qlUqrKycs+ePVVUVHr27KmkpNSzZ09lZWXFNkh5Rx2g0+nnz58fNGgQQRDx8fHm5uZkJwIA4DoUQAA/gEKhVFdXM5nM9+/fa2hoTJkyxczMrG2H0tLSd+/elZSUMJnMkpKSkpKSBw8elJSUfPjwgf2ysrJSUVFR4Uvk5eXl5eXl5OTk5eVlZWWpVGrXvKmwsLBhw4bR6XQqlcpkMi9evNg1+wUAIBEKIIAf8PTp09DQ0PPnz5uZmf3111/jx49/9uxZ2zke2aVMByM0NTV9+PChtLT0w4cPZWVl7I2SkpK8vLyy//n48WN5eXlDQ4OsrKycnJxsG3JycjQaTUZGRkZGhkaj0Wg0aWlpOTk5aWlpaWlpcXHxTrwpeXn5+/fvP3jwoL6+fsCAAZ0bBACge0EBBPADrly54u7uzp6FaPLkyVu2bMnJyTE1Nf3+EURERNjfr32zZ0NDA7sS+tybN28qKioqKysrKyurqqrKy8vZGy0tLVQqVU5OjkqlSklJSUlJycnJSUlJSUpKSktL02g0CQkJKSkpGo0mLi5OpVKlpaV79OhBo9EkJSUNDAxkZWX54JF+AIDvwcUCqKWlJSsrq7i4uKmpSV1dnU6nY3Vl6O5kZGTev3/P3m5qaiotLZWVleXSvsTExJSVlX9oVp6Ghobq6ury8vLq6urq6uqampry8vKampqamprq6uqKiooPHz7U1NRUVFTU1dXV1NRUVlbW19dXVVV9+vSpvr4+LCwsNDSUS28HAICncKsAunr16pw5c1RVVdkL4jCZTAaDERUV1bqgAUB3NHbs2PDw8CVLllhYWJw8efKXX37R1NQkO9S/xMTE2DcSkR0EAIDXcasAmjdvXmJioq6ubmsLg8Hw8PD42vo4dXV1tbW1bVuampq4lA2g02RlZW/fvr1jx47Y2FgHB4fAwECyEwEAQGdwqwBqbm5WVVVt26KsrMxisb7Wf9y4cffv32/bUlVVhZsxgQepqKj8/vvvZKeAb6ioqKDRaBQKJT09PScnx8LCwtjYmOxQAMBDuFUABQUFWVhYuLq6amhoUCgUJpN5/vz5oKCgr/X//MnbxYsXKykpcSkeAPCxffv2bd++PTs7e8uWLQcOHHBwcFixYsXKlSu7YFl4AOguuFUABQYGuri4JCcnM5lMgiC0tLSSkpJ46m4JAOBX69evz87OFhUV3b17d1ZWlry8fElJia2tLQogAGjFxafANDU1AwICuDc+AMAXKSkpNTQ0EAQhKyvLfvi07QpuAAAE5gECAP6zfv16W1vbsWPHGhgYODo6Dh8+PDExsYM1zi5evFhYWNi2paSkRFpamvtJAYA0vFsACQsL79q1q4PFBaurq1++fCkqKtqVqTrQ0NAgJiZGdor/4rUwoqKiPDLDXmNjo7CwMI9MSdXU1EShUHjk4kRzczOVSu3du/fXOjQ0NPBI1G8aMWLE/fv34+PjZWVljYyMlJSUYmNjjYyMvtY/MzOzoKCgbYuQkNDDhw9/+eWXr/0Kjj8d4LUwOP58EY4/lA6ezCJXeXn5ixcvOujw8OHDrVu3hoWFdVmkDpSXl4eFhUVERJAd5L+mT5++e/duHjkGLVu2bO7cuRoaGmQHIQiC2LFjx6BBg9hTOZMuOjpaQUFh5MiRZAchCIK4dOlSZWXlypUrO+gjJyeno6PTZZFIhOPPz8Dx52tw/Pkaco4/rG4rNTXV3t6e7BT/xWQyVVVVyU7xL0lJyZqaGrJT/Fffvn2zs7PJTvFfrq6u586dIzvFf4WGhm7fvp3sFP8VERERFBREdgquGDRoEMfHxPGnAzj+fA2OP19DyvGHJy7EAQBwT3NzM9kRAIDnoAACAD7n6+tLdgQA4DkogACAz2H6HwD4HAogAAAAEDjC4eHhZGfoJCqVqqCg0LdvX7KDEARBiIuLy8nJWVhYkB3kv+Tl5W1sbHjkyU8qlWpjY8Mjj4RISEj88ssvNBqN7CAEQRCSkpImJibKyspkByEIgqBSqb169ergMVRoC8efDuD48zU4/nwNKccf3n0MHgAAAIBL8BUYAAAACBwUQAAAACBwUAABAACAwEEBBAAAAAIHBRAAAAAIHBRAAAAAIHBQAAEAAIDA6a4F0J9//mliYtKvX7+rV6+SlWHt2rUmJiaampobN27knVRbtmxJSkoiPc/ff/9taWnZq1evffv2kR5m6dKl2traurq6x48fJzFMTU3Nnj17Wl9+nqErU7ULw5sfZp7FC38c3vwnw/Hnczj+fDMMaR/mLl59niMKCwuNjIxqamry8/MNDQ2bm5u7PsPly5f79+9fW1tbUlKipaWVnp7OC6lu3rwpJiZ2+PBhFql/pY8fP+rr65eUlFRUVBgaGlZUVJAY5vr163Q6va6u7u3btzIyMtXV1aSEefTokb+/v5eXF/vl5xm6MlW7MLz5YeZZvPDH4c1/Mhx/PofjzzfDkPhh7pZXgJKTk11dXSUlJfX09FRUVLKysro+Q2lp6axZs8TFxRUVFe3t7RkMBumpKisrly1b5uXlxX5JYp74+PixY8cqKirSaLTs7GxpaWkSwwgLC4uJiYmKioqLi4uIiBAk/WX++uuv8vLy1pefZ+jKVO3C8OCHmZfxwh+HB//JcPz5Ihx/vhmGxA9ztyyA3r59q6Ghwd7W0NAoLi7u+gyenp4BAQEEQWRnZ9++fdvR0ZH0VMHBwb/99pu8vDz7JYl5GAzGq1evTE1NtbS0tm7dSqFQSAxjZ2dnYGCgpqamqakZHh4uJSVFSpgNGzYEBga2vvw8Q1emaheGBz/MvIwX/jg8+E+G488X4fjzzTAkfpi7ZQHU0tLSdpm9pqYmUmKwWKwdO3a4u7ufPXtWRkaG3FSnTp2Sl5cfOnRoawuJeWpra/Pz82/cuJGZmbl///709HQSw6SkpOTn51+7di0pKWnHjh1MJpMXPj+fZyA3FU99tmXlVgAAIABJREFUmHkcj/xxeOqfDMefr8Hx53uQ9WEW4dK4XKWmplZYWMjeZjKZampqXZ+hubl5woQJsrKyDx48YC/tS26qw4cPv3jx4vLly2/fvj19+nR9fT2JeZSUlEaMGCEnJ0cQhIODQ25uLolhkpKSJk+ebGRkZGRkZG1tfePGDV74/HyegcRUvPZh5nG88MfhtX8yHH++BsefbyLzw8yle4u46vXr13379q2vr2cwGHp6eqTc7nfixImJEyfyWioWixUSEsK+CZHEPDk5OX369CktLX3//r2WltbTp09JDLNv377Ro0e33mGXkZFBVpjLly+33vf3eYYuTtU2DM9+mHkTL/xxePafDMefdnD8+WYYEj/M3fIKkJaW1ty5c+3t7QmCOHDggJAQCV/kpaamJicnq6qqsl/u379/zJgxpKdqi8S/kpGRUUBAgJWVFYvFWr58ubGxMUEQZIXx9/fPyMgwMTEhCGLRokV0Op3EMK0+/9ch8d+L9z/MPAXHn++B4w8bjj/fROKHmcJisbg0NAAAAABvwokdAAAACBwUQAAAACBwUAABAACAwEEBBAAAAAIHBRAAAAAIHBRAAAAAIHBQAAEAAIDAQQEEAAAAAgcFEAAAAAgcFEAAAAAgcFAAAQAAgMBBAQQAAAACBwUQAAAACBwUQAAAACBwUAABAACAwEEBBAAAAAIHBRAAAAAIHBRAAAAAIHBQAPGzly9fjhs3zs/Pb+/evWRnAQAA4CEUFotFdgbglsWLF3t7e5ubmw8cOPDGjRtkxwEAAOAVImQHAC5avXq1uLh4UVGRoqIi2VkAAAB4CL4C42cSEhJnzpwJCgr6888/ubeX69evjx8/vit/8Xs0NTV9+PDhiz/6+PHjq1evuLRfAADoFlAA8bPExMS7d+/GxcWpqKiQneWnZGRkWFtbq6mpzZgxo92PDA0NRUREREVFRUVFz54929q+YsWK0NBQgiBev35tY2NjZGS0a9cu9o9WrlxZX1/fZeEBAIAHoQDiZ3FxcXl5eZ6enlOmTPme/o2NjV/c7rgnt9XV1bm5ue3fv7+goKCwsPDixYutP2KxWGVlZU1NTY2NjY2NjW5ubuz269evHzx4kL29cePG8PDwJ0+eHDlypLKy8tmzZwRBGBoadll+AADgQSiA+MeSJUssLCzq6uoIgoiIiAgICDhw4EBCQsLp06ePHTvWtmdERISOjo6+vv7ChQtbWlpu3brl4+Pj7e29Zs2attsEQaxfv15XV1dHR2fx4sXten4eoN2wI0eOPHfuHPtHlpaWN27caNfhO99XQkLCgAEDTE1Ne/ToceHCBUdHx9YfvXnzRlVVtV3/8vLyRYsWrVq1iv1SRkbm6dOnDAajrq5OVFR0/fr1K1as+M5dAwAAv0IBxD9+//13U1PT6Ojourq6nTt3/vbbb1/sdv369aioqHv37mVmZjKZzI0bNxIEkZiYGB4evnbt2rbbf//997FjxzIyMnJycu7cuRMdHd2uZ8fDenl5sb+TevnyZVlZWUtLy+f7bSsxMdG6jVmzZrHbCwoK6uvrLSwsVFRUpk2b1tzc3PorBQUF5eXlTk5OKioq/v7+tbW1BEEEBgaGhYX17NmT3WfFihV5eXkzZ878z3/+8+TJE3V1dTU1tZ/+YwMAQPeGp8D4yuzZsxcuXFhZWenh4aGurv7FPleuXKmoqPD09CQIoqSkpKqqytHRkU6nt34r1LqdkpLi4eFBo9EIgvD29k5JSdHR0Wnbs+Nho6Ojly5d2tjYeOrUKT8/v5SUlHYd7O3t247g7Ozs7Oz8+cg1NTX379+/efNmz549fX19N2/evHr1avaPZGRkQkJCgoODKysrfXx8Nm/erKenR6VSXVxcTp06xe5Do9H27NnD3p4wYcL+/fsvXLgQGRlpamq6bt26Hj16/OgfGQAA+AAKIL5ibW1dWloaERGRnp7+tT5UKnXWrFlLly4lCKKurq6hoSE7O1taWrq1Q+s2i8WiUCjsbWFhYfall7Y9Ox6WRqNZW1tfvXo1Jibm7Nmzp0+fbtchIyOj7QiJiYltLyyZmZmx529UUVEZPny4jo4OQRDu7u5nzpxp7WNiYtKnTx9hYWEFBQUfH5+4uLjMzMx79+5pa2vX1NR8+vTJ2dk5MTGR3fmff/4ZMGBAc3NzeHh4XFzcgQMHIiMjFyxY8L1/XAAA4CP4Cozf9O7d29bWVk5O7msdhg8ffuTIkZKSkvr6end391u3bn2t55AhQ2JjY6urqxsaGqKjo4cMGdLBfr84rJeX18aNG+Xk5LS1tb+5X2dn57tttM5ePWbMmMuXL798+bK6uvrkyZO2trYsFisvL6+pqSkyMtLZ2bmurq6ysvL06dP29vZxcXFMJrOgoCAiImL8+PGt1Q+LxYqMjAwODq6pqZGTk9PU1DQ1NS0tLf3OvyoAAPAZFEB85datW6WlpSkpKR085k2n00NCQmxsbHR1dXV1dUeNGvW1niNHjpw0aZK5ubmxsfGAAQN8fX072PUXh3Vxcbl//76fn98P7bcdDQ2N7du3jxgxwtjYWFNTc/78+Q0NDUZGRm/evJk7dy6dTu/Tp0/fvn319PTmzZv3tUFiYmLGjh0rLi6upaVlb2/fp0+fiIiIDvoDAAB/w1IYfMXOzm7Tpk07d+6cOHHihAkTyI4DAADAo3AFiH+cO3eORqM5ODhMnz699bZfAAAA+BwKID7R3Ny8fPnydevWEQQxcuRIJSUlOp1eUlJCdi4AAABehK/AAAAAQODgChAAAAAIHBRAAAAAIHBQAAEAAIDAQQEEAAAAAgcFEAAAAAgc3l0L7NOnT8XFxWSnAIB/UanUnj17kp2iK+D4A8BrOH784ZXH4AcPHnzt2rV2jeLi4mpqamTEAYD2mpqahIWFX758SXaQrrB8+fK9e/fKysqSHQQACII7xx9euQJ09erVdi1WVlZKSkoJCQmk5AHoAKu+oeZmOnWoDdlBulRRUZGNjaC85aampqVLly5atIjsIABAENw5/uAeIIAfVnHu0ofd0XXZz8gOAgAAnYQCCODHNH34WHUxVX6qe2lUDKu5hew4AADQGSiAAH7MxyNnpZ0H0cYMFlGQq750k+w4AADQGSiAgOuqLt1qqf5EdgrOqM99WZ//SsZlKEEQ8tPcy08nt1TVkB0KAAB+GAog4K66nBel+05+PMkXN7OzWGWHzsj5uVF6iBEEIaqhImlNLz+TTHYsAAD4YbzyFBjwJxbr4+FYhYCJ5acSpYfbi2l170kN6vNf1b8oLNkWVdKmUUi8h5zveIqwMGmxAADgx6EAAi6qunKbEBGRdrIjmlvKDp1R+W0e2Yl+Sg8D7d5n/iA7BQAAcAC+AgNuaamtK49Jkvf3ICgU6eH2LRXVn+4/IjsUAAAAQaAAAu4pP50s2b9vD10tgiAIISH5aR5lR86yGhvJzgUA/8fencdFVfV/AD+zDzMM+zYM4IILbohbiCi5lbuWUEGWlrmkoWFlpdmqaTs9j9nj405lmvvCUqbinqgIhAviggoMyDIbzL7c3x/0Q55RUJGZOwyf96s/Zg6Hez7AdP3ee889FwBwCQxshKLUx86aFaqagyfvNjIY+oIifq8u9MUCAAAgBAUQ2AqDEbxuOd0hAFoxi1ZnuFHM79GZ7iAAzgmXwAAAHJHit7Q7n/9oqpLTHQTAOdnwDJDFYsnLyysvLzeZTBKJJCIigslEvQUA8GDG0jvqY2dFIwfJf97tu2A63XEAnJCtCqDMzMw5c+aIxWKJREIIkUqlxcXF69evj4mJsdGI4FAog9GsUJkVKrOy1lKrNteoLTVqi0Zr0Wgprd6iN1jUWspgoIwmSq+njGZCiEWnI408WovB4zLYbEII04VHWCwGh83gcpkuPAaHzXThMwR8Jo/HFPAZLnym0IUpFLBcBUxXIdNVwBQJmXyeXX9ygJYg27TTffIo0dODS5OW6S5d5XfHhTCAFmarAmj+/PlpaWmhoaH1LcXFxXFxcVlZWTYaEWhhVtSY7lQay6tMldXmSrmpSm6SKcwyBWU0sTzcWB5uLDdXpkjIFAlZrkK2twdT6MLg85l8LlPowuBwGFwOg8dlcNiEEAaf19hygpTeQJlMhBCLVk/MZspgpIxGi0ZHGU0Wnd6i0VI6g0WrtShrTNIKs1pjUWssNZq6wotYLEyRK8tNyPJwY7mLmG6ubC8PpruI7enO8nRjebkzBS52/ZUBPIgm+4Kpolo0JobBYnlOmSTbsDPwq3cJzqADtChbFUBms1ksFjds8ff3pyjKRsOBPVCUqVJmuCU1FpcZS8qM0gpj6R0Gm80O8GEH+LD9vLmd2wmi+rC9PVhe7kyhoAVHZvC4dU+faMZmKaPRrKy1qGrNyhqzqsairDXJFJaiYpNcaZarzNUKymJhe3uwPN3ZPp4sb0+2lzvL14vt68X29mQKURuBvVEmszxlt9ersXXHA8LovjV/HK85/JdoZDTd0QCciq0KoMTExH79+k2aNCkoKIjBYEil0n379iUmJtpoOLAJijKWVeqv3TJcu2W4WWIoKmEKXTghgdwQMT88TDR2KEfs5/glAoPDYft4Eh/PxjpQeoOpSm6WKU0yublKYbgtNedcMlXKTJUywmCwfb3Yfl5sP2+2nw/bz5vt783x96mrxgBsQXMmz1hWUb1uW32LRaM17fkTBRBAy7JVATR37twJEyZkZGRIpVJCSEhISHp6enBwcGP9i4qKZDJZwxa1Wu3p2eg/WmAjlNGkv3pTf/m6ruC6vrCIKRTwOrfjhrbzGNCL2yGY6dqS53UcBIPH5Uj8ORL/e79kUWtMFTJTpcxUUWW6U6n7u8B4p8p0p5op5LP9fTkBvmyxDyfAlx3gyxH7tuxJL2izBJERQT98bNXIwFQ2gJZmw7vAAgICpk6dyufzlUrl2bNnjU0uAfzpp59euHChYcvNmzcFAvyLYhcWi/5GsTb3su5CoeH6bU5wAC8sVPTUYJ/El1nuIrrD0YkpFHA7CLgdgqzazTKlsbzSdKfKWFapOfO3sbzSVFbJYLPZYl+OxJ8j9uME+rED/Thiv7rpTUCvNWvWzJo1i+4UD4vBYrL9fehOAeD8bLV3/uWXX+bMmePm5vbWW2+tWbOmU6dO58+fX7Zs2WuvvXbf/ps2bbJqiYyM9PX1tVE8IIRYatSanEvac/na/CtsLw+X3mHuE0fwuoUyXfh0R3N0LC93lpc76d6pYaNZUWOU3jGVVRrL7tQeO2OUVpgqqlle7hxJAFfiz5b4cyUBnKAApzyL5mjS09MrKyvr33722Wc8Ho8QMm3aNPpCAYBjsVUB9OGHH+bl5bVr1y4yMvKdd96ZOXNmVVVVVFRUYwUQ2I2pSq45k6c5nWcoKuaHdxX07en1aizL053uXK0ey0PE8hA1rIoos8VUUW0sKTeW3tFfKao9dMpYeofBYXOCxBxJADc4gBMUwAkJZLm50hjbKV28eHHJkiXx8fFeXl6EEJ1Ol5ubSxovgMaNG/fXX381bKmtrY2Kilq4cKEd0gIALWxVALHZ7I4dOxJCXnnllZEjRxJCPDw8bDQWPAyzQqU+laM+cc5UVunSv6fbpBEu4V0ZHA7duZwZg8XkiH05Yl8yoFd9o1mmNJSUG0vKDbek6pPZhuIyBoPJaRfICQrghgRyQsTcYDGmEz2mhQsXDho0aNGiRVOnTh0xYsTBgweTk5Ob6L99+3a9Xt+wZeTIka6uKEwBnJmtCqCYmJhJkya9++67dXd+Xbp0acWKFYMGDbLRcNAYymjUnPm79kiWvvCmoH8vj+fG8MPDGCwsKEIblpe7i5e7S3jX+hazQmUsLjPcLjPcKK49ctpYXM4U8DnBYm47CSdEzA2RcIIC7DCXSLnnT2NZhc+cKbYeyD6io6NTU1PnzZuXmppqVdzcSyAQWM04ZLPZDAbDlgEBmklXcL3yu42S5A8c/yZcB2erveratWt3795tsfyzsO+NGzf69+//+uuv22g4uJfhtrT2z5O1J87xOoa4Do30WziTwcX5npZBmc26C1ddeoe1yNbqVozk97pbEpkqqo3FZYbbUm3uZdXeQ8bySrafN7ddILedhBMSyG0nYft6tcjQd0esVqj2HiIctr7gBi+sY8tunC5ubm4pKSlbt24tKCigOwtAC6Eo2YYdLFeBYnu61yuxdKdp3WxVADGZzNjYu3+b8ePH22ggsEKZzZrTuaqMo6YKmWhEVODX77MbXwIHmqcm45gsZVfgl+9yOza6ssPjYPt5s/28Xfr1rHtLmc3GknLjbanhlrTmj+OGW6WU3sBtF8RtH8hpJ+G2D+IGix+zupX/vEc0JoYj9qtevy3wq/eIE538iI+Pj4+PpzsFQMuoOfwXg8Px/zCxNGmZaMQgTrD4wd8DjcA9us7DUqupOXBc9fsxTqC/+4QRLv174VKXLZhrapW7/nCfMFy2cUfAZ0l2qBUYLBa3nYTbTiIccjeD8Wap4Wap/vL1moxjRukdtp83t30Qt72E2z6I2z6I5fEI6xfoC4v0Bdc9nhtDGYxMPq/2SJbrsIE2+UmcjkWnL1/8rU/iS9yOIXRnAedn0eoUv6X5vTebKRK6P/u0bOMO/4/m0R2qFUMB5AxMVXJV6uHaI1mCJ3r7fzCX205CdyJnpvg1VRgzwPPlZ8re/1p9KkcY3df+GVgiV1avrvVXzSiz2VhcbrhVYrhZqtz7p6GohMFmczsEcTsEcTsEc9sHcQJ8Gi3UKEq2YYfnS5Oq12w1Vcp8k16p+HqdILI3HpH2MJS7DlgMhuoNO8RLFzjTaTNwTMrtGYI+PXihIYQQ0Zgnaw6d0mZfqD9VDI8KBVDrZqqoVu78Q30mTzQ8SvLdBywvx7ub3WKpTN7o8fxY5zhVaywu05zJk3y/hDAYXtPjKpM3Cvr3pP3JGAwWi9tewm0vIU/+02KqkhuKSgw3S9Qnzsl/3mOpUXPbSbgdg7jtg7gdgjnBYgb7n+fO1h7JMpbe0WZfNErvMEVCxbZ0Boet3P2n55SJtP08rYTpTlXtnyfF37xf+fVa9Yls4ZD+dCcCZ2a6U6VKPyIcMkD+y966FpaHmyxlt6RPdzwot3lQALVWpiq5Ynu69ky+aPSQoJUfO+zyejWH/9LmXzHXagI+doZTtbJNOz1eGMcUCQkhvK4deV07KPce9Hh+LN25rLF9PNk+noL/v/3eotYYikoMRSW6C4Wq1ExjeSVH4s/tEMztEETMZtHwgTWZWS59ezCFLurMLOGwSOeoVm1N9tNutwnD2d4eXtPjKr9dL3ginPZSGJwYg8vxiP+f2bQuvcMYbDZOPTYbCqDWx6ysUe78o/b4WbdRMZJVHzvypQqLVqfYlu7/wdzq1Vs0WXmCyN50J3os+sIibV6BUVqh3HuwroXS67XZF91jR9U9uNthMYUCfs8u/J5d6t5SBqPhttRQVGy4UWwoKtHfLGFyOYRBhNH9mK4Cc5XcNWYAvYEdn+5CoaGoxDfpFUIIr0sHXrdOjlkKg9Ngebq7P/MU3SmcCgqg1oQyGFWpmar9h4VPDpD860PHX0FYsS1d0K8nr3N7r1diq1b/6tK3e6teepHXuX3Qqk+sGhkctoNXP/dicDm8Tu14ndoRQswyRelbX/i++bJZUcPgcjwmjyp9c6mu4Do/LJTumA6MomQbd3KCAlRpR+oaWCKhcs+frsOjcNMlQGuBAqjVUJ88L/95N69ze/GKd9gBreBZicbySvWRM4HJiwkh/F5duO0CVamZ7s8+TXeux8BgON9TKlW/H6P0+up12+tbLGqtan8mCqCmUJTgiXDKaLLUauoaGDyu29ihtGYCgEeDAqgVMNyWyjbssKg1PkmvtKJ/luQbd/K6dtBdvl73lte5vXLXAdehkXjumEPxiB0tGmG9RHvdJCdoFJPp8cI4ukMAwGNBAeTQKL1BsS299kiWx/NjRU9Ft6ap/hTFFAkpvUFz6nx9m0tEN7OiBgWQQ2HwuM53WgsA4IFQADkubc6l6rW/8buFBiZ/4PjTfawxGD6JL9MdAgAA4P5QADkii1oj27BTV3Dd+/UXGz41EwAAAFoECiCHo82+UPXfrcLI3pLvFmNZEQAAAFtwlAJo//79ZWVlDVsqKiqEwjYwE5OiLBotUygghFB6gyxllzb3sm/SNH73znQnA4AHqzlwQhgzgMnn0R0EAB6NoxRAV69evXLlSsMWjUZjNBrpymM3NX8cV+4+IPn3R8bS8srvU3id2wV+u4jpwqc7FwA8mDbvcvXa30xVMs8XH/TkEIq6s2K125gnXfp0t0s0AHgARymA3nrrLauW3Nxcd3cnv13IUqtRbM/gSPwrvl5ruFHs9drztDxZEwCagTJb5Cm7vV9PUPyyz3VYFEfs20Rn9clsw81S2YbtgclL6h/EBgA0aj23VTsjxW9pggHhDB5Xl3/F750ZqH4AWpGa348x3UWiEYPcJgyrfz7lfVEGo3zzPt8Fr3IkAaq0TLslBIAmoACijbGkvPbEOd3fBRyxn9vkUaqMo3QnAoCHZanVKHf/4fVqLCHEbeJI4+1Sbd7lxjordx/ghXXkdwv1eiVWteegWa60Y1IAuD8UQLSp/H4jMZjcnx/rNT3OY/LT+uu3dZeu0h0KwElkZ2fv3btXqbxbaqSmprbg9uVb9gsH9eOGBBJCGGyW50uTZBt3UmbLvT1N1Yqa34/XTRJiB/i4Dh8o39qSSQCgeRxlDlDbQlGV320w3CwVRPbWX7mhv3KDEMJ0Fch/2iP+YiHd4QBavRUrVqxevXrw4MFJSUm7du3q06cPIWTmzJlWt5o2G2U2q4+dtWh1qvQjd1sZDMP1W7wuHaw6y3/azZb4aXP/OT/EdBcpUw+7jYrhdgxukTAA0DwogOzNotNX/TvFLFN6Tpvc8NZZbscQlocbjcEAnMaqVavy8vK8vb3z8/NjY2Ozs7NFIlET/bVarU6na9hiMpkoimqsP4PFCvn5m4cMw5H4MwV8w43b9S2i4VGUyfSQ3w4ANoICyK5MVfKKFat5XTr4vv0ag4U7QQBsQiQS1a0i1qtXrzlz5iQlJa1fv76J/pMnT87KymrYUlNT01JhPJ4f21KbgmYwSu9wxH6EwaA7CDgczAGyH8PN0vIPvnUdFuk9Ox7VD4DtJCQkREdHr1mzhhCSlJQkk8ni4+M1Gk1j/TMyMmT/q2/fvmKx2I6RwSbMMoV04Zc1f56kOwg4IpwBshNd/pXK7zd5TY8TRvdr2G4sr9ScynGf/DRdwQCcz0cffRQdHV1dXU0IYTAYO3bs2LZtm6enJ925wN7kv+wV9O2h2JYmjO7HFLrQHQccCwoge1CfOCfbtNN34Qx+WKjVl2Trt+vyr/C7d+KFdaQlG4BTGjFiRP1rFouVkJCQkJBAYx6wP31hke7iVcm/P5Kl7FJsT/d6JZbuROBYcAnM5moOnJD9tMf/w8R7qx9t9gVTRbX36y9Wr99GGp9xCQAAj4aiZBt2eL40icHjer44UX38nLG4Ze4BBKeBM0C2pdz5R23mafGyBWw/b6svUWazLGW316uxLhHdag+dqs087To8ipaQALSbM2fOfdv/85//2DkJOIfaI1mEyRAO7k8IYboK3J99WrZxh/9H8+jOBQ7kPgVQbm5ubm7ulClTysvLg4OxUkXzybfs157ND1i24L73t9ekHWGLfeuejOg1/bk7y/8jGBjBFOAqNbRF48ePpzsCOBGKkm9JZbrw7iz94Z8Gk0V36aru0jV+9070RgPHYV0Abdq0acWKFVqtNi4ubsiQIUuWLJkxY0bzNm2xWPLy8srLy00mk0QiiYiIYDLb0BU3xdZUbfZF/0/ns0Su937VrKxR7PzDN+kV050qQghTwOd1DFbu+sPzpWfsnhSAfuPGjbNqsVgsS5Ysubcd4MEYDL93XrPoDA3b3MloLD4JDVkXQMnJyadPn546daqrq2tubm7//v2bVwBlZmbOmTNHLBZLJBJCiFQqLS4uXr9+fUxMTAukdniylF26i9cCPpnPdBXct4P+6k2mq6B63baGjRa94b6dAdqIr7766vvvv5fL5f7+/hUVFcOGDaM7EbRW9y7JDWDFugBSq9V8Pr/utUAg4HA4zdvu/Pnz09LSQkPvTvstLi6Oi4uzWm3MwWlzLyt3Hwj4ZP4jLaIl37xPd/FawMfzmrjrUtC/l6B/r5bICOA8fvjhh8uXLy9evHju3LkqlWrjxo10JwIAe9BmXzBVyUWjhthzUOtrUlOmTJk4cWJRUdGGDRtGjhw5ZcqU5m3XbDZbLSPm7+/fxNLyDogym2UbdpgqqmsO//Xw36XYlq7NvhDw4RtYcwLgUfF4PDc3tz59+mRlZUVFReXl5dGdCABsjtIbqtf+Jv91v6lSZs9xrc8Affrpp2lpaSdOnCgtLV2+fPngwYObt93ExMR+/fpNmjQpKCiIwWBIpdJ9+/YlJiY+dmD7qZuk7Bk//s7nPwqj+jzM9GTl3oPqU+cDPn2TKRLaISGAk4mJiZk+ffqiRYuee+65O3fuuLu7050IAGxOuftPXrdOnEA/+c97fN+abrdxrQug6urqgQMHDhw4sP6tt7f1/dsPY+7cuRMmTMjIyJBKpYSQkJCQ9PT0Ju4pW7JkSWFhYcOWa9euseh7XoRZWaPc82fAsgWcQH9Bv57KHb97Tn226W+pOXSq5sAJ8dIFLPemHrsIAI1Zt25dYWFhly5dvvvuu+zsbFwCA3B6pip5zR/HxV+/x3JzLU1aprt0ld+9s32Gti6ARo8eXfdCqVRevXq1d+/eubm5zdt0RUWFv7//Cy+8UH8Yl5qa2tjNrqNHj+7du3fDlpycHIHg/jOI7UCxJVU4NJIT6E8I8XhxgnTB564jozmBfo3115z9W7E1LeDT+SwvHLMCNJNMJvPx8amuro6IiIiIiKA7DgDYnDzySz8SAAAgAElEQVRll2jck2wfT0KI55RJsg07A796l9jlnnHrAujs2bP1r/Pz86dOndq87a5YsWL16tWDBw9OSkratWtXnz59CCEzZ84sK7v/Wpz3Xmv75ptv6qdj25nhZmlN5l+eCRPrH6HHaRco/2m33/uz79tfd6GwevUW/w/fqCuYAKB5WvAADAAcn+7SNW1egWjcUMON24QQdoAPZTDUHP5LNDLaDqM3tRJ0r169jEZj87a7atWqvLw8b2/v/Pz82NjY7Oxskaj1XBiiKNHwKNOdyvoGToDvfRczJIQYi8sqv9/o+/Z0bvsge+UDcE4tdQAGAK2C4VYpO8BXtmFnfQuDzzdJK+wzunUBtGDBgvrXBQUFXbt2bd52RSKRUCgkhPTq1WvOnDlJSUnr169vdko743YI8p79UM9NNMuVd1as9nz5WbtdswRoIx7nAAwAWgW3MU+6jXmSrtGtC6CGl6LGjRvX7IXIEhISoqOjZ8+ePWvWrKSkpMmTJ8fHx2s0muYndTwWnf7O8v+Inh7s+uQTdGcBcAYtdQBGG4oyVcnZvl505wCAB7tbANUtURgU9D/Xcc6dOxcZGdmM7X700UfR0dHV1dWEEAaDsWPHjm3btnl6ej5eWkdCUVXfb+J1auf+zFN0RwFwEi11AEaX2szTsg07JP/+kOXlQXcWAHiAuwXQ999/TwiRSqXHjx8PDw93dXXNysp69tlnt23b1vi3N2XEiBH1r1ksVkJCQkLCQ11XahXkm/da1Brfd5r5oDQAaKhlD8BoYdHq5FtS+b3D5Jv3+czD7CUAR3e3ANqyZQsh5Kmnnrp8+XLdmeeioqJZs2bRFs2B1R49o/4rV/zFOww2bSsVATiTFj8Asz/l9gyXvt29pz9X+uZSXcF1fljog78HAAghhBhulpiVtS69w+w5qPWt9rdu3erSpUvd6/bt29ctYwgN6QuL5D/v9l/0+n0f8w4AzbBly5YtW7ZwudzLly/n5uaeOHGisLBQLpfTnethmcqrajOzPOPHM3hcj4QJsg07yWM/+cdSo9ac/btF4gE4NIul6odfqv6dYlFr7TmsdQEUFRUVFxe3a9eunTt3PvfccwMGDLBnGsdnVtRUfrfBZ84UTlAA3VkAnE0LHoBZLJacnJyMjIz9+/efP3/eYrG0UMb7k23a6fbMSJanOyHENWYAk8upPXrmMbcp37yv8tv1RnvdEgxAl5qDp5h8nmBAuGJHhj3Htb4LbM2aNatXr/7111+ZTOaTTz6JS2ANUWZL5XcbXEcOcunXk+4sAE6o7gBsypQpFEVt2bKl2QdgmZmZc+bMEYvFEomEECKVSouLi9evXx8TE9Oief+hv3JDcy7fotHq8i7XtZhVtfJf97k++QRhMJq3TcPNEs25fPdnnpal7PJf9HrLhbVGmc0M+h46BGDR6hTbM/wXzWb5eEnfXCYaHsUJFj/421rC3QJo3rx5S5cu/fDDDwkhdQ9yLywsfOedd1auXGmfKI5PnrKTKeB7xI6mOwiAc2qpA7D58+enpaWFht6dhVNcXBwXF1c31fpeRUVFMtn/PIZapVKxWKzs7Oy6t2FhYXULmykUiuvXr9d3q2vnhATyF0y7UVxc395l1BCRuxthMO7bv7HtNGyvXr3FpV+PPmOHqJf8S5tzSd8hsHnbabo9WKE17j8S+M0iZW2NLbaPdrQ/sF21/3CnXl24HUMIIWTUoCPLv/ea8fy9/fPz800mE2lRdwugkSNH8vn8+qXowYr6ZLY255L4y3ebfUgHAI1p2QMws9lct5F6/v7+VOOTcj799NMLFy40bCkqKiouLp49+5+n33z22Wdjx44lhKSlpSUnJ9d3q2tnuvAP3Sy0bo+MaKz/A9sttRqzQsXJCfisnf/Qac/KNu3M7N/++3/961G384B2iprv32VkWK+ajKNpqrLmbwftaG9uO2U0Gksrlv9nZd2T0o8alV/v+oV1NI0p4Fv1//LLL1UqFWlRDKudQkFBwcmTJ1977bUlS5YcOXLkiy++uPcpXfYRGRnp6+ubmppKy+j1LGotU+hiLK8s/+A7/yVvcDvgeRfQRpWUlERFRRU3OM/Rgvbu3Ttq1KhDhw5ZtY8bN64ZW/vxxx9Xrlw5adKkoKAgBoMhlUr37duXmJhYX9A8kC32P4Ybt9mB/kw+r+lulMFYmrTMJ/Gl+vXl73z+o0vvbm7jW3hVJMXWVFNFtfvzY8sXfRuYvLixp/0A2I5ia6pix+9WjcKoPr5vv2bVaIv9j/UcoDfffHP69Olnzpw5ceLEF1988fnnn2dk2HVSkkOpOXBCufP3wG8XVX673iN+PKofABuZNGkSISQ0NNTqAKx5W5s7d+6ECRMyMjLqplGHhISkp6cHBwe3ZOJHZFbWlH+y0nX4QK9XYpvuqTmXb6qUVXy5tr6FMpmMZRUtWwCZquQ1f5wQf/0e28dTOPQJxW9pD/n8H4AW5BE/3iN+PF2jWxdA5eXlL7zwwvvvv//SSy8NHjz4xo0btMRyBBa1RrEtjSP2u/PZDxxJgOgpezycFqAta8EDsODgYIe6h0Px635BZG/1sbOipwZzJP5N9BQO6nvvaigMDqdl88hTdonGDWX7eBJCPJ4fW/rmUtfrt3mhIS07CoAjs74NvmvXrm+//fbmzZvHjx+fnJzs5tZ2T4oqfksXRPYRxgww3Cxxn/w03XEAnF/dAdiuXbuc7ADMUFSiyb7g9Wqs++SnZesfvLQjUyiw+o/BbckCSFdwXX0616Kqlf+yV/7LXuXOP9heHvKf97TgEACOz7oAWrVqlUQi2bNnT0BAQHV1dd3y0G2QsaRcfeKc69PRit/ShCMGKbe33euAAHbjrAdgso07PBLGMwUuojFDTXKVNvvCg7/HllhuIs8pE1le7kxXQd1/goERrkOw6hu0LdaXwHx9fYcPH56bmxseHj579mx6r5rTSJayyz12lHzTTrfxw93GDi1NWqa7dLV+TiIA2MKqVat+/vlnJzsAU5/Mtmh1omEDCSEMFtPr1TjZmq2B4WEMjvXu1244gX54ijOA9f+BmzZtWrFihVarjYuLGzJkyJIlS2bMsMfzPmNjY3Nzcxu2lJaW9uxJz3qD2rwCbe5li7LGpFBxAnxlG3cwXQXylN3iL9+lJQ9AG+GUB2CKHb+bZcqSeZ/Vt5gqqtV/5bjG4IwLAJ2sC6Dk5OTTp09PnTrV1dU1Nze3f//+9imA1q5dq1AoGrZMnjzZ37+pqYK2ww0Wu08eVZN+xC12FEsoIIRwO4aw3EW0hAFoO+g6ALOpgI/nU3q9VSPLy52WMABQz7oAUqvVfD6/7rVAIOC09K0HjfHy8vLy8mrYwuPxGDQtOcgUCbVn//aeFS8c0p+WAABtE10HYDbF8hARgsMnAIdjPQl6ypQpEydOLCoq2rBhw8iRI6dMmUJLLHoptqayA3xQ/QDYGV0HYADQBlmfAfr000/T0tJOnDhRWlq6fPlyupaBppH+yo3ao2cCv11EdxCANqfuAKysrGzDhg2bNm1qmwdgAGAf97kNYdy4cfXLz9fU1IhEbejkLWUwVq36xXt2PGb8ANgfDsAAwG7uFkDnz59ftGiRTCYbO3bsvHnzFi1adPv27by8vPLychrz2ZnitzRuxxDBgHC6gwC0UW35AAwA7OnuHKBXX301KioqOTlZKpWGh4dLJJJFixadPn2axnB2ZigqqT16xuvVBzypBwBa3Pnz50eNGjVgwICPP/64qqpq5syZo0aN6twZK28BgK3cPQMkl8s/+eQTQki3bt327t1b97rtoMyWqh9/8XplMi5+Adjfq6+++uyzz44cOTIlJSU8PHzWrFlTpkxp37493bkAwGndLYBYLFbdC29vb6FQSFMe2ih3H2B7ewoH484vABq08QMwALA/69vg2yZjeaUqLdNrehzdQQDaqDZ+AAYA9nf3DNCtW7fq5xuq1er61zU1NTTksi/Z2t88Ykex/bzpDgIAAAD2cLcA0mq1NOawD/kve4SD+nE7/s8DhtTHz5kVNaIxT9KVCgDa8gEYANDibgHE4/FozGEHuguFqoxjukvXxJ+/Tf7/IRsWjVb+yx7ft19j/P8ZeACwvxY/AMvOzi4pKRk6dKi7+z9P3UpNTR0/fnzLjgIArZcN5wBZLJacnJyMjIz9+/efP3/eYrHYbqyHSSPbuNN33lRCkdpjZ+ub5b/ud+kfzuvSgcZoAMBrRPO2tmLFismTJ2/bti0iIiInJ6eucebMmS2XFwBavfusBN0iMjMz58yZIxaLJRIJIUQqlRYXF69fvz4mJsZGIzat5s+TTJFQMDCC5eVe+e16QWRvJp9nuC3VnM6RJC+hJRIA2MiqVavy8vK8vb3z8/NjY2Ozs7OxoCIAWLFVATR//vy0tLTQ0ND6luLi4ri4uKysLBuN2ASLWqPYnuH/4RuEEF6XDvwenVV7/vSIHy9bt90jfjxThFtOAJyKSCSqu5WsV69ec+bMSUpKWr9+Pd2hAMCx2KoAMpvNYrG4YYu/vz9FUTYarmmKbenc9hJiNhtu3CaECKL6Vn6/kSkSWnQ60fAoWiIBgO0kJCRER0fPnj171qxZSUlJkydPjo+P12g0jfU/fvy41TN/ZDIZThoBODdbFUCJiYn9+vWbNGlSUFAQg8GQSqX79u1LTEy00XBNM8tVZpW6avXW+ha22F+583e/918nTKyEBOBsPvroo+jo6OrqakIIg8HYsWPHtm3bPD09G+t/4MCBK1euNGyRyWR+fn42DwoA9GHY7qxMcXFxRkaGVColhIjF4rFjxwYHBzfWWSaTKRSKhi2TJ0+WSCRpaWm2yCbfst9crfBJfNkWGwdwSiUlJVFRUcXFxXQHsYfIyEhfX9/U1FS6gwAAIbbZ/9jqDBAhJDg4eNasWQ/ZeebMmbm5uQ1bSktL2WybxDNVK2oPnBR//Z4tNg4Ajmbo0KFHjhyhOwVAizFVyti+XnSnaPVsWAA9kp07d1q11B2B2WIs+U+7ReOeZPs0ej4cAJyJ2WymO4KjsljMChXLy4PuHPAI1Cezq374RfKvJXh6wWOyVQG0evXqo0eP3tu+ZcsWG434MCijyVBUrL9yw2fuFBpjAIA9TZ06le4IDkqZmqnae1Cy8iOmwIXuLPBQKL1B/vMel77d5T/v8X37NbrjtG62KoBefvnlrKwslUo1d+5cGw3xqCizWbpwBWUyeSRMYPC4dMcBADvBEoj3ZVbUqHb/ye0Uotjxu9fUZ+mOAw9FufcgLyzU540ppW8u0+YVuPQOoztRK2arAkgoFC5YsGDr1q0jRoyw0RCPqibjGKU3mmRKQf+edGcBAKCZ/Nd9riMGuk0YIV3wuWjkIE6gP92J4AFM1YqajGPir95lcDieLz8j27Qz8JtFDBbuZW4mG/7iwsPDly9fbrvtPxJLrUa5+wAhlEt4mHLnH3THAQCgk6GoRJtzyX3yKJa7yO2Zp+Qpu+lOBA8m/2m3aExM3fRnYVQflptr7eFTdIdqxRxlErStyTfv4waLCZvlk/iSdMHnrjjcAYA2TLZxh2f8+LqpP25jh0oPntLmXHLp053uXNAow80S9anzvK4d9AXX61rMyhrFllTX4VF4mHfztIkCyFhcps7KZTCZ/ovn1B/u+C16ne5cAAA00JzN112+znRz1eZe+qeJxZRv3osCyJFxxH7+H1ovJsxgs1D9NFubKIBkG3cwWCymUKA+dV596jxlMmvOX9T+fcUlvCvd0QAA7I0X1sH3rekNWwSEsNxd6coDD4PB4+LfrJbVJgogl8hw3ZUi1ycHMF0FdS2eUyayfbD0BQC0RSyRqzCqD90pAGjWJgogc4XcdehAz5eeoTsIAAAAOATnv33OrFDVHDrlPvlpuoMAAACAo3D+M0DKXQdchw5ke+OCF9Bmx44dX3zxBd0p7u/VV19944036E4BAGBvTl4AmWVK9fFzgckf0B0E2rQLFy7069fv4Z8NbDd79+7NycmhOwUAAA2cvABSbM9wHTmI5SGiOwi0dWKxuF+/fnSnsJabmyuVSulOAQBAA2eeA2SqkqtP57hPdJRncQAAAICDcOYCSLn7gOipaKZISHcQAAAAcCy2KoDMZnP965ycnM2bN1+5csUWA1Fmc8WX/zWWV1oHkCvVp867jRtmi0EBAACgVbNVARQQEFD3Ijk5efz48b///vuoUaN++umnFh+oJuOo7vJ1ecquf95TlFF6h9Sd/hkexXLH7B+Atm7NmjV0RwAAh2PzSdDffvvtuXPnxGJxZWVldHT01KlT79vt+PHj5eXlDVtkMplI9IDyxayqVe46EPBZUuU367S5l10iutUcOiXbuDNg2Vu1R89KvsfNXwBtUXp6emXl3bPCn332GY/HI4RMmzaNvlAA4FhsXgD5+PiIxeK6F6zGn9l24MABq2tkMpnMz8+v6Y0rtuwXPvkENyTQ8+VnZZt2BixNUmxNE/TtUfX9RtdhkSxP9xb5EQCgdbl48eKSJUvi4+O9vLwIITqdLjc3l6AAAoAGbHUJjMPh9OjRY8yYMSqVKiUlhaKoGTNmDB48uLH+S5cu3fa/OnXq5Onp2cQQhpslmrP57pOeKv/4XxyxL9vHs/LbDYIBvbxeiTWWVfC7dbbBjwVgK6Y7VZXJG4nF0ozvzcrKYjAYqamp9S1Lly5lMBh6vb7urVqtXr16dcsEbQ0WLlx4+PDhoqKi8ePHJycni8Xi5OTk5ORkunMBgAOx1RkgqVSq1Wpv3bp18+bNgIAAiqI6duz41ltvteAQsg07+d1CZSk7DTdLKr/fxOseWpt+zHvG8zWHTrl076zYkS54ohdhMFpwRADbkaXs0uUX1hz6S/RUdDO+3cfHZ8eOHePHj697m5qaWn8FOT8///vvv9doNK+//nqLxXV40dHRqamp8+bNS01NrS8EG7N58+aSkpKGLVKp1MXFxZYBAYBmNrwN3sXFJSwsbPTo0REREUwm84MPPmjBHQplNrN9PCijSfNXLq9ze1NFlfrYWU67QMpoqjlw3GvmC0w+rzbzdEsNB2BTuvxCwy2p/0eJit/SLGptM7bQt2/fnJwco9FICLlx44avr69Q+M8CEFu2bFEoFC0Zt5Vwc3NLSUmJjIwMDQ1tuqdWq5X/Lw8Pjw4dOtgnJ9iU/Je9pio53SnAEbXWlaAZLJbP/GlVP/zsNnG454sTlWmZ8o07mcHiyu82EEJk67eZVbXyLamuwwbiJBA4DkpvMCtU1q0WS/Xare7PjGS5ufJ7dpb9tMtj8ijrPkwm28eziQ8zk8kcNmzY4cOHR40atWvXrtjY2Ozs7LovLV++/NChQ+vWrWvJn6T1iI+Pj4+Pb7rPjBkzrFosFouvr6/NQoGdaLLyVKmHTZXVvgum050FHE5rLYAIIYYbt7V5BZJ/f0gIEY0YpDlxniUJMEorPJ4fyxb7EUKYPA6qH3AoqrTMmkN/WTVaNFpKb1DuOajcc5CyUJZqmTavgHHPHQO+86bywjo2sfG4uLiUlJRRo0bt3bt3z549ixcvbuH0rdbQoUOPHDlCdwqwN8pklv+yx/ft12Qbd+ouXeV3x8RQ+B+ttgCiqOp12/k9O2tzL9c1CPr3VOz8nde5vWjUEHqjATTGffIo9/89u2Op1ZTM+9TzxYmckMC6ltrDf1k0Wv/Fcx5144MGDZozZ87t27d5PJ63t3fLJHYKDddlhbZDte8QJyRQMCCcMphkG3YGfvUuYTrzww/gUbXWAoiyUByJP6XTa06d//8misFmC4dF0ZoL4NGYquS8DkHa8xe05y/UNzK4XGKxPOrOmslkxsTEvPnmm5MnT27pmK1bY8uPgRMzK2pU+w8HLH+bECKM7lvzx7GazNOiEYPozgUOpLUWQAwW0+eNlxq2aHMuGaUVoicH0BUJoBm47SX+H81rqa3FxcWNGDFi1apVLbVB5zBz5ky6I4C9yTfv5ffsTGm1hhu3CSHCmAGKLfuFUX2YAtzcB/9orQXQvVT7D7lNGolJP9AGRUZGZmRkEEKGDRtm+f+VhMrKyuo7jBgxYsSIEfSEA7A/ijLdqbbo9FWrt9a3sbw8jbfLmp5IB22KkxRAhttSQ8kdv+h+dAcBAAC6MRgBn71JdwhwdE4yI0y195Db2CcZ7EYftQEAAABQzxkKILNMqTmXj9ltAAAA8JCcoQBS/XHMNWYAUySkOwgAAAC0Dq2+AKKMptpDf4nGxNAdBAAAAFqNVl8AqU+c43YM5gT60x0EAAAAWo1WXwCpMo65jXmS7hQAAADQmtjwNniLxZKXl1deXm4ymSQSSd0z4Vt2CN2la5RW6xLRrWU3C9DiTp48+eWXX9KdwtrZs2c9PDzoTgEAQANbFUCZmZlz5swRi8USiYQQIpVKi4uL169fHxPTkpN1ajKOisYOxeKH4ODGjBmj0+nkcjndQax16tRp+PDhdKcAAKCBrQqg+fPnp6WlhYaG1rcUFxfHxcVlZWXdt/9333135cqVhi03b97k8XhNDGGWK7X5V7znTmmRwAC2ExkZGRkZSXcKAAC4y1YFkNlsFovFDVv8/f0pimqsf+fOnV1dXRu2XLhwoUePHk0MwXDh+709g+nCf8yoAAAA0NbYqgBKTEzs16/fpEmTgoKCGAyGVCrdt29fYmJiY/0nTJhg1XLt2jVfX98mhmDyefxeXVomLgAAALQltiqA5s6dO2HChIyMDKlUSggJCQlJT08PDg620XAAAAAAD8+Gd4EFBwfPmjXLdtsHAAAAaB6Hfhr8tWvXDh482NhXTSbTqVOnrGYO0Ugul3t6etKd4h+OFsbDw4PhGDfrqVQqoVDIYjnEc3M1Gg2LxWp6sr/dGAwGiUTSxGnaqqoqe+ahHfY/zeZoYbD/uS/sfxhNTEym19atW9evX99Eh8rKykuXLolEIrtFagJFUUql0nGWVHGoHZBSqXR1dXWQ/+dra2t5PB6Hw6E7CCGEaDQaJpPJ5zvERH69Xs/n8/v06dNEn44dO/73v/+1WyQaYf/zOLD/aQz2P42hZ/9DtVrHjx8fPHgw3Sn+IZVKxWIx3SnuEggEarWa7hT/6NmzZ35+Pt0p/jFp0qQ9e/bQneIfCxYs+O677+hO8Y+VK1cmJibSnaLVwP6nCdj/NAb7n8bQsv9p9Y/CAAAAAHhUKIAAAACgzUEBBAAAAG0OCiAAAABoc1pxAcTn8x1k+johhMvlCgQCulPcJRKJHOSuB+JgfylHC+Pi4kJ3in+4uLg4zm/G8TnUBwn7nyY41F/K0cK08f2P494G/zD0er2DrGFAEKZxCNMYo9HIZDId5J8Ki8ViNpsd5AbdVsGhPksI0xiEaQz2P627AAIAAABohlZ8CQwAAACgeVAAAQAAQJuDAggAAADaHBRAAAAA0OagAAIAAIA2BwUQAAAAtDmttQD68ccfu3fvHh4enpmZSVeGpUuXdu/ePTg4eMWKFY6T6uuvv05PT6c9zx9//DFgwIB27dqtWbOG9jDvv/9+hw4dQkNDf/nlFxrDqNXq1atX17+9N4M9U1mFccwPs8NyhF+OY/7JsP+5F/Y/DwxD24fZzk+fbxG3b98OCwtTq9VXr17t2rWr2Wy2f4aDBw/27dtXq9VWVlaGhIRkZ2c7QqoTJ05wudxNmzZRtP6W5HJ5586dKysrlUpl165dlUoljWGOHj0aERGh0+nKysrc3d1ra2tpCfP3339Pnz49Pj6+7u29GeyZyiqMY36YHZYj/HIc80+G/c+9sP95YBgaP8yt8gxQRkbGpEmTBAJBp06dAgIC8vLy7J+hurp69uzZfD7fx8dn8ODBxcXFtKdSqVSLFi2Kj4+ve0tjnr17906cONHHx8fNzS0/P18kEtEYhsVicblcDofD5/PZbDah6TezZcsWhUJR//beDPZMZRXGAT/MjswRfjkO+CfD/ue+sP95YBgaP8ytsgAqKysLCgqqex0UFFReXm7/DM8///ysWbMIIfn5+adOnRo6dCjtqebNm/fxxx97eXnVvaUxT3Fx8c2bN3v16hUSEvLNN98wGAwaw0RHR3fp0iUwMDA4OPiTTz4RCoW0hFm+fPncuXPr396bwZ6prMI44IfZkTnCL8cB/2TY/9wX9j8PDEPjh7lVFkAWi4XBYNS/NZlMtMSgKCo5OTk2NnbXrl3u7u70pvrtt9+8vLxGjBhR30JjHq1We/Xq1WPHjuXm5q5duzY7O5vGMIcPH7569eqRI0fS09OTk5OlUqkjfH7uzUBvKof6MDs4B/nlONSfDPufxmD/8zDo+jCzbbRdmwoMDLx9+3bda6lUGhgYaP8MZrP5ueee8/DwOHfunJubG+2pNm3adP369YMHD5aVlW3fvl2v19OYx9fXd9SoUZ6enoSQIUOGFBQU0BgmPT39pZdeCgsLCwsLGzhw4LFjxxzh83NvBhpTOdqH2cE5wi/H0f5k2P80BvufB6Lzw2yjuUU2devWrZ49e+r1+uLi4k6dOtEy3e/XX3994YUXHC0VRVFJSUl1kxBpzHP58uUePXpUV1dXVFSEhIRcunSJxjBr1qwZO3Zs/Qy7nJwcusIcPHiwft7fvRnsnKphGIf9MDsmR/jlOOyfDPsfK9j/PDAMjR/mVnkGKCQk5I033hg8eDAhZN26dUwmDRfyjh8/npGRIRaL696uXbt2/PjxtKdqiMbfUlhY2KxZsyIjIymKWrx4cbdu3QghdIWZPn16Tk5O9+7dCSELFy6MiIigMUy9e/86NP69HP/D7FCw/3kY2P/Uwf7ngWj8MDMoirLRpgEAAAAcEw7sAAAAoM1BAQQAAABtDgogAAAAaHNQAAEAAECbgwIIAAAA2hwUQAAAANDmoAACAACANgcFEAAAALQ5KIAAAACgzUEBBAAAAG0OCiAAAABoc1AAAQAAQJuDAggAAADaHBRAAAAA0OagAAIAAIA2BwUQAAAAtEyf67oAACAASURBVDkogAAAAKDNQQEEAAAAbQ4KIAAAAGhzUAA5sxs3bjzzzDPTpk3773//S3cWAAAAB8KgKIruDGAr7777bkJCQp8+fWJiYo4dO0Z3HAAAAEeBM0DO7NNPP42IiCgpKfHx8bHdKEePHn322Wft+Y0PQ6lUNvYluVx+8+ZNG40LAACtAgogZ+bi4rJjx47ExMQff/yR7iyP5b333gsODg4MDHznnXfqWvbs2RMWFhYaGvrvf//73v4qlSomJqbu9a1bt6KiosLCwup7LlmyRK/X2yc5AAA4JhRAziwtLe306dO7d+8OCAh4mP5Go/G+r5vuaWsHDhw4fPhwYWHhpUuXUlNT//jjD5lMNm/evAMHDuTk5Pz4448XL15s2H/58uWRkZEqlaru7YoVKz755JOLFy+mpKSoVKrCwkJCSNeuXe2WHwAAHBAKIOfx3nvv9evXT6fTEUJWrlw5a9as3bt3X7ly5fnnn3/55Zcb9ly5cmXHjh07d+789ttvWyyWkydPTpkyJSEh4bPPPmv4mhDy+eefh4aGduzY8d1337XqeW8Aq82OHj16z549dV8aMGDAsWPHrDo85M/l6uq6dOlSFxcXDw+Pnj17ajSatLS0oUOHhoSEuLm5JSQkbN++vWH/UaNGLV++vP6tu7v7pUuXiouLdTodh8P5/PPPP/jgg0f5vQIAgDOiwIlMmzZt3bp1Wq02NDS0pKTkvn2OHDnSu3fvioqK2tra+Pj4ZcuWnThxwt3dvaCggKKohq9///33rl27KpVKnU43ePDgn376qeFXG27wmWeeuXezGzdufPnllymKun79eseOHTMzM6061H1j/XZSU1MjG5g1a1bDUXbu3Dl27Njx48frdLq6IqaufcOGDTNmzLD6GUtKStq3b1/3WqlUzp49e+TIkQcOHDh79uyiRYse5zcMAADOgU13AQYt6fXXX3/77bdVKlVcXJxEIrlvn0OHDimVyueff54QUllZWVNTM3To0IiIiPqrQvWvDx8+HBcX5+bmRghJSEg4fPhwx44dG/ZserObN29+//33jUbjb7/9Nm3atMOHD1t1GDx4cMMtjBs3bty4cY39aFFRURaLZdmyZZmZmVZnj5q+Hufm5rZ69eq6188999zatWv379+/atWqXr16LVu2jMfjNfG9AADgrFAAOZWBAwdWV1evXLkyOzu7sT6urq6zZ89+//33CSE6nc5gMOTn54tEovoO9a8pimIwGHWvWSyW2Wxu+NUHbtbNzW3gwIGZmZnbtm3btWvX9u3brTrk5OQ03EJaWtrSpUvr3/bu3btu+aK6OUxRUVFxcXG3b9/eu3dvZGTkkSNH6rqVlJQEBwc/zC/nwIEDTzzxhNls/uSTT3bv3r1u3bpVq1a99dZbD/O9AADgZDAHyNm0b99+0KBBnp6ejXV4+umnU1JSKisr9Xp9bGzsyZMnG+s5fPjwnTt31tbWGgyGzZs3Dx8+vIlx77vZ+Pj4FStWeHp6dujQ4YHjjhs37nQD9Ys3lpaWfvDBByqVSqVSpaend+vWbdy4cQcPHqysrNTpdL/99tuzzz5LUdSVK1dMJlNj8SiKWrVq1bx589RqtaenZ3BwcK9evaqrq5v4iQAAwImhAHIqJ0+erK6uPnz4cBO3eUdERCQlJUVFRYWGhoaGho4ZM6axnqNHj37xxRf79OnTrVu3J554YurUqU0Mfd/NTpgw4cyZM9OmTXukca3MnTs3PDw8PDy8R48eYWFhc+fO9fX1TU5Ofuqpp/r37z9jxoy+ffsaDIawsLDS0tLGNrJt27aJEyfy+fyQkJDBgwf36NFj5cqV8+fPf8gMAADgZLAStFOJjo7+4osv/vWvf73wwgvPPfcc3XEAAAAcFM4AOY89e/a4ubkNGTLktddeq5/2CwAAAPdCAeQkzGbz4sWLly1bRggZPXq0r69vREREZWUl3bkAAAAcES6BAQAAQJuDM0AAAADQ5qAAAgAAgDbHVgVQ3aJ5dXJycjZv3nzlyhUbjQUAAADwSGxVANU/fjw5OXn8+PG///77qFGjfvrpJxsNBwAAAPDwbDUJ2tfXt+4WpKCgoLNnz4rF4srKyujo6MLCQlsMBwAAAPDwbD4HyMfHRywW171gsVi2Hg4AAADggWz1MFQOh9OjR4+QkBCVSpWSkjJ16tQZM2ZYPf27aQcPHty+fbuN4gFAM7Rr127x4sV0p7AH7H8AHE2L739suA6QVqu9devWzZs3AwICwsPDV6xY8dZbb7m4uNy3c2xsbG5ubsOWsrKyoKCgd955x0bxAOCRyOXyH374obi4mO4g9vDuu+/euHHj6aefpjsIABBim/2PoyyEKJfL5XJ5w5bJkydLJJK0tDS6IgG0KZdr/5YZqqK9hjfWoaSkJCoqqu0UQL6+vgsXLqQ7CAAQYpv9j60ugT0qT09PT0/Phi08Ho/BYNCVB6CtkeqKy/XSaNJoAQQA4ExsVQCtXr366NGj97Zv2bLFRiMCAAAAPCRbFUAvv/xyVlaWSqWaO3eujYYAgMdXUJt/R19GCLmmLlCaFEerDxBCRGy3vu4D6Y4GAGBDtiqAhELhggULtm7dOmLECBsNAQCPr1wvvaW9TgipNlZqzZqb2muEEB+uH925AABsy4ZzgMLDw8PDw223fQB4fEO9R9W9OFSVVq6XTpHMpDcPAIB94GGoAAAA0OY4yl1gAEAvN7aH1qyhOwUAgJ2gAAIAQggZ4BFNdwQAAPvBJTAAAABoc1AAAQAAQJuDAggAwBFJdcVrbyfTnQLAaaEAAgBwRLXmmmpDJd0pAJwWCiAAAABoc3AXGACAA7lQk5Ot/IsQojTJ7+jLUkp+JITwmPxY8UscBpfudADOAwUQAIADcWN7tHfpRAi5Yygr05bUvWYymCzCojsagFNBAQQA4EBCXDqEuHQghBSqL91QX3nS+2m6EwE4J8wBAgAAgDYHBRAAOLk1a9bQHaE5PDneXV170p0CwGnhEhgAOJv09PTKyrs3kH/22Wc8Ho8QMm3aNPpCPTJfrv+zAS/SnQLAaaEAAgBnc/HixSVLlsTHx3t5eRFCdDpdbm4uabwAmj59+t9//92w5fr165GRkQsXLrRDWgCgBQogAHA2CxcuHDRo0KJFi6ZOnTpixIiDBw8mJze1pPKSJUvkcnnDltZ1rggAmgEFEAA4oejo6NTU1Hnz5qWmpur1+qY7d+zY0apFKBSy2dg9AjgzTIIGAOfk5uaWkpISGRkZGhpKdxYAcDgogADAmcXHx2dkZNCdAgAcDgogAAAAaHNQAAEAAECbgwIIAAAA6GSmzAbLA25WaHEogAAAAIBOJ2SHdpVvtvOgKIAAAACATibKaKbMdh4UBRAAAAC0OVjpC6Ctq9CX3dAU3tBc7erao597FN1xAKCt2CbdVFCbTwhRmZQmynhdXUAI8ecFzm73th1GRwEE0OZozZobmsIbmsLrmsIbmkIek9dR0KWjoEs7F+sFkQEAbGek7/gnPIcQQs4qTlQZKsf4PUsIcWG62Gd0FEAAbUKFvuyapuCauuCauqDScKedS2gnYdcYr6deDXrDneNJdzoAaIu8OD5eHB9CyHV1gcFiaO9i10XbUQABOCcLZbmtK7pae6lQfemauoDJYHYWduskDIvxeirEpQOTwaI7IAAAnVAAATgPE2Uq0lwtVF8srL10TVPgxfHpLOze1z0yPnC6N9eX7nQAAPfX1bWnmB9s50FRAAG0bmbKfENTWFCbf0V98Ya6MIAv6SrsMcxn9EzBAle2iO50AAAPFsRvZ/9BUQABtD4WynJLe/1ybf7l2r+vq68E8CVhwp5P+0zo3K67C0tAdzoAgFbAhgWQxWLJy8srLy83mUwSiSQiIoLJxLJDAM1Xpi+5VJN3ufbvK7UXPThe3V17j/AZO6fdQgFLSHc0AIBWxlYFUGZm5pw5c8RisUQiIYRIpdLi4uL169fHxMTYaEQAp1RrqrlUm3exJvdiTS6Dweju2ru/e/TUoDlubA+6owEAtGK2KoDmz5+flpYWGnr3lrbi4uK4uLisrCwbjQjgNCyU5brmyoWanAs1OeX60q7Cnj1Evcf4PRvAk9AdDQDASdiqADKbzWKxuGGLv78/RVE2Gg7ACSiN8vya8/k15y/V5Hlz/XqJ+j4vfqWTMIyFW9adgokysRmYdgngKGz1f2NiYmK/fv0mTZoUFBTEYDCkUum+ffsSExNtNBxAK0UR6oam8G9V9t+q7CrDnR6iiHBRvymSmbjC9Ziys7NLSkqGDh3q7u5e15Kamjp+/HgaI/1486txfrGhwq40ZgCAerYqgObOnTthwoSMjAypVEoICQkJSU9PDw5u9C5/vV6v0WgatphMJhtlA6CdzqK9UJOTpzr3t+qcO9sz3K1/guS1ToIwJgM3CrSAFStWrF69evDgwUlJSbt27erTpw8hZObMmWVlZTSmMlB6vUVHYwAAaMiG52ODg4NnzZr1kJ0nTZp05syZhi01NTU2CAXQAjRm9Rbp+teC5z/qN8qMVbnKMzmqMzc0haGCrhFuAyb5x/tw/WwRsi1btWpVXl6et7d3fn5+bGxsdna2SNTUkkjffvttYWFhw5abN29yuVwbxwQAOjnKBenff//dqiUyMtLXF2vXgiOq0lecVZx8+AKoRHcrR5mVo8yqMlb2dus/1HtUYvv3eUy+TUO2ZSKRSCgUEkJ69eo1Z86cpKSk9evXN9G/a9euVhXSgQMHOBzO4ycxUcYNxSvNlJkQUqK7tf/OtqOyA4SQCLcBUZ5DH3/7ANBstiqAVq9effTo0Xvbt2zZYqMRAexGT+nM1AMu0VKEuqYuOK88naPMshCqn/vAeMn0ToJuuMhlBwkJCdHR0bNnz541a1ZSUtLkyZPj4+OtLrI3dO/coPXr1wsELbCkJJvBifIcqjNrCSFl+pIuwh5BLu0IIR0EnR5/4wDwOGxVAL388stZWVkqlWru3Lk2GgLAzsr0JTnKM4SQKkOFhbKkV+wihHAYnCHeI/lMl7o+FspcUHvhnPKvHFWWiOXez33gG+3fD3ZpT2PsNuijjz6Kjo6urq4mhDAYjB07dmzbts3Tk56H3vcS9a17cVR2oKtrj+6i3rTEAAArtiqAhELhggULtm7dOmLECBsNAWBnRZprOcrThBCdRUsIqXtNEUa4qD+Hx71c+/dZxckc1Rlfrn8/96hFnVb4cQNoTtyGNdzzsFishISEhIQEGvMAgKOx4Ryg8PDw8PBw220fwM46Cjor3AcSQqoMFWW6kj7uAynKUm2s3F+x7W/VOX9+YH/3QRP9X8Bz15thzpw5923/z3/+Y+ckttPXLdKfF0h3CgD4x30KoNzc3Nzc3ClTppSXlzdx4zpAWxPAk4z1m0wIuVJ78ZjsQKWh/LzytB9PPMA9erJ4ihfHh+6ArRi9K/TYx3CfsXRHAIC7rAugTZs2rVixQqvVxsXFDRkyZMmSJTNmzKAlGYCjoQh1Q12YpTh2Wn6cIpQ/L/Cjzt/ifE+LGDdunFWLxWJZsmTJve0AAC3C+oaU5OTk06dP9+7d29XVNTc394svvqAlFoBDKdXd3ln2y3uXZ28s+UHEdn+74yfRnsNH+z6D6qdlffXVV4GBgS4uLu3bt3d1dc3Ly6M7EQA4LeszQGq1ms//Z3kSgUDQIithALRScmN1luL4X/IjalPtE55D5rVfXH8/13TBPFqjOacffvjh8uXLixcvnjt3rkql2rhxI92JAMBpWRdAU6ZMmThxYllZ2YYNGzZt2jRlyhRaYgHQSGfRZiv++kt+5La2qJ9H1BTJzM7C7gzCoDuX8+PxeG5ubn369MnKypo+fXpSUhLdiQDAaVkXQJ9++mlaWtqJEydKS0uXL18+ePBgWmIB2J+FslyqzTslP/K36lxX157DfMb0duvPZuAkqP3ExMRMnz590aJFzz333J07d+qfY+ocUit2POEe7ccT0x0EAAi5twCqrq4eOHDgwIED6996e3vbPRWAXUl1xSflh/+SH/Xi+AzyHPpi4AxXdlOPjgIbWbduXWFhYZcuXb777rvs7GwnuwRWUJvf0aUzCiAAB2FdAI0ePbruhVKpvHr1au/evXNzc+2eCsAeNGZ1luL4CdkhhVEW5Tl0YehnYl4Q3aHaNJlM5uPjU11dHRERERERQXccAHBm1gXQ2bNn61/n5+dPnTrVvnkAbI4i1OXa/OPVf+bXnO8p6vNMQEIP1wg8oqvKcKfWVNOe1mdUOd8BGEWov+RHjRYDIURurP5blV1puEMI6SzsFsjHKmsAdGpqJehevXoZjUa7RQGwNZmx6oTs0AnZISFLNMRrxEtBs4UsV7pDOYo81blyvZTeAqg1HoBtlW54JiCh/mFwVsyU+abmqpEyEkJ0Fm25vlRP6QghIrY7CiD7uKYu6CDozGKw6A4CDse6AFqwYEH964KCgq5du9o3D0DLM1PmXNWZY9V/FmmuRnrGJLZfFOLSge5Q8ACt5QDsvPL0SJ/xfO79CyA2g/2iZGbd629ufPy070Q8DNXOUkp+fL3dOxJ+CN1BWpLWrHFhCehO0epZF0ANb/saN27csGHD7JsHoCVV6MuOyf48KcsU8yUxXk8ldnifw+DSHcqxGCmDwWIghBgseiNlVJtrCSFsBpvH5Ns/DA7AAB5IYZR99X/t3XlYVNfdB/Az+7DvyLCpgIBLKG6vC2g0mpBGq41i1OprqlGMRqJtlooaNTFq2tiSJ75aqzXVp31jYxLqwhIR3CBGDIhIcENl02Fghp1h1jv3/WN80aKgRWbOnZnv5687x8udLzOX4+/OnHvO7Q3bonfTDmLzHhRABQUFhJDg4H8bBFpYWDhmzBhrhwJ4NkbWWNxScLbhxF1t9Xjvyb+L2BqARSi78fe7fyluKSCEGFg9y5Ki5vOEkCBp6NqIbdYPYysXYHqTvqD5nIk1EUK0Jk1BU575tsFhbsN7mBzcU+jlLMRXrvCsWg0tjQYV7RT24EEB9NlnnxFC5HJ5Xl5eTEyMq6trQUHBq6++evjwYXrxAP4zSn3d2Ybs7xtPyaTBk3wSRniMFfJ6GugGS0KSSUgyISRXlaHQyRf8//c1VmZbF2BaU0dFRzlLWEKI0WS8q6uUGpwIIYHSkB4KoKWhmNfRSlqNze9dTTKwevPDD268bd54PXjl8z4v0cvVN3SslmEZ2inswYP/Gw4dOkQIefHFF69du2b+5LmioiIpKYlaNICnZmJNV9qKTquyKjW3xntN/l3ExwGSINqh4D9gWxdg7kLPRcErzNtlbZdnB/y3r9ifEKIzac83nR7vxdEPrhyHu9DzLzH3z5wPbrxtH2OASloLjyoOEUK0Jq2JZT66+Q4hRMAXruj/nrfIl3Y6m9T14riqqioyMtK8PWDAALlcbvVIAP+BVmPLucaTZxtOeIq8J/v8HKN8eo3PE1Bc7sM+LsAa9MrM+jQUQGAJEr4kyKk/IaTV2FKnk5u3CSHCHu/mhh50feHGjRuXmJi4YMEClmUPHTo0evRoKrEAnuiW+vqphszS1kujPMcnD1iHG7ue0VjPiXqTjm6GPrwAM5lMJSUlCoXCaDQGBQXFxsby+RaZ6slb5EtlwDg8JbtZxc9L5GOeqVXEE/MIMW8LeUKJAKdfL3UtgPbu3btnz54vv/ySz+c///zzNncFBnbPPP40V5WpN2kn+/x8YdByZ4EL7VD2wEngTP3G2r66ADt9+vSKFStkMllQUBAhRC6X19TU7N+/f+LEiX2alxBC1kZs05m0lZrbhBClrs7A6s3bIp7IDr52sQNLQpJldvGFeD9J4Cv+swgh5epreY0nzdvwLB4UQMnJyVu2bPnggw8IITKZjBBy8+bNd999d+fOndTSATxEpa8/3ZCV35gb7hw1R7ZoiNvP7ObaDsz66gLs7bffzsjICA8P72ypqalJTEw0D7V+VHZ2dmVl5cMtcrm8tbV179695ofTpk0z11K3bt06depU527m9sutPx768cCdC9UMy7QxrdeEb0ZOCvMIcH97wDplVcOj+3d3HLSj/WnaXQXu+nvM3gt7OZLHOu3ffvutVqslfepBATR16lSpVNo5FT0Ad1xvL81RZdxsL4vzfmHDoE/9xP1oJ4I+1rcXYAzDmA/SqV+/fizLdrd/UVFRlwKovb3dYDAUFRWZH8bExJg74oqKis7GzvYxnhNaifabu990MO21bfWDPcMSRYvHDhpLCLlYUfjo/t0dB+1of5p2T5GXX1Mwd/JYp720tLTPZ0bldekUrl+//v3337/xxhsbNmw4c+bMJ5988vDMHNY0ZswYPz+/9PR0Ks8OXGBg9Reazp1UppsIM9V32jivSRhsQdHdu3fHjRtXU1NjiYMfPXo0ISEhNze3S/u0adN6cbTdu3fv3Llz5syZwcHBPB5PLpcfO3Zs1apVy5cvf8oj9K7/kWtrdlf94eMom/nUXMN0lKuvxbiPpB0E4Aks0f90HRW4evVqV1fXixcv5ufnf/LJJ1u3bu3DJwN4Ss2GxjTF/753NelSy4V5gYu3RH0+yedlVD92bObMmVKpNDw8XKFQTJs27Ycffti+fbuHh0fvjrZy5crs7OywsDCVSqVUKkNDQzMzM5+++uk1Pk/AJ7a05lSNtiKzPo12CgA6ug6CVigUc+fOXbt27cKFC+Pj4+/cuUMlFtiT399e/0bIavNEKU9U2XErW3W8tLVonNfz6yK2+0tkT/4ZsBerV69esmTJwxdgWVlZvTtUSEiI9e/hCJAEvh/+kZWfFAB6p2sBFBUV9c477xw+fPjHH39MTU11d3enEgvsSbuxTWd6wuA1lrDFLQXZymMNBuVU3+n/HbSc+h1JYH12cAHmKrSBPlNv0rcYmwghzYYmA6tX6usIIRK+xF3oSTsagPV0LYB27dr197///ciRIwEBAQ0NDebZyQAsR2fS5jfmnlSluwncX/KbMdJjLJ9nS18iQB/CBZh1nGn4LleVQQgxsHq1Ub3j9kZCiJAn2hi5w/6+aN5XnTpH9rqnyJt2EOCcrgWQn5/fCy+8cPny5ZiYmOXLl4eEhFCJBbZOzbT/7719DGskhDQaVF/J/2b+ROd575eGuP3MvE+LoSlHlXG2MTvKZejSkNURLtE0EwMH4ALMOl7ym/GS3wxCyE11WVrtl2sj7HmsZ7WmQs20owCCR3UtgA4cOLB9+3aNRpOYmDhhwoQNGzYsXbqUSjKwaVK+00iPseblsm+rbwx1izWvVuMvCSCE3NNWf6c8crnl4liv5z/Abe3w/3ABBgBW07UASk1NvXDhwqJFi1xdXS9fvjxq1CgUQNALAp5gpMc48/axuq+GuQ03z4p7rf3K3+/+pVpbMcV32ieD97gIXKnGBG7BBZiV8Qifz7PD2UR1Jm1mfZr5E+hWY3O28pib0J0QMsozboBT+JN+GhxF1wJIrVZLpfe/A3Z2dhaJRFaPBHbIxJouNud/V/8vPatP8JuZ7JUi5OHUgq5wAWZl4c6RS0PX0E5hEU58JxNhCSE8wncSODsLXAnBzPHwb7oWQAsWLJgxY0Ztbe0XX3xx4MCBBQsWUIkFdkNv0vmI/D6v2Ooj9p8ZMD/GfSTWr4Du4ALMyvg8gfm7aTsj4Utf9n/VvP1D05kJ3lOxLhs8qmsB9OGHH2ZkZOTn59+7d2/btm20poEGO9BubDvVkHlKlTXIZfCbA94Nd46inQi4DhdgAGA1XQsgQsi0adM6p59va2tzc3OzbiSweY0G1Qnl0fONp0d6jFsbsS1AEkg7EdgGXIBBnwtzjnQV4H8xeIwHBdClS5dSUlIaGxtfeeWV5OTklJSU6urqkpIShUJBMR/YFrm2Jkv5r5KWH+N9pm6J+tx862mDXukt9sU3X/A0cAEGfWtxyCraEYCjHqwFtnjx4nHjxqWmpsrlcvMSrCkpKRcuXKAYDmzI7Y4bOyu3f3p7Yz9J4CeD97z20Mxjn1duVWjv0Y0HHHfp0qWEhITRo0dv2rRJpVItW7YsISFh0KBBtHMBgN168AlQU1PT5s2bCSGDBw8+evSoedvBaU2ayo5b0a7P0Q7CaWVtlzPqv1Xp6172++Xy0HfEfHGXHUysyURMVLKBrVi8ePGrr746derUgwcPxsTEJCUlLViwYMCAAbRzAYDdelAACQT31x/w8fFxcXGhlIdbqjUVRxSH1kagAHoMlrCXWi5k1H9rMOl/7j9rrOcELGEBvYYLMACwsscMgu4TDMN0VlTFxcVXr14dNWpUVBTuA7IHDMsUNJ/LqP9Wynea7j9nuMd/PTq+p9XYvP3WOhPLEEIaDQ2pdz4U8ISEkBn95sZ5v0AhNHAbLsAAwMoeFEBVVVWd4w3VanXndltbWy+OGxAQoFQqCSGpqak7dux44YUX1q9f/9FHHy1atOiZM1scwzJak4YQojVpGJZRM+2EED7hY31yI2vIa8zNqv+Xn7jfgqCkIa4x3e3pLvR8N+xD8zdff7rz4YKgZf0kgYQQT5GX9eICAAB040EBpNFoLPEEf/zjHwsLC2UymVKpjIuL664AqqqqUqlUD7d0dHSYTHQGjuSo0tPrviaEMIQxMPq1194khEgFTh9H7bS/pZKfks6kPdNwIlt5rL9T+Jv93wlzjnzij/iI/cwbQp7QS+SDBb+gB317AUYIKSoqunv37qRJkzw8PMwt6enp06dPf/aoAGAfHhRAEonEEk/g6+srk8nMG52fcj9q48aNZWVlD7dUVFR0zglrZQl+MxP8ZhJCbqqvptX+Y23ENioxOELNtJ9SZeaqMga7xqwJ+yBEOoB2IrBDfXsBtn379j179sTHx69ZsyYtLW348OGEkGXLltXW1vbhswCATbPUGCCRSDR06NDQ0NDW1taDBw8uWrRo6dKlPUxrdvDgwS4tY8aM8fPzs1A8eBptxtZs5bGzDSeGe4xJidjer7fzA6ziRgAAFRZJREFUGb4a8Cs/cUDfZgM707cXYLt27SopKfHx8SktLZ09e3ZRURHmEwKALixVAMnlco1GU1VVVVlZGRAQwLJsWFjYb3/7Wws9nYXwCI9nj0slP1GrsSVbefRsQ/bP3EdtGPQHf4nsWY42wmNsXwUDeBpubm7mkdTPPffcihUr1qxZs3//ftqhAIBbLFUAEUKcnJyio6Ojo6PND9evX2+557KQMOfI14Pfop3CqlT6+qz6tIvN+fHeUz6O2umBMctgg+bPnx8XF7d8+fKkpKQ1a9bMmjVr3rx5HR0d3e1fX1/f3t7+cItWq6U1BhEArMOCBZAdEPAEjrOOVZ1OnlH/7eWWi5N8ErZF73YTutNOBNBLGzdujIuLa2hoIITweLxvvvnm8OHDXl7dVvPJycmFhYUPt9y7d08s7jqlJwDYExRAQO5pqzPqv7naVjLFd9ong/c4CzALC9i8KVOmdG4LBIL58+fPnz+/u52/+uqrLi0Ygwhg91AAObRqzZ3jdV/fUl9/yW/GouAVUr4T7UQAfW/SpElnzpyhnQIAuAUFkIO603HzeN3hak3Fy/6/XBa6Rsy3yCQIAFzAMAztCADAOSiAHM5Nddnxuq/rdPJX/Ge9NeB3Qp6IdiIAy7KJCegBwMpQADmQsrbLx+u+bjE2TfOfPc5rkgBrl4JjWLZsGe0IAMA5KIDsH0vYK61Fx+sOa02a6f6J/+U5gc/j0w4FAABAEwoge8YS9lLLheN1XxPCTu83Z6THuEeXbQcAAHBAKIDsE0vYktbCI4pDhLDT/BNHeY5H6QMAANAJBZC9MbLGH5rOZNaneQi95sgWDXWLpZ0IAACAc1AA2Q+9SZ/XePI75RGZJHhxyKpIlyG0EwEAAHAUCiB7oDNpT6uyTqiOhTlHruz//kDnQbQTAQAAcBoKINumZtpzVRmnVJmDXWPeDfswSBpKOxEAAIANQAFkq1qMTdnKY3kNOcM9xqyL+MRfIqOdCLp15cqVrKws2ikeLy4uLj4+nnYKAHBo9XqF2thm5a8vUADZHpW+Lqv+Xxeb88d5Pb85KtVb5Es7ETxBWlpaTk4OB+uMK1eulJeXczAYADiU0tYihU6OAgi6dU9bnVn/bWnrpUk+CVujd7kLPWgngqc1derUzZs3007R1f79+3/44QfaKQAAKEABZBtuqa9n1qdVdJS/6Dd9YdByJ4Ez7UQAAADPRMN0tDNthJB2pk3DdCj1dYQQZ4GLi8DVCs+OAojTWMKWthZl1qc1Gxpf9v/ligHvinhi2qEAAAD6wNe1B8vaLhNCNEwHQ5hy9VVCiL9E9k7YZis8u2MVQEUtP0S6DHUTutMO8mQMyxQ0531X/y8+T/Bz/1dHe8RhAS8AALAni4JXmDdyVRkKnXxBkFXXLXasAuh0w3dOfOchbj+jHaQnWpPmbEN2jiq9nyRwbuBiTOUMAADQ5xyrAOK4JkNDjio9ryFnqFts8oCUUKcw2okA7MHevXuTkpJopwCAbvF5AusvWOkQBZBCd09r0hJCNCZNnU7uLHQlhPQTy7gzlLhaU5GtPFbS+uN4r8kbI//oK/annQjAhmVmZiqVys6HH330kUQiIYS8/vrr9EIBQLfGe00ymPRWflL7L4CMrGF/9ecMaySEKPTyk8rj0kYnQsh478lTfafTzWYe45ytPFaruzfVd/qvgpY6C1zoRgKK2o1trkI32insQVlZ2YYNG+bNm+ft7U0I0Wq1ly9fJiiAALhKwpdK+FIrP6n9F0BCnmj9oN+bt3fc2fSK3ywujAHSmbTnm86cVB7n8/hCnvAPg/cKeALaoYAmlb4utWLL1qj/6cXPFhQUjB079vjx49On36/pt2zZsnHjRq1WK5FItmzZcujQoba2tpUrV6akpPRpao567733xo8fn5KSsmjRoilTpuTk5KSmptIOBQDcYv8FENeo9PWnVJn5TaeiXIYuDlnFElNa7ZeofoBhGRNr6vWP+/r6fvPNN50FUHp6upubGyEkNzf3yJEjly5dam9vHzlyZEJCwogRI/omMbfFxcWlp6cnJyenp6frdLqed3755ZcvXrz4cEtbW9vIkSMtGRAAKHOsAkhABHxKpQZL2GttV041ZN5UX4v3nrJp0B99xH6EkJvqMip5wM6MGDGiuLjYYDCIRKI7d+74+flVV1cTQhoaGpYvXy6VSqVSaXx8fE1NjYMUQIQQd3f3gwcP/vOf/7x+/XrPex45ckSj0Tzcsnnz5uDgYEumAwDKHKsAeiP0bTerLx+hYTrON50+1ZAlIMIpvq8sC/2NhC+92Jx/uPYAIaSdaZFr7/656lNCiIAn+HXwKjEfUx3arRPKowVN5x5uUTPtrYZmQghLWANrWHFlLiGE5ZF+4gAB79/+PBcGLw9zjuzuyHw+f/LkyadOnUpISEhLS5s9e3ZRUREh5LXXXjPvUFpaev78+d27d/f5L8Vx8+bNmzdvXs/7mAvEh1vMg6YBwI45VgHkLvS05tNVam6fbThR2Hx+mNvw14NXRLoM7fynAc4R5lv+anU1zfqWUR7jCSFSgZOIL7JmQrCyeO8pUa7DHm4xmPRKQz0hpMXQlK08Oifw14QQEU/oK+738E2hfMIPlob2fPDExMSDBw8mJCQcPXr0yJEj69atM7ezLPvZZ5/9+c9/TktL8/BwxPXjJk2adObMGdopAIBbHKsAOqXK/Jn7aPN3T5ajNWkKmvLONma3G1sn+ry0Nfp/Hi28/MUB/uIAQshNddnVttLRnnEWjQQc4SJwdXHqusbNIDKYEFKnk59ryB7v+XyvDz5+/PgVK1ZUV1dLJBIfHx9zI8Mwc+bM8fT0LCwsdHe3gTnQLYFhGNoRAIBzHKsAutRaECAJslwBdFt941zjyUstF6Jch70a8KthbsOtP7MTOCw+nz9x4sTVq1fPmjWrs/Hw4cNisfiLL76gGIy6RYsW0Y4AAJzjWAWQhTQbGs83nfm+8RRL2AneUx/7kU93/MQBIzzGWDQe2AQvke90/8RnPEhiYuKUKVN27drV2ZKXl5eVlSWTycwP9+3b13mnmONYtsyqCwwBgE1wiAKotO2SltEQQlqNzTfay9RMOyFkoHOEr7jfsxxWZ9IWtxScbzpT0VE+ynP84pBVES7R/+lBvEQ+L/nNeJYYYB/EfHGc9wu9+9kxY8ZkZWURQiZPnmwy3b+Xvra2lhCye/duBxz4DADwRPZQAKmZdhdB13EVnYys4YemMwzLEEJajS031WUK/T1ze+8KIIZlytqKLzTnXWktjHAZHO89JXlgioiHW7cAAABshj0UQFvK3/td+MdeIp/H/quQJ0oK/a15+1lmgmZY5lr7lR+bvy9uvSiTBI3xnDg/8A03oYOOKgUAALBp9lAAmVjG/AGPJehNup/aii+1XChpLZRJg0d7xP0yYH53xRYAAADYBHsogJ7e0y+31qBXlrZdKmn98ab6aphz5HD3MYmyRZ4ib0snBAAAACuw1QKoQa/cVfV7wrKEkGZD4+cVW4U8ISEkwf+XWkbzvM9Lj/2plf3f72HVLZ1Je6P9p7L2krK2y+3G1mFuw8d7TU4K/a2TwNlCvwUAAABQYasFkLfY99fBb5mIiRDyecXWWbKFniJvPuGxLNlXk9pdAfRo9aNhOm51XL/RXnZTXVajqQxzjhzsFrM0ZHV/53BM4QMAAGCvLFgAmUymkpIShUJhNBqDgoJiY2P5fH5fHZxHeKFOA83bQp4wWNrfV+xPCJFra3r+Qb1Jf1dbWaW5U9lx63bHjUa9aqDzoEjXIbMCFoQ5R2EdLrCQpqamO3fu0E7RlUqloh0BAIAOSxVAp0+fXrFihUwmCwoKIoTI5fKampr9+/dPnDjREk+nYToqNbcJIUpdnYHVm7f5hOcp8lHp6+p0coVOXqu9e1db1WhQyiQh/Z3CBjhHTPGdFiwNpbU+PDiOiIiITZs2paen0w7yGG+++SbtCAAAFFiqAHr77bczMjLCw8M7W2pqahITEwsKCh67f35+vnnetk6NjY09L11kZA0lrYUm1hTlOixXlfFT6yUTMRlZQwfT8YdbGxjWyBCjE9/FTxLQTyzzl8hGe8bNlM4LkAT1MAwIwBIWLly4cOFC2ikAAOABSxVADMN0zr5v1q9fP5Zlu9s/KyurvLz84RaNRuPq2u30hoQQtbH9x+bvWcJK+U48Hu+/vCZIBFKjyXi+6fTS0NWeQm8fsZ+YL3nGXwQAAADsj6UKoFWrVo0cOXLmzJnBwcE8Hk8ulx87dmzVqlXd7b9169YuLe+//76fX0+rlnqIvN7s/26XRrm2pri1YLBrTK+TAwAAgN2zVAG0cuXKX/ziF1lZWXK5nBASGhqamZkZEhJioafrJOQJzffDAwAAAHTHgrVCSEhIUlKS5Y7/WP4S2e/Cu36YBAAAAPAwTn9Y0tjY2POdwzdv3nRycrJanp6p1WoXFxfaKe5DmO50dHQ4OTnxeJyY5Emn0wkEAqGQE3+GDMN4eHh4eXl1t4NCobBmHurQ//QawnQH/U93qPQ/nPjNH2vgwIE7duw4fPhwdzt0dHTU1dVx5M1jWZZhGI6EIYQYDAahUMiRPzOj0SgQCLgThs/n9+GUVM+CYRgej8eRMCaTSSwWd7l3oYuYGEcZXYf+51mg/+kO+p/u0Ol/WJuVl5cXHx9PO8V9crlcJpPRTvGAs7OzWq2mneK+YcOGlZaW0k5x38yZM48cOUI7xX2/+c1v/vSnP9FOcd/OnTtXrVpFO4XNQP/TA/Q/3UH/0x0q/Q8nSj8AAAAAa0IBBAAAAA4HBRAAAAA4HBRAAAAA4HBsuAASCoUCAVdW9eJUGEKISCTiyNh+Qgh37rQkCNM9rp3DHMepl4tTYQj6n+4hTHeonMM8tvv1uTjOZDIpFIrAwEDaQe67d+9eUFAQ7RT3cS1MYGAgR25Draur8/Hx4ciffVNTk1gs5sgkJR0dHVqt1tvbm3YQ24D+pwdcC4P+57HQ/9hwAQQAAADQO1z5lBIAAADAalAAAQAAgMNBAQQAAAAOBwUQAAAAOBwUQAAAAOBwUAABAACAw7HVAmj37t1DhgyJiYk5ffo0rQxbtmwZMmRISEjI9u3buZPq008/zczMpJ7nxIkTo0eP7t+//969e6mHWbt27cCBA8PDw//xj39QDKNWq/fs2dP58NEM1kzVJQw3T2bO4sKLw823DP3Po9D/PDEMtZPZyqvP94nq6uro6Gi1Wl1eXh4VFcUwjPUz5OTkjBgxQqPRKJXK0NDQoqIiLqTKz88Xi8UHDhxgqb5KTU1NgwYNUiqVLS0tUVFRLS0tFMOcPXs2NjZWq9XW1tZ6eHi0t7dTCXPlypUlS5bMmzfP/PDRDNZM1SUMN09mzuLCi8PNtwz9z6PQ/zwxDMWT2SY/AcrKypo5c6azs3NERERAQEBJSYn1MzQ0NCxfvlwqlfr6+sbHx9fU1FBP1drampKSMm/ePPNDinmOHj06Y8YMX19fd3f30tJSNzc3imEEAoFYLBaJRFKp1DwBK5Uwhw4dam5u7nz4aAZrpuoShoMnM5dx4cXh4FuG/uex0P88MQzFk9kmC6Da2trg4GDzdnBwsEKhsH6G1157LSkpiRBSWlp6/vz5SZMmUU+VnJy8adOmzqnEKeapqamprKx87rnnQkNDd+zYwePxKIaJi4uLjIwMDAwMCQnZvHmzi4sLlTDbtm1buXJl58NHM1gzVZcwHDyZuYwLLw4H3zL0P4+F/ueJYSiezDZZAJlMpocXdjEajVRisCybmpo6e/bstLQ0Dw8Puqm++uorb2/vKVOmdLZQzKPRaMrLy8+dO3f58uV9+/YVFRVRDHPq1Kny8vIzZ85kZmampqbK5XIunD+PZqCbilMnM8dx5MXh1FuG/qc76H+eBq2TmRNLsv2nAgMDq6urzdtyuZzKeoQMw8yZM8fT07OwsNDd3Z16qgMHDty+fTsnJ6e2tvbrr7/W6XQU8/j5+SUkJHh5eRFCJkyYcP36dYphMjMzFy5cGB0dHR0dPXbs2HPnznHh/Hk0A8VUXDuZOY4LLw7X3jL0P91B//NENE9mC40tsqiqqqphw4bpdLqampqIiAgqw/2+/PLLuXPnci0Vy7Jr1qwxD0KkmOfatWtDhw5taGior68PDQ29evUqxTB79+595ZVXOkfYFRcX0wqTk5PTOe7v0QxWTvVwGM6ezNzEhReHs28Z+p8u0P88MQzFk9kmPwEKDQ1966234uPjCSF//etf+XwKX+Tl5eVlZWXJZDLzw3379k2fPp16qodRfJWio6OTkpLGjBnDsuy6desGDx5MCKEVZsmSJcXFxUOGDCGEvPfee7GxsRTDdHr03aH4fnH/ZOYU9D9PA/2PGfqfJ6J4MvNYlrXQoQEAAAC4CRd2AAAA4HBQAAEAAIDDQQEEAAAADgcFEAAAADgcFEAAAADgcFAAAQAAgMNBAQQAAAAOBwUQAAAAOBwUQAAAAOBwUAABAACAw0EBBAAAAA4HBRAAAAA4HBRAAAAA4HBQAAEAAIDDQQEEAAAADgcFEAAAADgcFEAAAADgcIS0A4DNmzVrVnFxsVarbW5uDggIIITMnz+/ra1t586dtKMBgJ1D/wO9xmNZlnYGsAdnz579+OOPT548SQjRarVardbT05N2KABwCOh/oBfwFRj0vdzc3A8++CA3N3fq1KnPPfdcSEjIp59+OmPGjMGDB2/ZssW8z+bNmyMiImJjYw8cOEA1LADYFfQ/8JTwFRhY0I0bN65evXrt2rWxY8f+9NNP/fv3Dw8PX7du3YkTJ/Lz869cuaLT6caOHRsbGxsbG0s7LADYFfQ/0DMUQGBBkydPdnNzi4qKioyMHDJkCCHEy8vLaDTm5eUpFIq5c+cSQjQaTXFxMTogAOhb6H+gZyiAwIIEAoF5QyQSPdzu5OT0zjvvLF68mBDS3t4ukUgohAMAu4b+B3qGMUBAweTJk//2t7/p9fqWlpZhw4YpFAraiQDAUaD/ATN8AgQUTJgw4cUXX4yJiTEYDOvXrw8JCaGdCAAcBfofMMNt8AAAAOBw8BUYAAAAOBwUQAAAAOBwUAABAACAw0EBBAAAAA4HBRAAAAA4HBRAAAAA4HBQAAEAAIDDQQEEAAAADgcFEAAAADgcFEAAAADgcFAAAQAAgMNBAQQAAAAOBwUQAAAAOBwUQAAAAOBwUAABAACAw/k/uBfMiaBLr+IAAAAASUVORK5CYII=" /><!-- --></p>
</div>
<div id="background" class="section level1">
<h1>Background</h1>
<p>The <code>mkin</code> package <span class="citation">(J. Ranke
2021)</span> implements the approach to degradation kinetics recommended
in the kinetics report provided by the FOrum for Co-ordination of
pesticide fate models and their USe <span class="citation">(FOCUS Work
Group on Degradation Kinetics 2006, 2014)</span>. It covers data series
describing the decline of one compound, data series with transformation
products (commonly termed metabolites) and data series for more than one
compartment. It is possible to include back reactions. Therefore,
equilibrium reactions and equilibrium partitioning can be specified,
although this often leads to an overparameterisation of the model.</p>
<p>When the first <code>mkin</code> code was published in 2010, the most
commonly used tools for fitting more complex kinetic degradation models
to experimental data were KinGUI <span class="citation">(Schäfer et al.
2007)</span>, a MATLAB based tool with a graphical user interface that
was specifically tailored to the task and included some output as
proposed by the FOCUS Kinetics Workgroup, and ModelMaker, a general
purpose compartment based tool providing infrastructure for fitting
dynamic simulation models based on differential equations to data.</p>
<p>The ‘mkin’ code was first uploaded to the BerliOS development
platform. When this was taken down, the version control history was
imported into the R-Forge site (see <em>e.g.</em> <a href="https://cgit.jrwb.de/mkin/commit/?id=30cbb4092f6d2d3beff5800603374a0d009ad770">the
initial commit on 11 May 2010</a>), where the code is still being
updated.</p>
<p>At that time, the R package <code>FME</code> (Flexible Modelling
Environment) <span class="citation">(Soetaert and Petzoldt 2010)</span>
was already available, and provided a good basis for developing a
package specifically tailored to the task. The remaining challenge was
to make it as easy as possible for the users (including the author of
this vignette) to specify the system of differential equations and to
include the output requested by the FOCUS guidance, such as the <span class="math inline">\(\chi^2\)</span> error level as defined in this
guidance.</p>
<p>Also, <code>mkin</code> introduced using analytical solutions for
parent only kinetics for improved optimization speed. Later, Eigenvalue
based solutions were introduced to <code>mkin</code> for the case of
linear differential equations (<em>i.e.</em> where the FOMC or DFOP
models were not used for the parent compound), greatly improving the
optimization speed for these cases. This, has become somehow obsolete,
as the use of compiled code described below gives even faster execution
times.</p>
<p>The possibility to specify back-reactions and a biphasic model
(SFORB) for metabolites were present in <code>mkin</code> from the very
beginning.</p>
<div id="derived-software-tools" class="section level2">
<h2>Derived software tools</h2>
<p>Soon after the publication of <code>mkin</code>, two derived tools
were published, namely KinGUII (developed at Bayer Crop Science) and
CAKE (commissioned to Tessella by Syngenta), which added a graphical
user interface (GUI), and added fitting by iteratively reweighted least
squares (IRLS) and characterisation of likely parameter distributions by
Markov Chain Monte Carlo (MCMC) sampling.</p>
<p>CAKE focuses on a smooth use experience, sacrificing some flexibility
in the model definition, originally allowing only two primary
metabolites in parallel. The current version 3.4 of CAKE released in May
2020 uses a scheme for up to six metabolites in a flexible arrangement
and supports biphasic modelling of metabolites, but does not support
back-reactions (non-instantaneous equilibria).</p>
<p>KinGUI offers an even more flexible widget for specifying complex
kinetic models. Back-reactions (non-instantaneous equilibria) were
supported early on, but until 2014, only simple first-order models could
be specified for transformation products. Starting with KinGUII version
2.1, biphasic modelling of metabolites was also available in
KinGUII.</p>
<p>A further graphical user interface (GUI) that has recently been
brought to a decent degree of maturity is the browser based GUI named
<code>gmkin</code>. Please see its <a href="https://pkgdown.jrwb.de/gmkin/">documentation page</a> and <a href="https://pkgdown.jrwb.de/gmkin/articles/gmkin_manual.html">manual</a>
for further information.</p>
<p>A comparison of scope, usability and numerical results obtained with
these tools has been recently been published by <span class="citation">Johannes Ranke, Wöltjen, and Meinecke
(2018)</span>.</p>
</div>
</div>
<div id="unique-features" class="section level1">
<h1>Unique features</h1>
<p>Currently, the main unique features available in <code>mkin</code>
are</p>
<ul>
<li>the <a href="https://pkgdown.jrwb.de/mkin/articles/web_only/compiled_models.html">speed
increase</a> by using compiled code when a compiler is present,</li>
<li>parallel model fitting on multicore machines using the <a href="https://pkgdown.jrwb.de/mkin/reference/mmkin.html"><code>mmkin</code>
function</a>,</li>
<li>the estimation of parameter confidence intervals based on
transformed parameters (see below) and</li>
<li>the possibility to use the <a href="https://pkgdown.jrwb.de/mkin/reference/sigma_twocomp.html">two-component
error model</a></li>
</ul>
<p>The iteratively reweighted least squares fitting of different
variances for each variable as introduced by <span class="citation">Gao
et al. (2011)</span> has been available in mkin since <a href="https://pkgdown.jrwb.de/mkin/news/index.html#mkin-0-9-22-2013-10-26">version
0.9-22</a>. With <a href="https://pkgdown.jrwb.de/mkin/news/index.html#mkin-0-9-49-5-2019-07-04">release
0.9.49.5</a>, the IRLS algorithm has been complemented by direct or
step-wise maximisation of the likelihood function, which makes it
possible not only to fit the variance by variable error model but also a
<a href="https://pkgdown.jrwb.de/mkin/reference/sigma_twocomp.html">two-component
error model</a> inspired by error models developed in analytical
chemistry <span class="citation">(Johannes Ranke and Meinecke
2019)</span>.</p>
</div>
<div id="internal-parameter-transformations" class="section level1">
<h1>Internal parameter transformations</h1>
<p>For rate constants, the log transformation is used, as proposed by
Bates and Watts <span class="citation">(1988, 77, 149)</span>.
Approximate intervals are constructed for the transformed rate constants
<span class="citation">(compare Bates and Watts 1988, 135)</span>,
<em>i.e.</em> for their logarithms. Confidence intervals for the rate
constants are then obtained using the appropriate backtransformation
using the exponential function.</p>
<p>In the first version of <code>mkin</code> allowing for specifying
models using formation fractions, a home-made reparameterisation was
used in order to ensure that the sum of formation fractions would not
exceed unity.</p>
<p>This method is still used in the current version of KinGUII (v2.1
from April 2014), with a modification that allows for fixing the pathway
to sink to zero. CAKE uses penalties in the objective function in order
to enforce this constraint.</p>
<p>In 2012, an alternative reparameterisation of the formation fractions
was proposed together with René Lehmann <span class="citation">(J. Ranke
and Lehmann 2012)</span>, based on isometric logratio transformation
(ILR). The aim was to improve the validity of the linear approximation
of the objective function during the parameter estimation procedure as
well as in the subsequent calculation of parameter confidence intervals.
In the current version of mkin, a logit transformation is used for
parameters that are bound between 0 and 1, such as the g parameter of
the DFOP model.</p>
<div id="confidence-intervals-based-on-transformed-parameters" class="section level2">
<h2>Confidence intervals based on transformed parameters</h2>
<p>In the first attempt at providing improved parameter confidence
intervals introduced to <code>mkin</code> in 2013, confidence intervals
obtained from FME on the transformed parameters were simply all
backtransformed one by one to yield asymmetric confidence intervals for
the backtransformed parameters.</p>
<p>However, while there is a 1:1 relation between the rate constants in
the model and the transformed parameters fitted in the model, the
parameters obtained by the isometric logratio transformation are
calculated from the set of formation fractions that quantify the paths
to each of the compounds formed from a specific parent compound, and no
such 1:1 relation exists.</p>
<p>Therefore, parameter confidence intervals for formation fractions
obtained with this method only appear valid for the case of a single
transformation product, where currently the logit transformation is used
for the formation fraction.</p>
<p>The confidence intervals obtained by backtransformation for the cases
where a 1:1 relation between transformed and original parameter exist
are considered by the author of this vignette to be more accurate than
those obtained using a re-estimation of the Hessian matrix after
backtransformation, as implemented in the FME package.</p>
</div>
<div id="parameter-t-test-based-on-untransformed-parameters" class="section level2">
<h2>Parameter t-test based on untransformed parameters</h2>
<p>The standard output of many nonlinear regression software packages
includes the results from a test for significant difference from zero
for all parameters. Such a test is also recommended to check the
validity of rate constants in the FOCUS guidance <span class="citation">(FOCUS Work Group on Degradation Kinetics 2014,
96ff)</span>.</p>
<p>It has been argued that the precondition for this test, <em>i.e.</em>
normal distribution of the estimator for the parameters, is not
fulfilled in the case of nonlinear regression <span class="citation">(J.
Ranke and Lehmann 2015)</span>. However, this test is commonly used by
industry, consultants and national authorities in order to decide on the
reliability of parameter estimates, based on the FOCUS guidance
mentioned above. Therefore, the results of this one-sided t-test are
included in the summary output from <code>mkin</code>.</p>
<p>As it is not reasonable to test for significant difference of the
transformed parameters (<em>e.g.</em> <span class="math inline">\(log(k)\)</span>) from zero, the t-test is
calculated based on the model definition before parameter
transformation, <em>i.e.</em> in a similar way as in packages that do
not apply such an internal parameter transformation. A note is included
in the <code>mkin</code> output, pointing to the fact that the t-test is
based on the unjustified assumption of normal distribution of the
parameter estimators.</p>
</div>
</div>
<div id="references" class="section level1">
<h1>References</h1>
<!-- vim: set foldmethod=syntax: -->
<div id="refs" class="references csl-bib-body hanging-indent">
<div id="ref-bates1988" class="csl-entry">
Bates, D., and D. Watts. 1988. <em>Nonlinear Regression and Its
Applications</em>. Wiley-Interscience.
</div>
<div id="ref-FOCUS2006" class="csl-entry">
FOCUS Work Group on Degradation Kinetics. 2006. <em>Guidance Document on
Estimating Persistence and Degradation Kinetics from Environmental Fate
Studies on Pesticides in EU Registration. Report of the FOCUS Work Group
on Degradation Kinetics</em>. <a href="http://esdac.jrc.ec.europa.eu/projects/degradation-kinetics">http://esdac.jrc.ec.europa.eu/projects/degradation-kinetics</a>.
</div>
<div id="ref-FOCUSkinetics2014" class="csl-entry">
———. 2014. <em>Generic Guidance for Estimating Persistence and
Degradation Kinetics from Environmental Fate Studies on Pesticides in EU
Registration</em>. 1.1 ed. <a href="http://esdac.jrc.ec.europa.eu/projects/degradation-kinetics">http://esdac.jrc.ec.europa.eu/projects/degradation-kinetics</a>.
</div>
<div id="ref-gao11" class="csl-entry">
Gao, Z., J. W. Green, J. Vanderborght, and W. Schmitt. 2011.
<span>“Improving Uncertainty Analysis in Kinetic Evaluations Using
Iteratively Reweighted Least Squares.”</span> Journal. <em>Environmental
Science and Technology</em> 45: 4429–37.
</div>
<div id="ref-pkg:mkin" class="csl-entry">
Ranke, J. 2021. <em>‘<span class="nocase">mkin</span>‘:
<span>K</span>inetic Evaluation of Chemical Degradation Data</em>. <a href="https://CRAN.R-project.org/package=mkin">https://CRAN.R-project.org/package=mkin</a>.
</div>
<div id="ref-ranke2012" class="csl-entry">
Ranke, J., and R. Lehmann. 2012. <span>“Parameter Reliability in Kinetic
Evaluation of Environmental Metabolism Data - Assessment and the
Influence of Model Specification.”</span> In <em>SETAC World 20-24
May</em>. Berlin. <a href="https://jrwb.de/posters/Poster_SETAC_2012_Kinetic_parameter_uncertainty_model_parameterization_Lehmann_Ranke.pdf">https://jrwb.de/posters/Poster_SETAC_2012_Kinetic_parameter_uncertainty_model_parameterization_Lehmann_Ranke.pdf</a>.
</div>
<div id="ref-ranke2015" class="csl-entry">
———. 2015. <span>“To t-Test or Not to t-Test, That Is the
Question.”</span> In <em>XV Symposium on Pesticide Chemistry 2-4
September 2015</em>. Piacenza. <a href="https://jrwb.de/posters/piacenza_2015.pdf">https://jrwb.de/posters/piacenza_2015.pdf</a>.
</div>
<div id="ref-ranke2019" class="csl-entry">
Ranke, Johannes, and Stefan Meinecke. 2019. <span>“Error Models for the
Kinetic Evaluation of Chemical Degradation Data.”</span>
<em>Environments</em> 6 (12). <a href="https://doi.org/10.3390/environments6120124">https://doi.org/10.3390/environments6120124</a>.
</div>
<div id="ref-ranke2018" class="csl-entry">
Ranke, Johannes, Janina Wöltjen, and Stefan Meinecke. 2018.
<span>“Comparison of Software Tools for Kinetic Evaluation of Chemical
Degradation Data.”</span> <em>Environmental Sciences Europe</em> 30 (1):
17. <a href="https://doi.org/10.1186/s12302-018-0145-1">https://doi.org/10.1186/s12302-018-0145-1</a>.
</div>
<div id="ref-schaefer2007" class="csl-entry">
Schäfer, D., B. Mikolasch, P. Rainbird, and B. Harvey. 2007.
<span>“<span>KinGUI</span>: A New Kinetic Software Tool for Evaluations
According to <span>FOCUS</span> Degradation Kinetics.”</span> In
<em>Proceedings of the XIII Symposium Pesticide Chemistry</em>, edited
by Del Re A. A. M., Capri E., Fragoulis G., and Trevisan M., 916–23.
Piacenza.
</div>
<div id="ref-soetaert2010" class="csl-entry">
Soetaert, Karline, and Thomas Petzoldt. 2010. <span>“Inverse Modelling,
Sensitivity and Monte Carlo Analysis in <span>R</span> Using Package
<span>FME</span>.”</span> <em>Journal of Statistical Software</em> 33
(3): 1–28. <a href="https://doi.org/10.18637/jss.v033.i03">https://doi.org/10.18637/jss.v033.i03</a>.
</div>
</div>
</div>



<!-- code folding -->


<!-- dynamically load mathjax for compatibility with self-contained -->
<script>
  (function () {
    var script = document.createElement("script");
    script.type = "text/javascript";
    script.src  = "https://mathjax.rstudio.com/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML";
    document.getElementsByTagName("head")[0].appendChild(script);
  })();
</script>

</body>
</html>

Contact - Imprint