Hi Jeff,
>From your words (if our words really describe us), I hope you don't expect
me to teach you that this is "r-help" room. I don't expect you to tell me
that I am a layman. I already know it and that is why I am here seeking a
help. If I am an expert of things I asked here, there is no need for me to
come here with such "stupid (in your understanding)" question. I know there
are people on this planet who consider themselves advanced only when they
meet stupid people like me - you may be one of them. It is always helpful
for others if you could write a single line with a helping mind than
writing tones of useless junk words. R is applied statistical language, to
be used in different professions. So, don't expect all R users to be
"experts" like you are. If you can't help people, at least keep yourself
away from helping environment.

Cheers

Bukana


On Tue, Dec 24, 2013 at 5:23 PM, Jeff Newmiller <jdnewmil@dcn.davis.ca.us>wrote:

> Once you start describing your code in language that indicates you have
> read some of the recommended materials, we can make some progress. For
> example, you are still posting in HTML format so your code is getting
> mangled (see Posting Guide). You don't seem to understand what tryCatch
> does (see ?tryCatch; hint... it does not alter the result returned by
> ca.jo, nor does it guarantee that any result will be returned). For that
> matter, you don't seem to understand how to give valid inputs to ca.jo,
> but then you don't seem to understand how to provide a reproducible example
> either (see for example
> http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example
> ).
>
> In short, if you ask us why one divided by zero doesn't give 3, we have to
> wonder if you don't belong in some other educational forum, because this is
> not a mathematics theory support group... it is about R. Please read or
> otherwise absorb the recommended background materials mentioned here and
> then ask clear questions here about R, or find some one-on-one help offline.
> ---------------------------------------------------------------------------
> Jeff Newmiller                        The     .....       .....  Go Live...
> DCN:<jdnewmil@dcn.davis.ca.us>        Basics: ##.#.       ##.#.  Live
> Go...
>                                       Live:   OO#.. Dead: OO#..  Playing
> Research Engineer (Solar/Batteries            O.O#.       #.O#.  with
> /Software/Embedded Controllers)               .OO#.       .OO#.  rocks...1k
> ---------------------------------------------------------------------------
> Sent from my phone. Please excuse my brevity.
>
> mamush bukana <mamushbukana@gmail.com> wrote:
> >Dear Pat,
> >I tried your suggestion:
> >for(i in 1:N1){
> >for(j in 1:N2){
> >co<-tryCatch(ca.jo <http://ca.jo>(data.frame(cbind(y2[
> >i,j,],y1[i,j,])),type="trace", K=2,
> >spec="transitory",ecdet="const",season=NULL,dumvar=NULL),error=function(e)
> >NaN)
> >}}
> >
> >and the earlier error does not show up. However, when I try to extract
> >some
> >results from the "co.jo" function, it does not work the way it does
> >without
> >the "tryCatch" thing and there appears another error. For example:
> >
> >for(i in 1:N1){
> >for(j in 1:N2){
> >co<-tryCatch(ca.jo <http://ca.jo>(data.frame(cbind(y2[
> >i,j,],y1[i,j,])),type="trace", K=2,
> >spec="transitory",ecdet="const",season=NULL,dumvar=NULL),error=function(e)
> >NaN)
> >
> >cor1<-cajorls(co,r=1)
> >}}
> >
> >Error in cajorls(vecm, r = 1) :
> >Please, provide object of class 'ca.jo' or 'cajo.test' as 'z'.
> >
> >So it seems "tryCatch" is disabling "ca.jo" from running properly. My
> >intention is to run this function ("ca.jo") and extract some estimates
> >for
> >further computation.
> >
> >regards
> >
> >Bukana
> >
> >
> >
> >
> >
> >
> >
> >On Mon, Dec 23, 2013 at 5:09 PM, Patrick Burns
> ><pburns@pburns.seanet.com>wrote:
> >
> >> There is a fundamental problem with your
> >> code, and there is the problem that you
> >> have (sort of) identified.
> >>
> >> The fundamental problem is that you are
> >> only going to get the results of the last
> >> call to 'ca.jo' that is done -- assuming
> >> it were to run.  You presumably want to
> >> save some information from each of the
> >> computations.
> >>
> >> You can get the loops to run even when you
> >> run into an error with some combinations
> >> by using 'try' or 'tryCatch'.  There is an
> >> example in Circle 8.3.13 of 'The R Inferno'.
> >>
> >> http://www.burns-stat.com/documents/books/the-r-inferno/
> >>
> >> If you have a question related to the actual
> >> function as opposed to general problems with
> >> R, then the r-sig-finance mailing list would
> >> be appropriate (you need to subscribe before
> >> posting).
> >>
> >> I'm not sure if this is enough of a hint for
> >> you or not.  If not, then trying to formulate
> >> a more explicit question might help.  (There
> >> are some suggestions in Circle 9 of 'The R
> >> Inferno'.)
> >>
> >> Pat
> >>
> >>
> >>
> >> On 23/12/2013 17:07, mamush bukana wrote:
> >>
> >>> Dear all,
> >>> I fit co-integration function between two integrated variables(y1
> >and y2)
> >>> over different grid points:
> >>>
> >>> for(i in 1:N1){
> >>> for(j in 1:N2){
> >>> co<-ca.jo(data.frame(cbind(y2[i,j,],y1[i,j,])),type="trace", K=2,
> >>> spec="transitory",ecdet="const",season=NULL,dumvar=NULL)
> >>> }}
> >>>
> >>> I have already extracted grid points with integrated time series.
> >However,
> >>> when I run the above function, there happens an error
> >>>
> >>> Error in solve.default(t(V) %*% SKK %*% V) :
> >>>    system is computationally singular: reciprocal condition number =
> >>> 1.10221e-35
> >>>
> >>> May you suggest me how to fix this problem please?
> >>>
> >>> Thanks in advance
> >>>
> >>>         [[alternative HTML version deleted]]
> >>>
> >>> ______________________________________________
> >>> R-help@r-project.org mailing list
> >>> https://stat.ethz.ch/mailman/listinfo/r-help
> >>> PLEASE do read the posting guide http://www.R-project.org/
> >>> posting-guide.html
> >>> and provide commented, minimal, self-contained, reproducible code.
> >>>
> >>>
> >> --
> >> Patrick Burns
> >> pburns@pburns.seanet.com
> >> twitter: @burnsstat @portfolioprobe
> >> http://www.portfolioprobe.com/blog
> >> http://www.burns-stat.com
> >> (home of:
> >>  'Impatient R'
> >>  'The R Inferno'
> >>  'Tao Te Programming')
> >>
> >
> >       [[alternative HTML version deleted]]
> >
> >______________________________________________
> >R-help@r-project.org mailing list
> >https://stat.ethz.ch/mailman/listinfo/r-help
> >PLEASE do read the posting guide
> >http://www.R-project.org/posting-guide.html
> >and provide commented, minimal, self-contained, reproducible code.
>
>

	[[alternative HTML version deleted]]

