. program anova
args n1 m1 s1 n2 m2 s2 n3 m3 s3
scalar k=3
gen xbar=(n1*m1+n2*m2+n3*m3)/(n1+n2+n3)
gen s2w=((n1-1)*s1*s1+(n2-1)*s2*s2+(n3-1)*s3*s3)/(n1+n2+n3-k)
gen s2b=(n1*(m1-xbar)*(m1-xbar)+ n2*(m2-xbar)*(m2-xbar)+n3*(m3-xbar)*(m3-xbar))/(k-1)
gen F=s2b/s2w
display "F-value==" F
gen P= invF(k-1,N-k, .95)
display "P-value of the F-statistics is" P
end
I am trying to calculate P-value of F-statistics obtained by the formulas above. I want to give 9 arguments which will be scalar n1 m1 s1......n3 m3 s3. So if initiate my program and type >anova 1 2 3 4 5 6 7 8 9 it would display me the F-value and P value. This is my first time I am doing something like this and I cant figure out. Can someone help me? I can elaborate if my questions are not clear.
Related Posts with Writing a simple program to calculate ANOVA directly
STATA command to make the PSTR and PSTAR models.Hello I am a Master's student in economics and I would like you to help me with the STATA command to…
mlogit: change y = 0.5 to other options if opt-out is selectedHi All, My dependent variable has three options: "No preference", "Prefer A", "Prefer B". I would l…
How do I calculate the number of workers in the household?Hi, I am working on panel data, and I want to know how many workers are in the household. An exampl…
High coefficients values of the interaction termsHi experts and researchers, I am working on a Panel data model using system GMM, xtabond2 command. …
shapley returns "variable __weigh not found"Dear all, I wanted to replicate the example in the help file of the shapley.ado. But it returns "va…
Subscribe to:
Post Comments (Atom)
0 Response to Writing a simple program to calculate ANOVA directly
Post a Comment