search for: jpnolan

...quot;,i==.(i) )) title( str1 ) } #------------------------------------------------- ........................................................................... John P. Nolan Math/Stat Department 227 Gray Hall American University 4400 Massachusetts Avenue, NW Washington, DC 20016-8050 jpnolan at american.edu 202.885.3140 voice 202.885.3155 fax http://academic2.american.edu/~jpnolan ...........................................................................

...ot;)' pattern before without seeing this warning, so I wondered if the behavior had changed, and that's what prompted me to write. I know I can squelch the warning by wrapping the require() in suppressWarnings(). Dan ----- Original Message ----- > From: "John P. Nolan" <jpnolan at american.edu> > To: "Dan Tenenbaum" <dtenenba at fredhutch.org>, "R-devel" <r-devel at r-project.org> > Sent: Thursday, December 8, 2016 12:37:02 PM > Subject: RE: require(..., quietly=TRUE) does not suppress warning > Well, it says "most ofte...

Hi, The `quietly` argument of `require` is documented as follows: quietly: a logical. If ?TRUE?, no message confirming package attaching is printed, and most often, no errors/warnings are printed if package attaching fails. However: > require(foo, quietly=TRUE) Warning message: In library(package, lib.loc = lib.loc, character.only = TRUE, logical.return = TRUE, :

...on() output. I'd appreciate any suggestions. John ........................................................................................................... John P. Nolan Math/Stat Department 227 Gray Hall, American University 4400 Massachusetts Avenue, NW Washington, DC 20016-8050 jpnolan at american.edu voice: 202.885.3140 web: academic2.american.edu/~jpnolan ........................................................................................................... ----------------- Makevar.win ------------------------------------------ ifeq "$(WIN)" &qu...

...as.integer(method), as.raw(buf) , output=double(1) ) return(a$output) } John ........................................................................... John P. Nolan Math/Stat Department 227 Gray Hall American University 4400 Massachusetts Avenue, NW Washington, DC 20016-8050 jpnolan@american.edu 202.885.3140 voice 202.885.3155 fax http://academic2.american.edu/~jpnolan ........................................................................... [[alternative HTML version deleted]]

...be nice if there was a fast, general way to do this. John ............................................................................................ John P. Nolan Math/Stat Dept., American University Gray Hall, 4400 Massachusetts Ave, NW Washington, DC 20016-8050 Phone: 202-885-3140 E-mail: jpnolan at american.edu Web: http://fs2.american.edu/jpnolan/www/

...boundaries at another place: allowing gi==0 and interpreting log(gi) as -Inf. John ........................................................................... John P. Nolan Math/Stat Department 227 Gray Hall American University 4400 Massachusetts Avenue, NW Washington, DC 20016-8050 jpnolan at american.edu 202.885.3140 voice 202.885.3155 fax http://academic2.american.edu/~jpnolan ...........................................................................

...LSE,quietly=TRUE),silent=TRUE) and then check to see if aaa is of class "try-error" and check for failure John ??????????????????????????????.. John P. Nolan Math/Stat Dept., American University Gray Hall, 4400 Massachusetts Ave, NW Washington, DC 20016-8050 Phone: 202-885-3140 E-mail: jpnolan at american.edu Web: http://fs2.american.edu/jpnolan/www/ -----Original Message---- From: R-devel [mailto:r-devel-bounces at r-project.org] On Behalf Of Dan Tenenbaum Sent: Thursday, December 8, 2016 2:43 PM To: R-devel <r-devel at r-project.org> Subject: [Rd] require(..., quietly=TRUE)...

Hello,   Could someone help me please and to tell how to get the probability from empirical distribution (not parametric) for each data value (R function). For example, for normal distribution there is such a function like:   “pnorm(q, mean = 0, sd = 1, lower.tail = TRUE, log.p = FALSE)”   I need the same function only for the empirical distribution (which does not correspond to any typical

...s documented, or if there is an example somewhere, please let me know. Thanks, John ........................................................................... John P. Nolan Math/Stat Department 227 Gray Hall American University 4400 Massachusetts Avenue, NW Washington, DC 20016-8050 jpnolan@american.edu 202.885.3140 voice 202.885.3155 fax http://academic2.american.edu/~jpnolan ........................................................................... [[alternative HTML version deleted]]

Hello,   Could someone help me please and to tell how to get the probability from empirical distribution (not parametric) for each data value (R function). For example, for normal distribution there is such a function like:   “pnorm(q, mean = 0, sd = 1, lower.tail = TRUE, log.p = FALSE)”   I need the same function only for the empirical distribution (which does not correspond to any typical

...t;, x) } That I think it is a bit verbose but readable and could be made fast. And I think binding objects would be useful in other ways, as they are essentially a "named object". For example, when iterating over an environment. Michael On Thu, Jan 8, 2015 at 6:03 AM, John Nolan <jpnolan at american.edu> wrote: > Adding an optional argument to get (and mget) like > > val <- get(name, where, ..., value.if.not.found=NULL ) (*) > > would be useful for many. HOWEVER, it is possible that there could be > some confusion here: (*) can give a NULL because either...

Dear all, I am trying to use the rstable(n, alpha, beta, gamma = 1, delta = 0, pm = c(0, 1, 2))) function to generate positive stable random numbers. For positive stable distribution, beta==1 and alpha is in (0,1), which defines random variables with support (0, infinity). So, I used rstable(100, 0.5, 1) for an example. I found that this gives me some negative numbers. For example, >

Ben Bolker gives some reasons why Numerical Recipes may be problematic as a starting point for R codes. CUP did a masterful job of marketing, but the license is restrictive as the links he gives points out. In some tests, I've also noted that some of the algorithms are less than stellar e.g, convergence tests in one or two optimization routines. Should we have a wiki item to help people find

...nding on some decision made at run time, but 1. or 2. might be useful for the standard case. John .............................................................. John P. Nolan Math/Stat Department 227 Gray Hall, American University 4400 Massachusetts Avenue, NW Washington, DC 20016-8050 jpnolan at american.edu voice: 202.885.3140 web: academic2.american.edu/~jpnolan .............................................................. From: G?bor Cs?rdi <csardi.gabor at gmail.com> To: Gabriel Becker <gmbecker at ucdavis.edu>, Cc: "r-devel at r-project.org&...

In the help for ?integrate: >When integrating over infinite intervals do so explicitly, rather than just using a large number as the endpoint. This increases the chance of a correct answer ? any function whose integral over an infinite interval is finite must be near zero for most of that interval. I understand that and there are examples such as: ## a slowly-convergent integral integrand

Hello: The example "integrate(dnorm,0,20000)" says it "fails on many systems". I just got 0 from it, when I should have gotten either an error or something close to 0.5. I got this with R 2.12.0 under both Windows Vista_x64 and Linux (Fedora 13); see the results from Windows below. I thought you might want to know. Thanks for all your work in creating

Where does the value 2.2e-16 come from in p-values for chisq tests such as those reported below? > Anova(cm.mod2) Analysis of Deviance Table (Type II tests) Response: Freq LR Chisq Df Pr(>Chisq) B 11026.2 1 < 2.2e-16 *** W 7037.5 1 < 2.2e-16 *** Age 886.6 8 < 2.2e-16 *** B:W 3025.2 1 < 2.2e-16 *** B:Age 1130.4 8 < 2.2e-16 *** W:Age 332.9 8 < 2.2e-16 *** --- Signif.

On Wed, Apr 29, 2015 at 1:45 PM, Gabriel Becker <gmbecker at ucdavis.edu> wrote: > Gabor, > > To play devil's advocate a bit, why not just have the package formally > import the functions it wants to use (or the whole package if that is > easier)? > This is exactly my goal. And to facilitate this, R CMD check could remind you if you forgot to do the formal import.

...etIfExists, or .get, or get0, or whatever seems fine. Adding an argument to get() with missing giving current behavior may be OK too. Rewriting exists and get as .Primitives may be sufficient though. Best, luke > Michael > > > > > On Thu, Jan 8, 2015 at 6:03 AM, John Nolan <jpnolan at american.edu> wrote: > >> Adding an optional argument to get (and mget) like >> >> val <- get(name, where, ..., value.if.not.found=NULL ) (*) >> >> would be useful for many. HOWEVER, it is possible that there could be >> some confusion here: (*) ca...