search for: epidem

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 sugges...

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

.... are recreated which takes quite a while. Is there a way for R CMD check to compile only the help files that have changed since the last check (the way that R CMD INSTALL does)? Thanks, Frank -- Frank E Harrell Jr Prof. of Biostatistics & Statistics Div. of Biostatistics & Epidem. Dept. of Health Evaluation Sciences U. Virginia School of Medicine http://hesweb1.med.virginia.edu/biostat -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help&...

...e 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 -Frank -- Frank E Harrell Jr Prof. of Biostatistics & Statistics Div. of Biostatistics & Epidem. Dept. of Health Evaluation Sciences U. Virginia School of Medicine http://hesweb1.med.virginia.edu/biostat -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help&...

...(w$d) <- c('labelled',class(w$d)) class(w[['d']]) [1] "labelled" "POSIXt" "POSIXct" Is this a bug? Does anyone know of a workaround? Thanks. -- Frank E Harrell Jr Prof. of Biostatistics & Statistics Div. of Biostatistics & Epidem. Dept. of Health Evaluation Sciences U. Virginia School of Medicine http://hesweb1.med.virginia.edu/biostat -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help&...

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....

...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...

...e example inputs that each should result in the string or expression "x": x; x^2; g(x); h(g(x)); h(g(x^3)); j(h(g(pmin(x,3)))) h(g(x)+g(y)) To anyone who meets the challenge, thanks! -- Frank E Harrell Jr Prof. of Biostatistics & Statistics Div. of Biostatistics & Epidem. Dept. of Health Evaluation Sciences U. Virginia School of Medicine http://hesweb1.med.virginia.edu/biostat -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help&...

...year 2002 month 11 day 01 language R Package: foreign Priority: recommended Version: 0.5-10 Date: 2003-03-01 --- Frank E Harrell Jr Prof. of Biostatistics & Statistics Div. of Biostatistics & Epidem. Dept. of Health Evaluation Sciences U. Virginia School of Medicine http://hesweb1.med.virginia.edu/biostat

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....

...ply(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 Harrell Jr Prof. of Biostatistics & Statistics Div. of Biostatistics & Epidem. Dept. of Health Evaluation Sciences U. Virginia School of Medicine http://hesweb1.med.virginia.edu/biostat -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help&...

...using pdflatex rather than including postscript files with latex so that I can use the hyperref style for pdflatex. Has anyone had this problem or have a suggested solution? Thanks in advance. -- Frank E Harrell Jr Prof. of Biostatistics & Statistics Div. of Biostatistics & Epidem. Dept. of Health Evaluation Sciences U. Virginia School of Medicine http://hesweb1.med.virginia.edu/biostat -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help&...

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...

...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 f...

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