I have no idea. A data set that generates the error would be very helpful to me. What is
the role of the last line BTW, the one with "1%" on it?
Looking at the code I would guess that the vector "tied" has an NA in it, but how that
would happen I can't see. There is a reasonable chance that it will become obvious once I
have an example.
Terry Therneau
On 03/06/2015 05:00 AM, r-help-requ...@r-project.org wrote:
I don?t understand an error message from a thielsen function call within a
dplyr do function call.
by.CaseNo <- group_by(as.data.frame(MAP.dt), CaseNo)
>MAP.thielsen <- by.CaseNo %>%
+ do(model = thielsen(noninvMAP ~ invMAP, symmetric = T, data = .,
+ x = T, y = T, model = T, nboot = 1000))
| | 1% ~40 m remaining
Error in if (any(tied)) { : missing value where TRUE/FALSE needed
______________________________________________
R-help@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.