s-news
[Top] [All Lists]

Re: Another S+ Controlling Question

To: Cheng Peng <chpeng@math.utoledo.edu>
Subject: Re: Another S+ Controlling Question
From: John Fox <jfox@mcmaster.ca>
Date: Sat, 16 Nov 2002 09:23:23 -0500
Cc: s-news@lists.biostat.wustl.edu
In-reply-to: <002701c28d33$106d6b10$4365b783@utftyatg7fcqkx>
Dear Cheng Peng,

At 11:43 PM 11/15/2002 -0600, Cheng Peng wrote:

I am using S+ to bootstrap the p-value of a proposed test statistic without analytic distribution. For each bootstrap sample, I have to find solutions to the score equations. Here comes the problem:

I use Quasi-Newton algorithm, for some of the bootstrap samples, when the s+ calls Fortran function, I got following message from the system and program cracked down:
-----------------------------------------------------------------------------
Error in .Fortran(if(!cmplx) "dqr" else "zqr",: subroutine dqr: 25 missing value(s) in argument 1
-----------------------------------------------------------------------------
When I use try() function and try to terminate the current loop to generate next bootstrap sample, the system gives another message:
--------------------------------------------------------
the object cmplx cannot be found.
--------------------------------------------------------
Is there any control statement or function can be used here in case we confront such error, we can terminate the current loop and go to next bootstrap sample without re-running the whole program.


It's a bit hard for me to tell from your message what the nature of the problem is, but I'm going to guess that if the current bootstrap sample isn't processed successfully, for some reason emplx doesn't exist. Is that right? If so, you can use exists("emplx") to decide whether to continue.

John

-----------------------------------------------------
John Fox
Department of Sociology
McMaster University
Hamilton, Ontario, Canada L8S 4M4
email: jfox@mcmaster.ca
phone: 905-525-9140x23604
web: www.socsci.mcmaster.ca/jfox
-----------------------------------------------------


<Prev in Thread] Current Thread [Next in Thread>