search for: epidemic

Hi all, following is my R -code and shows the error given below > n <- 100 > k<-2 > x1 <-c(1, 3);x2 <- c(2,5) > X <- matrix(c(0,0), nrow = 2, ncol = n) > for(i in 1:k) + X[i, ] <- mh1.epidemic(n,x1[i],x2[i]) Error in X[i, ] <- mh1.epidemic(n,x1[i],x2[i]):   number of items to replace is not a multiple of replacement length   > psi <-t(apply(X, c(1,2), cumsum)) > for(i in 1:nrow(psi)){ + psi[i,] <psi[i,]/(1:ncol(psi)) + print(Gelman.Rubin(psi)) + + } Can some one suggest...

Hi R community, I have some data set and construct the likelihood as follows likelihood <- function(alpha,beta){ lh<-1 d<-0 p<-0 k<-NULL data<-read.table("epidemic.txt",header = TRUE) attach(data, warn.conflicts = F) k <-which(inftime==1) d <- (sqrt((x-x[k])^2+(y-y[k])^2))^(-beta) p<-1 - exp(-alpha*d) for(i in 1:100){ if(i!=k){ if(inftime[i]==0){ lh<-lh*(1-p[i]) }...

Given the recent spam epidemic, I would like to hear people's opinion on switching this list to subscribed-post only (all non-subscribers would have their posts held for moderation.) -hpa

Hi i am doing an environmental research The equation is as follow: gam(y1 ~ x1 + s(x2) + s(x3) + s(x4), family = gaussian, fit = true) I would like to obtain the beta coefficient and 95CI of x4 (or s(x4)), what should I do? Thanks, Lung -- View this message in context: http://r.789695.n4.nabble.com/mgcv-beta-coefficient-and-95-CI-tp3320491p3320491.html Sent from the R help mailing list

I am unclear on whether to run R CMD check pgkname as user or as root on Linux. When running as user, after all the latex, html, and man files are created, I get the following error message: Rdconv(): Couldn't open '': Permission denied Has anyone dealt with that message? When I run R CMD check a second time, all latex, html, etc. are recreated which takes quite a while. Is there

In a .First.lib I want to issue two require()s to insure that two other packages are loaded. But I want the package being loaded by .First.lib using library.dynam("mypackage",pkb,lib) to be higher in the search order than the two required packages, because I want to have a couple of functions from the two required packages overridden. What is the best way to do that? Thanks in advance

Under platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 1 minor 5.1 year 2002 month 06 day 17 language R I am having the following problem. d <- Sys.time() w <- data.frame(d)

Following is my code, can some one help on the error at the bottom? > mh<-function(iterations,alpha,beta){ + data<-read.table("epidemic.txt",header = TRUE) + attach(data, warn.conflicts = F) + k<-97 + d <- (sqrt((x-x[k])^2 + (y-y[k])^2)) + p <- 1-exp(-alpha*d^(-beta)) + p.alpha<-1 - exp(-3*d^(-beta)) + p.beta <- 1 - exp(alpha*d^(-2)) + iterations<-1000 + mu.lambda <- c(0,0);s.la...

...ation] ]~ ******* ------>>> Which is more important, $$$ or your daughters [Question] <<<------ ******* ~~~RED ALERT~~~ {~See CDC quote---76% 'skyrocket' rise is real danger all across The Planet, [internet incited] if 30,000 to 50,000 girls commit suicide in a 2 day "epidemic", like the 1986 youth suicide epidemic in Omaha, NE that SUICIDE VACCINE stopped~} * {Another equal danger "The most dangerous thing in the world is a young, unemployed man!" ~Taxi driver in Nairobi, Kenya~[during 2007-2008] {Is America listening?!? signed....CaptainChurch} * This is...

I am trying to replicate part of the function of the S-Plus terms.inner function. I don't need an entire model terms object but R expressions or character strings containing the innermost variable name given a vector of character strings or expressions. Here are some example inputs that each should result in the string or expression "x": x; x^2; g(x); h(g(x)); h(g(x^3));

I am using read.spss in the foreign package to read an SPSS save file. For date variables I get huge values such as 11489990400. Does anyone know how to convert these values to R POSIXct date objects? Thanks in advance -Frank platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 1 minor

Hello R-user I want to compute a multiple regression but I would to include a check for collinearity of the variables. Therefore I would like to use a ridge regression. I tried lm.ridge() but I don't know yet how to get p-values (single Pr() and p of the whole model) out of this model. Can anybody tell me how to get a similar output like the summary(lm(...)) output? Or if there is

Hi all I have a function of log defined by y = log(1- exp(-a)), when exp(-a) is greater, 1, it produce NaN. How can I remove this in R? [[alternative HTML version deleted]]

Hi all, Our company has deployed a Mitel 3300 system (only about 2.5 years ago) and we are experimenting with setting up Asterisk in our head office (for business continuity, ie we have a bird flu epidemic and no-one can come in, therefore use SIP softphones at home to co-ordinate activity) and at a remote site in the Isle of Man (connected via 2Mbps SDSL) Ideally we'd like anyone on either Asterisk servers (IOM and London) to be able to dial anyone internally on the Mitel 3300 and vice-versa. W...

Does anyone have something like tapply that is extremely fast for matrices when there is a very large number of levels of the grouping variable? I'm referring to, for example, tapply(x, grouping.variable, function.operating.on.submatrix) where x is a matrix and the submatrix is a subset of the rows of x. The grouping variable's length equals the number of rows of x. -- Frank E

I am having good success in using pdf() to produce pdf graphics files for inclusion in documents processed by pdflatex [I'm using R 1.6.0 on Linux]. I am having one small problem consistently though, which I believe I also had with S-Plus's pdf.graph(): wasted space at the top of the graph which causes pdflatex to put the graph on a later page, when the current page has plenty of space.

I am using GLM to calculate logit models based on cross-sectional data. I am now down to the hard work of making the results intelligible to very average readers. Is there any way to calculate a psuedo analoque to the R^2 in standard linear regression for use as a purely descriptive statistic of goodness of fit? Most of the readers of my report will be vaguely familiar and more comfortable with

...and many other companies: http://blog.revolutionanalytics.com/2017/06/applications-earl-sf-2017.html A demo of real-time predictions from a model created with Microsoft R, at a rate of one million predictions per second: http://blog.revolutionanalytics.com/2017/06/real-time-predictions.html The R Epidemics Consortium is a coalition of researchers developing epidemiology resources for R: http://blog.revolutionanalytics.com/2017/06/r-epidemics-consortium.html Syberia is a on open-source framework for orchestrating R scripts in production: http://blog.revolutionanalytics.com/2017/06/syberia.html A fr...

...d that every effort to contain the virus and slow the spread saves lives. As EuroLLVM attracts attendees from all over the world, we feel that cancellation is an effort to contain and stop the virus. 1. We believe France will reach stage 3, which is the last stage of the Covid-19 epidemic in the next few days. Once stage 3 is deployed, the French government could implement new restrictions that impact those who attempt to travel to EuroLLVM 2020 or impact large gatherings such as EuroLLVM. In addition, several countries have imposed a mandatory quarantine upon return for...

Hi all, I'm currently working with a dataset that has quite a few missing values and after some investigation I figured that multiple imputation is probably the best solution to handle the missing data in my case. I found several references to functions in S-Plus that perform multiple imputation (NORM, CAT, MIX, PAN). Does R have corresponding functions? I searched the archives but was not