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/