Replay of saved session does not work
by Marco Dieckhoff
Hi!
I have a problem, a saved session will not replay what I did manually.
My step-by-step approach:
Open CSV File (Colon, Point)
Yes, use as time series
Time Series
Options: Yeary, Quarterly, Monthly, Weekly, Daily (5, 6, 7), Hourly, Decadic
(translated from german, may not be exact, also on following terms).
I'm missing "daily (1)"?
So I just choose "Other: 24" for my purposes. That interprets as
"Hourly, 1:01 to 1:24". But it's not the problem I wanted to discuss:
After importing, I do
Model - Time Series - ARIMA
temp as dependent, hour as non-dependend variable
All other settings on standard: AR 1, I 0, MA 1, Seasonals all on 0.
Use constant, do not show iteration details, do Hesse matrix, no
X-12-Arima, exact maximum likelihood
OK =>
function calculations: 91, gradient: 37, ...
In this window, I use the menu to:
Analysis, Prognosis, add 6
automatic prognosis
plot previous 24
OK
And I get the graph I want to see.
So far, so good.
Now I save this session in gretl with File - Session - Save Session.
I open the resulting file in WinRar and get the .INP file.
# Log started 2009-10-19 08:57
# recording of session commands. please not this must most likely
# be edited before executing
open C:\TempData1.csv
setobs 24 1:01 --special-time-series
# Modell 1
arima 1 0 1 ; 2 1 0
fcasterr 2:01 2:06
gnuplot temp --with-lines --time-series
And I get the error
? fcasterr 2:01 2:06
Observation Number out of boundaries (Beobachtungsnummer außerhalb der
Grenzen)
What am I missing, I thought the session should at least replay?
Thanks in advance for you help,
marco
"hour";"temp"
0;36.8850403940694
1;35.8449242163804
2;35.1328048127208
3;34.5454214903154
4;33.2165294935116
5;32.566773832396
6;31.908624056118
7;32.5068833835506
8;33.781475965867
9;34.823469472098
10;35.8511630975871
11;37.0685152179023
12;36.9298478305651
13;36.1227166129879
14;35.3174865764317
15;34.5371106629089
16;33.3070930219017
17;32.4677509971601
18;31.4553978170153
19;32.1524502276138
20;33.4864781189783
21;34.5596741333415
22;35.3637751201174
23;37.0507264496705
14 years, 11 months
concerning the commands for double bounded DC in contingent valuation
by Dorian Litvine
Hi,
My name is Dorian, I'm trying to estimate the data of a double bounded
dichotomous choice (contingent valuation) thanks to a maximum
likelihood estimation with a linear specification (From E. Flachaire).
But when I put the command, then GRETL gives several errors messages
as "missing values found". I think the model may not be specified for
linear. But it is strange because when I put initial values quite high
(5) then it converges but I don't know if results are "possible". May
someone help me? Here is the code
genr alpha=3
genr sigma=1
logl=(1-CD1)*(1-CD2)*log(Pnn)+(1-CD1)*CD2*log(Pny)+CD1*(1-CD2)*log(Pyn)+CD1*CD2*log(Pyy)
series Pnn=1/(1+exp(-(BID2-alpha)/sigma))
series
Pny=1/(1+exp(-(BID1-alpha)/sigma))-1/(1+exp(-(BID2-alpha)/sigma))+(BID2>BID1)*100
series
Pyn=1/(1+exp(-(BID2-alpha)/sigma))-1/(1+exp(-(BID1-alpha)/sigma))+(BID2<BID1)*100
series Pyy=1-1/(1+exp(-(BID2-alpha)/sigma))
params alpha sigma
And attached you have the results if I set the initial values at 5.
Thank you very much !!
--
Dorian LITVINE
Docteur en Sciences Economiques
Chercheur (associé) en Economie de l'Energie et de l'Environnement
Research Associate in Economics (PhD), specialized in
Energy/Environment, Green
Electricity, and pro-ecological decisions
Laboratory LASER-CREDEN (Energy Economics Institute)
Faculté des Sciences Economiques CS 79606
34960 Montpellier Cedex 2 -- FRANCE
Tel : (office) 00 33 (0)4 67 15 83 32
(mobile) 00 33 (0)6 89 15 15 74
E-mail : dorian.litvine(a)univ-montp1.fr / dorian.litvine(a)gmail.com
CREDEN website: www.creden.univ-montp1.fr
14 years, 11 months
Ralph M Rodriguez/PO/KAIPERM is out of the office.
by Ralph.M.Rodriguez@kp.org
I will be out of the office starting 10/20/2009 and will not return until
10/23/2009.
Hi All, I will be out of the office from tomorrow Tuesday Oct 20 through
Oct 22, 2009 returning to the office Oct 23, Monday. I will be taking my
computer during my trip and will be working on occassion during that time.
Please, if you have an immediate question regarding Cost Model or
Construction Economics, contact Bradly A Njus at 510 625 4595 or email us
with your questions.
Thanks much for your consideration.
Ralph
14 years, 11 months
Standard error for confidence intervals
by Francisco Sosa
Hi
I have sent an e-mail to the mailing list, but I haven't got so much
success, I am going to rewrite the question to make it more understandable
(I am not a native English speaker)
Basically, I have a dynamic prediction and I do not know why the standard
error grows up (see below), does anybody have the formula or any paper where
I can look into, thanks
----------------------------------------------------------------------------
---------------------
For 95% confidence intervals, t(4, 0.025) = 2.776
Obs Premium prediction std. error 95% interval
2004 302.00 301.15
2005 337.00 340.66
2006 375.00 372.89
2007 407.25 2.792 399.50 - 415.00
2008 442.78 2.846 434.88 - 450.68
2009 474.53 2.848 466.63 - 482.44
2010 508.07 2.848 500.16 - 515.97
----------------------------------------------------------------------------
-------------------------------------
The model is AR(1) with intercept, Prais-Winsten transformation
(Cochrane-Orcutt).
Thank you very much, any help would be great.
Kind regards
14 years, 11 months
SVAR
by Erhanfadli Azrai
Hello everyone.
Can the command 'system' and 'estimate' be used for structural var
equation? or alternatively, could anyone point me to a manual or
write-up to do svar using gretl?
I've read somewhere that i can do this in JMulti, but i prefer to code
in gretl.
Many thanks.
Erhan
Malaysia.
14 years, 11 months
Some details of gretl's boxplots
by Giorgio Gozzi
I am not able to obtain the printing of the median and quartiles next to the
boxplot and to obtain also a list of the outliers (points lying beyond 1.5
times the interquartile range from the central box). What type of statements
or options in the boxplot command I must write?
Thank you very much, any help would be great.
Kind regards
Giorgio Gozzi
Prof. Giorgio Gozzi
Professor of Economic Statistics
Department of Economics
Via J.F. Kennedy, 6
University of Parma
Italy
ISI Member
Phone: +390521032418
Fax: +390521032375
http://economia.unipr.it/DOCENTI/
14 years, 11 months
Ralph M Rodriguez/PO/KAIPERM is out of the office.
by Ralph.M.Rodriguez@kp.org
I will be out of the office starting 10/13/2009 and will not return until
10/19/2009.
Hi All, I will be out of the office from tomorrow Wednesday Oct 14 through
Oct 16, 2009 returning to the office Oct 19, Monday. I will be taking my
computer during my trip and will be working on occassion during that time.
Please, if you have an immediate question regarding Cost Model or
Construction Economics, contact Bradly A Njus at 510 625 4595 or email us
with your questions.
Thanks much for your consideration.
Ralph
14 years, 11 months
confidence interval for prediction PRAIS-WINSTEN
by Francisco Sosa
Hi
I have got an issue, I estimated a model, I used Prais-Winsten
autoregressive model (lag one and constant), I forecasted the value for the
next two years, see table below:
For 95% confidence intervals, t(4, 0.025) = 2.776
Obs Global_Sales prediction std. error 95% interval
2004 125.00 125.47
2005 139.00 139.01
2006 151.00 150.70
2007 166.82 0.693 164.89 -
168.74
2008 181.81 0.709 179.84 -
183.77
My question is; how the std. error is calculated, I could calculate the std.
error for 2007 (summarize of the squares of residuals, divided by n-2), but
I do not know how it is calculated for the next year (2008 - 0.709)
Does anybody know the formula or how it is calculated? I am interested in
any paper or book where I can find it.
Thank you very much
Kind regards
Francisco Sosa
dd: +44 (0) 20 7448 3875
email: francisco(a)digitab.uk.com
digitab
tel: +44 (0) 20 7588 6669
fax: +44 (0) 20 7256 7540
digitab 4th floor 6-8 bonhill street london ec2a 4bx
please consider the environment before printing this document
14 years, 11 months
Re: [Gretl-users] Gretl-users Digest, Vol 33, Issue 13
by Lars Pålsson-Syll
Re: Changing base in MN logit,
Thanks Allin.
Your suggestion, as always, worked just fine.
Maybe something for the "Cheat sheet"?
Lars
Lars Palsson Syll
Professor of Civics
Malmo University
Sweden
>>> <gretl-users-request(a)lists.wfu.edu> 09-10-12 18:05 >>>
Send Gretl-users mailing list submissions to
gretl-users(a)lists.wfu.edu
To subscribe or unsubscribe via the World Wide Web, visit
http://lists.wfu.edu/mailman/listinfo/gretl-users
or, via email, send a message with subject or body 'help' to
gretl-users-request(a)lists.wfu.edu
You can reach the person managing the list at
gretl-users-owner(a)lists.wfu.edu
When replying, please edit your Subject line so it is more specific
than "Re: Contents of Gretl-users digest..."
Today's Topics:
1. Re: Gretl-users Digest, Vol 33, Issue 12 ( Lars P?lsson-Syll )
2. Multinomial logit (wasRe: Gretl-users Digest, Vol 33, Issue
12) (Allin Cottrell)
----------------------------------------------------------------------
Message: 1
Date: Sun, 11 Oct 2009 22:18:22 +0200
From: " Lars P?lsson-Syll " <Lars.Palsson-Syll(a)mah.se>
Subject: Re: [Gretl-users] Gretl-users Digest, Vol 33, Issue 12
To: <gretl-users(a)lists.wfu.edu>
Message-ID: <4AD259AE020000F00008215B(a)biggs.mah.se>
Content-Type: text/plain; charset=US-ASCII
Hi,
I wonder if there is a possibility, somehow, in gretl to choose which alternative should count as the base category when
running a multinomial logit model. To me it seems as if gretl automatically always chooses the FIRST alternative. [In stata
one has the command baseoutcome() to enable you to choose for yourself.]
Best regards,
Lars
Lars Palsson Syll
Professor of Civics
Malmo University
Sweden
>>> <gretl-users-request(a)lists.wfu.edu> 09-10-09 18:04 >>>
Send Gretl-users mailing list submissions to
gretl-users(a)lists.wfu.edu
To subscribe or unsubscribe via the World Wide Web, visit
http://lists.wfu.edu/mailman/listinfo/gretl-users
or, via email, send a message with subject or body 'help' to
gretl-users-request(a)lists.wfu.edu
You can reach the person managing the list at
gretl-users-owner(a)lists.wfu.edu
When replying, please edit your Subject line so it is more specific
than "Re: Contents of Gretl-users digest..."
Today's Topics:
1. Re: descriptive labels (Shaianne Osterreich)
----------------------------------------------------------------------
Message: 1
Date: Thu, 8 Oct 2009 16:55:57 -0400
From: Shaianne Osterreich <shaianneo(a)gmail.com>
Subject: Re: [Gretl-users] descriptive labels
To: Gretl list <gretl-users(a)lists.wfu.edu>
Message-ID:
<3b32039f0910081355n45882f7fqc0c351423b369685(a)mail.gmail.com>
Content-Type: text/plain; charset="iso-8859-1"
Hi everyone.
Thank you all for your helpful comments. I got Allin's trick to work. Here's
a summary of what I did.
1. Open mydata
store mydata --traditional
2. I had an excel file with the variables and labels in two columns. Save
as a text file. Make sure that gretl has read the variables exactly the same
way you have them in the file. (My data file/ and description file had extra
00s which gretl wisely dropped - but this made it seem like the variables
were not the same). Open in Word and save as a text file again. This
provides a more extensive format check to make sure the spaces are spaces
(and not tabs) etc and hide the txt extension, adding .lbl. Save in the
same user folder as the data file.
3. Open mydata.dat again. This time it opens with all the variables.
So basically, Allin, it's what you said. I just had to be more careful with
what the final text file actually looked like.
Henrique - For future purposes I am still working on importing a script from
a text file with an .inp extension. I got an error which read 'error in
byte compatibility'. Not sure what that is about. I will try again.
Thanks!!
Shaianne
On Thu, Oct 8, 2009 at 5:33 AM, yinung at Gmail <yinung.cycu(a)gmail.com>wrote:
> Dear all,
>
> I tried Allin's solution by command line in Gretl 1.8.4 console under
> Windows XP.
>
> It did work to generate 3 files, e.g., *.dat, *.hdr, and *.lbl.
>
> After some experiments, I found that the working file must have at least
> _ONE_ variable with a descriptive label. Otherwise, gretl only generate two
> files (without *.lbl) so that you are not able to edit the *.lbl by general
> text editor thereafter.
>
> It is also confirmed that re-open the *,dat via GUI does not import the
> associated descriptive labels defined in *.lbl. You have to re-open the
> *.dat file trhough the gretl Console (command line).
>
>
> good luck
>
> Yi-Nung
>
>
> 2009/10/8 Allin Cottrell <cottrell(a)wfu.edu>
>
>>
>> On Wed, 7 Oct 2009, Shaianne Osterreich wrote:
>>
>> > Thanks for your advice Allin and Henrique.
>> >
>> > Allin - I've tried your suggestion. When the traditional file
>> > is created it does create 2 files - dat and hdr. Then I have
>> > saved a plain text file named mydata.lbl (in exactly the format
>> > you suggested) in the same user folder with the newly created
>> > ones. But, when I re-open mydata.dat it doesn't open with the
>> > descriptive labels - it's just the variable names as before. I
>> > think gretl isn't reading the text file. Do you have
>> > suggestions?
>>
>>
>
> _______________________________________________
> Gretl-users mailing list
> Gretl-users(a)lists.wfu.edu
> http://lists.wfu.edu/mailman/listinfo/gretl-users
>
14 years, 11 months
Re: [Gretl-users] Gretl-users Digest, Vol 33, Issue 12
by Lars Pålsson-Syll
Hi,
I wonder if there is a possibility, somehow, in gretl to choose which alternative should count as the base category when
running a multinomial logit model. To me it seems as if gretl automatically always chooses the FIRST alternative. [In stata
one has the command baseoutcome() to enable you to choose for yourself.]
Best regards,
Lars
Lars Palsson Syll
Professor of Civics
Malmo University
Sweden
>>> <gretl-users-request(a)lists.wfu.edu> 09-10-09 18:04 >>>
Send Gretl-users mailing list submissions to
gretl-users(a)lists.wfu.edu
To subscribe or unsubscribe via the World Wide Web, visit
http://lists.wfu.edu/mailman/listinfo/gretl-users
or, via email, send a message with subject or body 'help' to
gretl-users-request(a)lists.wfu.edu
You can reach the person managing the list at
gretl-users-owner(a)lists.wfu.edu
When replying, please edit your Subject line so it is more specific
than "Re: Contents of Gretl-users digest..."
Today's Topics:
1. Re: descriptive labels (Shaianne Osterreich)
----------------------------------------------------------------------
Message: 1
Date: Thu, 8 Oct 2009 16:55:57 -0400
From: Shaianne Osterreich <shaianneo(a)gmail.com>
Subject: Re: [Gretl-users] descriptive labels
To: Gretl list <gretl-users(a)lists.wfu.edu>
Message-ID:
<3b32039f0910081355n45882f7fqc0c351423b369685(a)mail.gmail.com>
Content-Type: text/plain; charset="iso-8859-1"
Hi everyone.
Thank you all for your helpful comments. I got Allin's trick to work. Here's
a summary of what I did.
1. Open mydata
store mydata --traditional
2. I had an excel file with the variables and labels in two columns. Save
as a text file. Make sure that gretl has read the variables exactly the same
way you have them in the file. (My data file/ and description file had extra
00s which gretl wisely dropped - but this made it seem like the variables
were not the same). Open in Word and save as a text file again. This
provides a more extensive format check to make sure the spaces are spaces
(and not tabs) etc and hide the txt extension, adding .lbl. Save in the
same user folder as the data file.
3. Open mydata.dat again. This time it opens with all the variables.
So basically, Allin, it's what you said. I just had to be more careful with
what the final text file actually looked like.
Henrique - For future purposes I am still working on importing a script from
a text file with an .inp extension. I got an error which read 'error in
byte compatibility'. Not sure what that is about. I will try again.
Thanks!!
Shaianne
On Thu, Oct 8, 2009 at 5:33 AM, yinung at Gmail <yinung.cycu(a)gmail.com>wrote:
> Dear all,
>
> I tried Allin's solution by command line in Gretl 1.8.4 console under
> Windows XP.
>
> It did work to generate 3 files, e.g., *.dat, *.hdr, and *.lbl.
>
> After some experiments, I found that the working file must have at least
> _ONE_ variable with a descriptive label. Otherwise, gretl only generate two
> files (without *.lbl) so that you are not able to edit the *.lbl by general
> text editor thereafter.
>
> It is also confirmed that re-open the *,dat via GUI does not import the
> associated descriptive labels defined in *.lbl. You have to re-open the
> *.dat file trhough the gretl Console (command line).
>
>
> good luck
>
> Yi-Nung
>
>
> 2009/10/8 Allin Cottrell <cottrell(a)wfu.edu>
>
>>
>> On Wed, 7 Oct 2009, Shaianne Osterreich wrote:
>>
>> > Thanks for your advice Allin and Henrique.
>> >
>> > Allin - I've tried your suggestion. When the traditional file
>> > is created it does create 2 files - dat and hdr. Then I have
>> > saved a plain text file named mydata.lbl (in exactly the format
>> > you suggested) in the same user folder with the newly created
>> > ones. But, when I re-open mydata.dat it doesn't open with the
>> > descriptive labels - it's just the variable names as before. I
>> > think gretl isn't reading the text file. Do you have
>> > suggestions?
>>
>>
>
> _______________________________________________
> Gretl-users mailing list
> Gretl-users(a)lists.wfu.edu
> http://lists.wfu.edu/mailman/listinfo/gretl-users
>
14 years, 12 months