Re: [Gretl-devel] regarding the new modprint command...
On Thu, 29 Jan 2009, Riccardo (Jack) Lucchetti wrote:
On Wed, 28 Jan 2009, andreas.rosenblad(a)ltv.se wrote:
>
>
> Regarding the new modprint command...
>
> Rlease
> 1. Document it in the user's manual.
A brief description is already in the reference manual and in the online
help. But you're right, an example would be helpful.
Have a look here:
<script>
set echo off
set messages off
function my_own_ols(series y, series x)
s2 = 1
b0 = 0
b1 = 0
mle ll = -0.5*(ln(s2) + e^2/s2)
series e = y - b0 - b1*x
params b0 b1 s2
end mle -q
printf "\n ------ My own ols (via ML!) --------------------------"
table = $coeff ~ $stderr
table = table[1:2,]
names = "constant,x,sigma"
matrix sigma = sqrt(s2)
modprint table names sigma
end function
nulldata 10
set seed 12345
x = normal()
y = normal()
ols y 0 x --simple-print
my_own_ols(y,x)
</script>
Is this too silly to package as a practice .inp file? Or maybe it could
find its way to the User's Manual?
Riccardo (Jack) Lucchetti
Dipartimento di Economia
Università Politecnica delle Marche
r.lucchetti(a)univpm.it
http://www.econ.univpm.it/lucchetti