Stata The Stata listserver
[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

st: L-estimator


From   "Andreas Aschbacher" <[email protected]>
To   [email protected]
Subject   st: L-estimator
Date   Tue, 16 Mar 2004 15:26:25 +0100 (MET)

Dear fellows !

I have  300 values in radioactive-measurement in one column :
now I want to compute :

 1)             - Tn = sum 1 to 300 [a(i).x(i)]             0 < alpha < 0.5
                               a(i) = 0                                     
              for i <= n1 or i >= n- n1 + 1
                               a(i) = (n1 + 1 -alpha*n)/((1 -2*alpha)*n)  
for i = n1 + 1 or i = n - n1
                               a(i) = 1/((1-2*alpha)*n)                     
      for n1 + 1 < i < n - n1
                                                                            
                 with : n1 = [alpha*n],n = 300
 / in europe we say L - estimator type I/

 2)             - Tn = sum 1 to 300 [a(i).x(i)]             0 < alpha < 0.5
                               a(i) = (1/n)*[x(n1 +1)/x(i)]                 
      for i <= n1
                               a(i) = 1/n                                   
              for n1 + 1 <= i <= n - n1
                               a(i) = (1/n)*[x(n - n1)/x(i)]                
       for i >= n - n1 +1
                                                                            
                  with : n1 = [alpha*n],n = 300
/ in europe we say L - estimator type II/

 3)             - Tn = 0.25 *[x(n1 + 1) + x(n - n1)]  + 0.5*median-value    
      0 < alpha < 0.5
                                                                            
                  with : n1 = [alpha*n],n = 300

/ in europe we say trimean   /

4)             - Tn = 0.3 *[x(n1 + 1) + x(n - n1)]  + 0.4*median-value      
    0 < alpha < 0.5

/ in europe we say L - estimator type III/

5)if I create (x(i) + x(j))/2 that I get 300^2 values,I want to compute
median of this -
 we say wilcoxon-appliance,I hope you will understand my terminologie
/ in europe we say L - estimator type IIII/

please help me a little bit in Stata-syntax to compute more quickly
is there a command for all these values like lv-command for
midsummaries,spreads... ??
sincerely,andreas aschbacher

















-- 
+++ NEU bei GMX und erstmalig in Deutschland: T�V-gepr�fter Virenschutz +++
100% Virenerkennung nach Wildlist. Infos: http://www.gmx.net/virenschutz

*
*   For searches and help try:
*   http://www.stata.com/support/faqs/res/findit.html
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/



© Copyright 1996–2024 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index