Hi all, If I use a function in R which takes some parameters as an input, how can I stop this function in the while loop and try another parameter in case the function takes long time or could not converge. Actually, I'm using "penalized" function in a loop for some fixed lambdas (pre-calculated), for some of them the function converged, but for some not. Now I decide to proceed in this way that if it takes longer than 3 mins, stop the "penalized" function and try the other lambdas. I need to do it in the loop since the loop is in the big program and I can't manually stop and begin again. Looking forward to your reply. NO [[alternative HTML version deleted]]
Have a look at evalWithTimeout() from the R.utils package
Best regards,
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature and
Forest
team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance
Kliniekstraat 25
1070 Anderlecht
Belgium
+ 32 2 525 02 51
+ 32 54 43 61 85
Thierry.Onkelinx at inbo.be
www.inbo.be
To call in the statistician after the experiment is done may be no more than
asking him to perform a post-mortem examination: he may be able to say what the
experiment died of.
~ Sir Ronald Aylmer Fisher
The plural of anecdote is not data.
~ Roger Brinner
The combination of some data and an aching desire for an answer does not ensure
that a reasonable answer can be extracted from a given body of data.
~ John Tukey
-----Oorspronkelijk bericht-----
Van: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org]
Namens n omranian
Verzonden: dinsdag 13 mei 2014 14:15
Aan: r-packages-owner at r-project.org; r-help at r-project.org; r-help-request
at r-project.org
Onderwerp: [R] stop a function
Hi all,
If I use a function in R which takes some parameters as an input, how can I stop
this function in the while loop and try another parameter in case the function
takes long time or could not converge.
Actually, I'm using "penalized" function in a loop for some fixed
lambdas (pre-calculated), for some of them the function converged, but for some
not. Now I decide to proceed in this way that if it takes longer than 3 mins,
stop the "penalized" function and try the other lambdas.
I need to do it in the loop since the loop is in the big program and I can't
manually stop and begin again.
Looking forward to your reply.
NO
[[alternative HTML version deleted]]
______________________________________________
R-help at 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.
* * * * * * * * * * * * * D I S C L A I M E R * * * * * * * * * * * * *
Dit bericht en eventuele bijlagen geven enkel de visie van de schrijver weer en
binden het INBO onder geen enkel beding, zolang dit bericht niet bevestigd is
door een geldig ondertekend document.
The views expressed in this message and any annex are purely those of the writer
and may not be regarded as stating an official position of INBO, as long as the
message is not confirmed by a duly signed document.
Hi All, I'm getting actually nuts. I don't understand the following lines in R: Here is the data, all rows are exactly the same!> ord_data[pid,]???????????? c0m2???? c0m4???? c0m8? c0m16??? c0m24??? c0m48???? c0p2???? c0p4???? c0p8?? c0p16??? c0p24??? c0p48?? c24m2??? c24m4??? c24m8? c24m16 13336382 5.632195 6.442133 5.818143 5.7683 5.862075 6.533181 5.807341 6.709774 5.664199 5.54022 5.385181 6.531977 5.29061 5.776121 6.136176 6.34699 13465320 5.632195 6.442133 5.818143 5.7683 5.862075 6.533181 5.807341 6.709774 5.664199 5.54022 5.385181 6.531977 5.29061 5.776121 6.136176 6.34699 13467455 5.632195 6.442133 5.818143 5.7683 5.862075 6.533181 5.807341 6.709774 5.664199 5.54022 5.385181 6.531977 5.29061 5.776121 6.136176 6.34699 13518680 5.632195 6.442133 5.818143 5.7683 5.862075 6.533181 5.807341 6.709774 5.664199 5.54022 5.385181 6.531977 5.29061 5.776121 6.136176 6.34699 ?????????? c24m24?? c24m48??? c24p2??? c24p4??? c24p8?? c24p16? c24p24?? c24p48 13336382 5.489161 6.043948 5.807802 5.761756 5.833559 5.293438 5.92068 5.759683 13465320 5.489161 6.043948 5.807802 5.761756 5.833559 5.293438 5.92068 5.759683 13467455 5.489161 6.043948 5.807802 5.761756 5.833559 5.293438 5.92068 5.759683 13518680 5.489161 6.043948 5.807802 5.761756 5.833559 5.293438 5.92068 5.759683 Now I found correlation:> pcor <- cor(t((ord_data[pid,]))) > pcor???????? 13336382 13465320 13467455 13518680 13336382??????? 1??????? 1??????? 1??????? 1 13465320??????? 1??????? 1??????? 1??????? 1 13467455??????? 1??????? 1??????? 1??????? 1 13518680??????? 1??????? 1??????? 1??????? 1 But, then I get this funny result !!! all(pcor[,1]==1) [1] FALSE> pcor[2,2]==1[1] TRUE> pcor[3,2]==1[1] FALSE Could anybody please comment on this?Many thanks. [[alternative HTML version deleted]]
Hi All, I'm getting actually nuts. I don't understand the following lines in R: Here is the data, all rows are exactly the same!> ord_data[pid,]???????????? c0m2???? c0m4???? c0m8? c0m16??? c0m24??? c0m48???? c0p2???? c0p4???? c0p8?? c0p16??? c0p24??? c0p48?? c24m2??? c24m4??? c24m8? c24m16 13336382 5.632195 6.442133 5.818143 5.7683 5.862075 6.533181 5.807341 6.709774 5.664199 5.54022 5.385181 6.531977 5.29061 5.776121 6.136176 6.34699 13465320 5.632195 6.442133 5.818143 5.7683 5.862075 6.533181 5.807341 6.709774 5.664199 5.54022 5.385181 6.531977 5.29061 5.776121 6.136176 6.34699 13467455 5.632195 6.442133 5.818143 5.7683 5.862075 6.533181 5.807341 6.709774 5.664199 5.54022 5.385181 6.531977 5.29061 5.776121 6.136176 6.34699 13518680 5.632195 6.442133 5.818143 5.7683 5.862075 6.533181 5.807341 6.709774 5.664199 5.54022 5.385181 6.531977 5.29061 5.776121 6.136176 6.34699 ?????????? c24m24?? c24m48??? c24p2??? c24p4??? c24p8?? c24p16? c24p24?? c24p48 13336382 5.489161 6.043948 5.807802 5.761756 5.833559 5.293438 5.92068 5.759683 13465320 5.489161 6.043948 5.807802 5.761756 5.833559 5.293438 5.92068 5.759683 13467455 5.489161 6.043948 5.807802 5.761756 5.833559 5.293438 5.92068 5.759683 13518680 5.489161 6.043948 5.807802 5.761756 5.833559 5.293438 5.92068 5.759683 Now I found correlation:> pcor <- cor(t((ord_data[pid,]))) > pcor???????? 13336382 13465320 13467455 13518680 13336382??????? 1??????? 1??????? 1??????? 1 13465320??????? 1??????? 1??????? 1??????? 1 13467455??????? 1??????? 1??????? 1??????? 1 13518680??????? 1??????? 1??????? 1??????? 1 But, then I get this funny result !!! all(pcor[,1]==1) [1] FALSE> pcor[2,2]==1[1] TRUE> pcor[3,2]==1[1] FALSE Could anybody please comment on this?Many thanks. [[alternative HTML version deleted]]
On Jan 30, 2015, at 2:20 AM, n omranian via R-help wrote:> Hi All, > I'm getting actually nuts. I don't understand the following lines in R: > Here is the data, all rows are exactly the same!Their print representation may be "exactly the same", but their internal values may differ by an amount that does not show up with your current print settings. We would need to see dput( ord_data[pid,] ) to comment further. -- David>> ord_data[pid,] > c0m2 c0m4 c0m8 c0m16 c0m24 c0m48 c0p2 c0p4 c0p8 c0p16 c0p24 c0p48 c24m2 c24m4 c24m8 c24m16 > 13336382 5.632195 6.442133 5.818143 5.7683 5.862075 6.533181 5.807341 6.709774 5.664199 5.54022 5.385181 6.531977 5.29061 5.776121 6.136176 6.34699 > 13465320 5.632195 6.442133 5.818143 5.7683 5.862075 6.533181 5.807341 6.709774 5.664199 5.54022 5.385181 6.531977 5.29061 5.776121 6.136176 6.34699 > 13467455 5.632195 6.442133 5.818143 5.7683 5.862075 6.533181 5.807341 6.709774 5.664199 5.54022 5.385181 6.531977 5.29061 5.776121 6.136176 6.34699 > 13518680 5.632195 6.442133 5.818143 5.7683 5.862075 6.533181 5.807341 6.709774 5.664199 5.54022 5.385181 6.531977 5.29061 5.776121 6.136176 6.34699 > c24m24 c24m48 c24p2 c24p4 c24p8 c24p16 c24p24 c24p48 > 13336382 5.489161 6.043948 5.807802 5.761756 5.833559 5.293438 5.92068 5.759683 > 13465320 5.489161 6.043948 5.807802 5.761756 5.833559 5.293438 5.92068 5.759683 > 13467455 5.489161 6.043948 5.807802 5.761756 5.833559 5.293438 5.92068 5.759683 > 13518680 5.489161 6.043948 5.807802 5.761756 5.833559 5.293438 5.92068 5.759683 > > Now I found correlation: >> pcor <- cor(t((ord_data[pid,]))) >> pcor > 13336382 13465320 13467455 13518680 > 13336382 1 1 1 1 > 13465320 1 1 1 1 > 13467455 1 1 1 1 > 13518680 1 1 1 1 > But, then I get this funny result !!! > all(pcor[,1]==1) > [1] FALSE >> pcor[2,2]==1 > [1] TRUE >> pcor[3,2]==1 > [1] FALSE > Could anybody please comment on this?Many thanks. > > > > > > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see > 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.David Winsemius Alameda, CA, USA