I think I answered my own
question…there are 66 0’s (zero) in pspd. Apparently, MLE ignores the 0’s, but
they are important. I’ll have to make a work-around. Sorry to bother.
Win
From:
s-news-owner@lists.biostat.wustl.edu
[mailto:s-news-owner@lists.biostat.wustl.edu] On Behalf Of Crawford.Winnie
Sent: Wednesday, February 11, 2009 10:12 AM
To: S-PLUS Newsgroup
Subject: [S] Robust Weibull limits?
All,
I am using S-PLUS8 on a Windows XP PC.
I am trying to determine the best-fit distribution for a very
large sample size. At this time I want to compare the output from the functions
weibullRob and weibullMLE. In this specific example, the sample size is 24,775
in a vector named pspd. The function weibullMLE(pspd) works and produces output
as expected, but when I use the function
>weibullRob(pspd)
I get the error
Problem in
.Fortran("s_tmadve",: subroutine s_tmadve: 66 Inf value(s) in
argument 1
Use traceback()
to see the call stack
This is the call stack:
> traceback()
9: eval(action, sys.parent())
8: doErrorAction("Problem in
.Fortran(\"s_tmadve\",: subroutine s_tmadve: 66 Inf value(s) in
argument 1", 1000)
7: .Fortran("s_tmadve",
6: weibullRob(pspd)
5: eval(i, local)
4: source(auto.print = auto.print, exprs =
substitute(exprs.literal))
3: script.run(exprs.literal = {
2: eval(_expression_(script.run(exprs.literal = {
1:
Message: Problem in .Fortran("s_tmadve",:
subroutine s_tmadve: 66 Inf value(s) in argument 1
Could it be that the sample size is too large for the Robust
Weibull? I do not know how to interpret the traceback output. Thanks for any
help.
Win
*****************************************************************
Winifred
C. Crawford Staff Scientist/Senior Meteorologist
ENSCO,
Inc.
Aerospace
Sciences and Engineering Division
1980
N. Atlantic Ave., Suite 830
Cocoa
Beach, FL 32931
VOICE:
321.853.8130 FAX: 321.853.8415
EMAIL:
crawford.winnie@ensco.com
AMU
Quarterly Reports are available online:
http://science.ksc.nasa.gov/amu
*****************************************************************
The
information contained in this email message is intended only for the use of the
individual(s) to whom it is addressed and may contain information that is
privileged and sensitive. If you are not the intended recipient, or otherwise
have received this communication in error, please notify the sender immediately
by email at the above referenced address and note that any further
dissemination, distribution or copying of this communication is strictly
prohibited.
The U.S. Export Control Laws regulate the export and re-export of technology originating
in the United States. This includes the electronic transmission of information
and software to foreign countries and to certain foreign nationals. Recipient
agrees to abide by these laws and their regulations -- including the U.S.
Department of Commerce Export Administration Regulations and the U.S.
Department of State International Traffic in Arms Regulations -- and not to
transfer, by electronic transmission or otherwise, any content derived from
this email to either a foreign national or a foreign destination in violation
of such laws.