search for: grave

...15:41, Avraham Adler wrote: > Floor (maybe round) of non-negative numerics, though. Poisson should > never have anything after decimal. > > Still think it?s worth allowing long long for R64 bit, just for purity > sake. > > Avi > > On Sun, Jan 19, 2020 at 4:38 PM Spencer Graves > <spencer.graves at prodsyse.com <mailto:spencer.graves at prodsyse.com>> wrote: > > > > On 2020-01-19 13:01, Avraham Adler wrote: >> Crazy thought, but being that a sum of Poissons is Poisson in the >> sum, can you break your ?big? simulation i...

...The development version on R-Forge uses this and returns ?Raul? under Windows and NA under Mac OS X (and presumably also Linux). The ?iconv? R code merely calls compiled code, which I?ve used very little in 30 years. Thanks, Spencer > On Nov 30, 2014, at 2:32 AM, Spencer Graves <spencer.graves at structuremonitoring.com <mailto:spencer.graves at structuremonitoring.com>> wrote: > > Wonderful. Thanks very much. Spencer > > > On 11/30/2014 2:25 AM, Milan Bouchet-Valat wrote: [[alternative HTML version deleted]]

....max.? On my Mac, .Machine$integer.max = 2147483647 = 2^31 > 1e9.? That still means that a Poisson distributed pseudo-random number just under that would have to be over 23000 standard deviations above the mean to exceed .Machine$integer.max. > > On Sun, Jan 19, 2020 at 1:58 PM Spencer Graves > <spencer.graves at prodsyse.com <mailto:spencer.graves at prodsyse.com>> wrote: > > ????? This issue arose for me in simulations to estimate > confidence, prediction, and tolerance intervals from glm(., > family=poisson) fits embedded in a BMA::bic.glm fit...

...t the repositories I installed lme4, maps, mapproj, CircStats, scatterplot3d, gregmisc, Hmisc without problems. To confirm, 'library(lme4)' produced the following error: Error: package 'Matrix' required by 'lme4' could not be found What do you suggest? Spencer Graves p.s. I get the same result using several different (US) mirrors. -- Spencer Graves, PhD Senior Development Engineer PDF Solutions, Inc. 333 West San Carlos Street Suite 700 San Jose, CA 95110, USA spencer.graves at pdf.com www.pdf.com <http://www.pdf.com> Tel: 408-938-4420 Fax: 408-280-7...

Hi, all, Does anybody got the most general expression of a unitary matrix? I found one in the book, four entries of the matrix are: (cos\theta) exp(j\alpha); -(sin\theta)exp(j(\alpha-\Omega)); (sin\theta)exp(j(\beta+\Omega)); (cos\theta) exp(j\beta); where "j" is for complex. However, since for any two unitary matrices, their product should also be a unitary matrix. When I

...> objects, and several of us are thinking about how me might get there. > But it isn't easy to get right, so it may take some time. I doubt > anything can happen for R 4.0.0 this year, but 2021 may be possible. > > I few notes inline below: > > On Sun, 19 Jan 2020, Spencer Graves wrote: > >> On my Mac: >> >> >> str(.Machine) >> ... >> $ integer.max????????? : int 2147483647 >> ?$ sizeof.long????????? : int 8 >> ?$ sizeof.longlong????? : int 8 >> ?$ sizeof.longdouble??? : int 16 >> ?$ sizeof.pointer?????? :...

????? This issue arose for me in simulations to estimate confidence, prediction, and tolerance intervals from glm(., family=poisson) fits embedded in a BMA::bic.glm fit using a simulate.bic.glm function I added to the development version of Ecfun, available at "https://github.com/sbgraves237/Ecfun".? This is part of a vignette I'm developing, available at "https://github.com/sbgraves237/Ecfun/blob/master/vignettes/time2nextNuclearWeaponState.Rmd". This includes a simulated mean of a mixture of Poissons that exceeds 2e22.? It doesn't seem unreasonable to m...

> On Dec 15, 2014, at 12:21 PM, Kurt Hornik <Kurt.Hornik at wu.ac.at> wrote: > >>>>>> Spencer Graves writes: > >> Hello, All: >> What would it take to make ?iconv? portable? > > >> I ask, because I want to convert accented characters to >> vanilla ASCII, thereby converting, e.g., ?Ra?l? to ?Raul?, and >> Milan Bouchet-Valet suggested on R-...

...;utils" "datasets" "methods" [7] "base" My work around is the following: pi2.2 <- (2*pi)^2 lines(c(0, 0), c(-pi2.2, pi2.2)) Why won't "lines(c(0,0), (2*c(-pi, pi))^2)" work for me? Thanks, Spencer Graves

Michael Grave just posted a question about surround conferences. http://www.facebook.com/notes.php?id=564633430#/note.php?note_id=5009726 3908&id=564633430&index=0 I didn't see it posted on the ast-list, what do you think? Does something like this have potential? I'd love to listen...

...3. The vignette mentioned below was extracted from a larger vignette fitting several models that seem to encounter convergence problems. I should probably switch to optimx using all the methods that offers for constrained optimization, including nminb. Best Wishes, Spencer Graves On 10/8/2016 5:00 PM, Mark Leeds wrote: > Hi Spencer: See the link below about L-BFGS-B below because I had > problems > with it a good while back (and I think the link description is the > cause but I can't prove it ) so eventually I moved to the Rvmmin(b) > package. I...

Dear all, I have a population with three age-classes, at time t=0 the population is: n.zero <- c(1,0,0) I have a transition matrix A which denotes "fertility" and "survival": A <- matrix(c(0,1,5, 0.3,0,0, 0,0.5,0), ncol=3, byrow=TRUE) To obtain the population at t=1, I calculate: A %*% n.zero To obtain the population t=2, I calculate: A %*% (A %*% n.zero) ... and so

Dear List, We are planning a genotyping study to be analyzed using false discovery rates (FDRs) (See Storey and Tibshirani PNAS 2003; 100:9440-5). I am interested in learning if there is any consensus as to how many features (ie. how many P values) need to be studied before reasonably reliable FDRs can be derived. Does anyone know of a citation where this is discussed? Bill Dupont William D.

again, under another subject: sorry, maybe an all too trivial question. But we have power data from J frequency spectra and to have the same range for the data of all our subjects, we just transformed them into % values, pseudo-code: power[i,j]=power[i,j]/sum(power[i,1:J]) of course, now we have a perfect linear relationship in our x design-matrix, since all power-values for each subject sum up

..."User guides, package vignettes and other documentation" -> "Ecdat::AverageIncomeModels". I've found other optimizers that will get around the problem in this case but none that performs as well as optim with many other problems. Thanks, Spencer Graves p.s. I've also tested bobyqa{minqa} or nloptr{nloptr}, recommended in a vignette in the lme4 package. These did better than optim in this example but worse in others I tried.

Dear R users; Ive got two questions concerning nlme library 3.1-65 (running on R 2.2.0 / Win XP Pro). The first one is related to augPred function. Ive been working with a nonlinear mixed model with no problems so far. However, when the parameters of the model are specified in terms of some other covariates, say treatment (i.e. phi1~trt1+trt2, etc) the augPred function give me the following

This error can arise when getOption("width") is too small. 80 seems to be the limit for me with R-3.3.2 on Windows. Bill Dunlap TIBCO Software wdunlap tibco.com On Wed, Mar 8, 2017 at 10:28 PM, Spencer Graves <spencer.graves at prodsyse.com> wrote: > Hello: > > > I tried "debug(help)" with the problem mentioned below. It stopped > with a call to "library", from which I generate the following simple > replication of this error: > > >> librar...

...arger integer range for R integer objects, and several of us are thinking about how me might get there. But it isn't easy to get right, so it may take some time. I doubt anything can happen for R 4.0.0 this year, but 2021 may be possible. I few notes inline below: On Sun, 19 Jan 2020, Spencer Graves wrote: > On my Mac: > > > str(.Machine) > ... > $ integer.max????????? : int 2147483647 > ?$ sizeof.long????????? : int 8 > ?$ sizeof.longlong????? : int 8 > ?$ sizeof.longdouble??? : int 16 > ?$ sizeof.pointer?????? : int 8 > > > ????? On a Windows 10 machi...

On Dec 15, 2014, at 1:37 PM, Spencer Graves <spencer.graves at prodsyse.com> wrote: > > >> On Dec 15, 2014, at 10:13 AM, Simon Urbanek <simon.urbanek at r-project.org> wrote: >> >>> >>> On Dec 15, 2014, at 12:21 PM, Kurt Hornik <Kurt.Hornik at wu.ac.at> wrote: >>> >>&...

> From: Spencer Graves [mailto:spencer.graves at pdf.com] > > Martin says, "This is another instance of S-Plus following R > behind and doing it incompatibly [with a reason?] ... ." > > This is one example of a major issue in "how to wage and win a > standards war",...