3:35 p.m.
I can't seem to make this work. Here is the example from the User's Guide
<hansl>
function matrix stata_reorder (matrix se)
scalar n = rows(se)
return se[n] | se[1:n-1]
end function
open data4-1
ols 1 0 2 3 --cluster=bedrms
matrix se = $stderr
foreign language=stata --send-data
regress price sqft bedrms, vce(cluster bedrms)
matrix vcv = e(V)
gretl_export vcv "vcv.mat"
end foreign
matrix stata_vcv = mread("(a)dotdir/vcv.mat")
stata_se = stata_reorder(sqrt(diag(stata_vcv)))
matrix check = se - stata_se
print check
<\hansl>
There is something fishy about gretl_export. I can't seem to get a matrix
written to the dotdir directory (or anywhere else). The data are created
and loaded into Stata, the regression runs in the background, but vcv.mat
is nowhere to be found. I'm using version Stata 12 on Windows.
Am I missing something obvious?
Lee
--
Lee Adkins
Professor of Economics
lee.adkins(a)okstate.edu
learneconometrics.com
4:10 p.m.
On Thu, 11 Apr 2013, Lee Adkins wrote:
I can't seem to make this work. Here is the example from
the User's Guide
<hansl>
function matrix stata_reorder (matrix se)
scalar n = rows(se)
return se[n] | se[1:n-1]
end function
open data4-1
ols 1 0 2 3 --cluster=bedrms
matrix se = $stderr
foreign language=stata --send-data
regress price sqft bedrms, vce(cluster bedrms)
matrix vcv = e(V)
gretl_export vcv "vcv.mat"
end foreign
matrix stata_vcv = mread("(a)dotdir/vcv.mat")
stata_se = stata_reorder(sqrt(diag(stata_vcv)))
matrix check = se - stata_se
print check
<\hansl>
There is something fishy about gretl_export. I can't seem to get a matrix
written to the dotdir directory (or anywhere else). The data are created
and loaded into Stata, the regression runs in the background, but vcv.mat
is nowhere to be found. I'm using version Stata 12 on Windows.
It works OK with Stata 12 on Linux. I don't have Stata for
Windows so I'm not in a position to test. But the
"gretl_export" command is implemented by gretl_export.ado
(which is written into dotdir). It looks like this:
<stata>
program define gretl_export
version 8.2
local matrix `1'
local fname `2'
tempname myfile
file open `myfile' using "`fname'", write text replace
local nrows = rowsof(`matrix')
local ncols = colsof(`matrix')
file write `myfile' %8.0g (`nrows') %8.0g (`ncols') _n
forvalues r=1/`nrows' {
forvalues c=1/`ncols' {
file write `myfile' %15.0e (`matrix'[`r',`c']) _n
}
}
file close `myfile'
end
</stata>
Perhaps you can experiment by running something like this
manually.
Note that there's no dotdir path in this function (sorry,
"program"!), but the writing to dotdir is (or should be)
achieved by gretl doing chdir() into dotdir before running the
stata commands under "foreign" -- this being the only way we
can control where stata writes its .log output.
Allin
4:14 p.m.
Excellent. I'll give it a go. Stata is horrendous to program I think.
Or maybe I'm just not much of a programmer! Same result, though....
Thanks,
Lee
On Thu, Apr 11, 2013 at 3:10 PM, Allin Cottrell <cottrell(a)wfu.edu> wrote:
On Thu, 11 Apr 2013, Lee Adkins wrote:
> I can't seem to make this work. Here is the example from
> the User's Guide
>
> <hansl>
> function matrix stata_reorder (matrix se)
> scalar n = rows(se)
> return se[n] | se[1:n-1]
> end function
>
> open data4-1
> ols 1 0 2 3 --cluster=bedrms
> matrix se = $stderr
>
> foreign language=stata --send-data
> regress price sqft bedrms, vce(cluster bedrms)
> matrix vcv = e(V)
> gretl_export vcv "vcv.mat"
> end foreign
>
> matrix stata_vcv = mread("(a)dotdir/vcv.mat")
> stata_se = stata_reorder(sqrt(diag(stata_vcv)))
> matrix check = se - stata_se
> print check
> <\hansl>
>
> There is something fishy about gretl_export. I can't seem to get a
matrix
> written to the dotdir directory (or anywhere else). The data are created
> and loaded into Stata, the regression runs in the background, but vcv.mat
> is nowhere to be found. I'm using version Stata 12 on Windows.
It works OK with Stata 12 on Linux. I don't have Stata for
Windows so I'm not in a position to test. But the
"gretl_export" command is implemented by gretl_export.ado
(which is written into dotdir). It looks like this:
<stata>
program define gretl_export
version 8.2
local matrix `1'
local fname `2'
tempname myfile
file open `myfile' using "`fname'", write text replace
local nrows = rowsof(`matrix')
local ncols = colsof(`matrix')
file write `myfile' %8.0g (`nrows') %8.0g (`ncols') _n
forvalues r=1/`nrows' {
forvalues c=1/`ncols' {
file write `myfile' %15.0e (`matrix'[`r',`c']) _n
}
}
file close `myfile'
end
</stata>
Perhaps you can experiment by running something like this
manually.
Note that there's no dotdir path in this function (sorry,
"program"!), but the writing to dotdir is (or should be)
achieved by gretl doing chdir() into dotdir before running the
stata commands under "foreign" -- this being the only way we
can control where stata writes its .log output.
Allin
_______________________________________________
Gretl-devel mailing list
Gretl-devel(a)lists.wfu.edu
http://lists.wfu.edu/mailman/listinfo/gretl-devel
--
Lee Adkins
Professor of Economics
lee.adkins(a)okstate.edu
learneconometrics.com
10:33 a.m.
Well, I'm not making much progress. The .ado file works manually, but I
don't think Stata can find it when the foreign command is executed.
Actually, gretl doesn't like the foreign language=stata in Windows at all.
So, I don't think the problem is with gretl_export.
In fact, this generates the error (Windows) "The System cannot find file
specified."
<>
foreign language=stata
end foreign
</>
It's possible that Stata is being confused by the forward slashes in the
adopath +? This is the only statement that appears in the gretltmp.do.
Stata normally looks for backslashes under windows. But I'm just
guessing. Anyone else with Windows using Stata?
Lee
On Thu, Apr 11, 2013 at 3:14 PM, Lee Adkins <lee.adkins(a)okstate.edu> wrote:
Excellent. I'll give it a go. Stata is horrendous to program I
think.
Or maybe I'm just not much of a programmer! Same result, though....
Thanks,
Lee
On Thu, Apr 11, 2013 at 3:10 PM, Allin Cottrell <cottrell(a)wfu.edu> wrote:
> On Thu, 11 Apr 2013, Lee Adkins wrote:
>
> > I can't seem to make this work. Here is the example from
> > the User's Guide
> >
> > <hansl>
> > function matrix stata_reorder (matrix se)
> > scalar n = rows(se)
> > return se[n] | se[1:n-1]
> > end function
> >
> > open data4-1
> > ols 1 0 2 3 --cluster=bedrms
> > matrix se = $stderr
> >
> > foreign language=stata --send-data
> > regress price sqft bedrms, vce(cluster bedrms)
> > matrix vcv = e(V)
> > gretl_export vcv "vcv.mat"
> > end foreign
> >
> > matrix stata_vcv = mread("(a)dotdir/vcv.mat")
> > stata_se = stata_reorder(sqrt(diag(stata_vcv)))
> > matrix check = se - stata_se
> > print check
> > <\hansl>
> >
> > There is something fishy about gretl_export. I can't seem to get a
> matrix
> > written to the dotdir directory (or anywhere else). The data are
> created
> > and loaded into Stata, the regression runs in the background, but
> vcv.mat
> > is nowhere to be found. I'm using version Stata 12 on Windows.
>
> It works OK with Stata 12 on Linux. I don't have Stata for
> Windows so I'm not in a position to test. But the
> "gretl_export" command is implemented by gretl_export.ado
> (which is written into dotdir). It looks like this:
>
> <stata>
> program define gretl_export
> version 8.2
> local matrix `1'
> local fname `2'
> tempname myfile
> file open `myfile' using "`fname'", write text replace
> local nrows = rowsof(`matrix')
> local ncols = colsof(`matrix')
> file write `myfile' %8.0g (`nrows') %8.0g (`ncols') _n
> forvalues r=1/`nrows' {
> forvalues c=1/`ncols' {
> file write `myfile' %15.0e (`matrix'[`r',`c']) _n
> }
> }
> file close `myfile'
> end
> </stata>
>
> Perhaps you can experiment by running something like this
> manually.
>
> Note that there's no dotdir path in this function (sorry,
> "program"!), but the writing to dotdir is (or should be)
> achieved by gretl doing chdir() into dotdir before running the
> stata commands under "foreign" -- this being the only way we
> can control where stata writes its .log output.
>
> Allin
>
>
> _______________________________________________
> Gretl-devel mailing list
> Gretl-devel(a)lists.wfu.edu
> http://lists.wfu.edu/mailman/listinfo/gretl-devel
>
--
Lee Adkins
Professor of Economics
lee.adkins(a)okstate.edu
learneconometrics.com
--
Lee Adkins
Professor of Economics
lee.adkins(a)okstate.edu
learneconometrics.com
1:42 p.m.
On Fri, 12 Apr 2013, Lee Adkins wrote:
Well, I'm not making much progress. The .ado file works
manually, but I
don't think Stata can find it when the foreign command is executed.
Actually, gretl doesn't like the foreign language=stata in Windows at all.
So, I don't think the problem is with gretl_export.
In fact, this generates the error (Windows) "The System cannot find file
specified."
<>
foreign language=stata
end foreign
</>
It's possible that Stata is being confused by the forward slashes in the
adopath +?
Possibly. Windows itself (cmd.exe) can handle forward slashes
OK, but maybe something is going wrong between Stata and
Windows. (Things work OK using forward slashes with R and Ox
on Windows, in the analogous context.) In CVS and snapshot I'm
trying putting the slashes the other way in the "adopath +"
command for Stata.
Allin
5:31 a.m.
New subject: --output=display in corrgm and pergm
I am missing an option '--output=display' in comands 'corrgm' and
'pergm', so that we can put them
in the menus built with a bundle-plot. Allin please, could you take a
look to see whether it would be possibly to add it?
--
Firma Arista
*Ignacio Díaz-Emparanza*
Zuzendaria/Director
ignacio.diaz-emparanza(a)ehu.es <mailto:ignacio.diaz.emparanza@ehu.es>
94 6013732
*EKONOMIA APLIKATUA III SAILA (EKONOMETRIA ETA ESTATISTIKA)/
DEPARTAMENTO DE ECONOMÍA APLICADA III (ECONOMETRÍA Y ESTADÍSTICA)
UPV/EHU*
Avda. Lehendakari Aguirre, 83 | 48015 BILBAO
*T.: +34 946013740* | *F.: +34 946013754*
*www.ea3.ehu.es* <http://www.ea3.ehu.es>
ERNE! Baliteke mezu honen zatiren bat edo mezu osoa legez babestuta
egotea. Mezuak badu bere hartzailea. Okerreko helbidera heldu bada
(helbidea gaizki idatzi, transmisioak huts egin) eman abisu igorleari,
korreo honi erantzunda. Kontuz! Mezua ez bada zuretzat, ez erabili, ez
zabaldu beste inori, ez kopiatu eta ez baliatu.
¡ATENCIÓN! Este mensaje contiene información privilegiada o confidencial
a la que sólo tiene derecho a acceder el destinatario. Si usted lo
recibe por error le agradeceríamos que no hiciera uso de la información
y que se pusiese en contacto con el remitente.
E-mail hau inprimatu baino lehen egiaztatu inprimatzeko beharra.
Antes de imprimir este e-mail piense bien si es necesario hacerlo.
1:07 p.m.
Good news. That fixed it in Windows. Thanks!
Lee
On Fri, Apr 12, 2013 at 12:42 PM, Allin Cottrell <cottrell(a)wfu.edu> wrote:
On Fri, 12 Apr 2013, Lee Adkins wrote:
> Well, I'm not making much progress. The .ado file works manually, but I
> don't think Stata can find it when the foreign command is executed.
>
> Actually, gretl doesn't like the foreign language=stata in Windows at
all.
> So, I don't think the problem is with gretl_export.
>
> In fact, this generates the error (Windows) "The System cannot find file
> specified."
>
> <>
> foreign language=stata
> end foreign
> </>
>
> It's possible that Stata is being confused by the forward slashes in the
> adopath +?
Possibly. Windows itself (cmd.exe) can handle forward slashes
OK, but maybe something is going wrong between Stata and
Windows. (Things work OK using forward slashes with R and Ox
on Windows, in the analogous context.) In CVS and snapshot I'm
trying putting the slashes the other way in the "adopath +"
command for Stata.
Allin
_______________________________________________
Gretl-devel mailing list
Gretl-devel(a)lists.wfu.edu
http://lists.wfu.edu/mailman/listinfo/gretl-devel
--
Lee Adkins
Professor of Economics
lee.adkins(a)okstate.edu
learneconometrics.com
4189
days inactive
4194
days old
6 comments
3 participants
participants (3)
-
Allin Cottrell -
Ignacio Diaz-Emparanza -
Lee Adkins