Re: [R] Number of Cores limited to two in CRAN

Dear Henrik, Thank you! This is quite helpful, especially your longer blog post. Best regards, Ravi From: Henrik Bengtsson Sent: Sunday, July 2, 2023 4:33 AM To: Ravi Varadhan Cc: R-Help Subject: Re: [R] Number of Cores limited to two in CRAN External

Re: [R] Number of Cores limited to two in CRAN

This is the specific error messsage from R CMD check --as-cran Error in .check_ncores(length(names)) : 16 simultaneous processes spawned Calls: prepost -> makeCluster -> makePSOCKcluster -> .check_ncores Execution halted Thanks, Ravi ____ From: Ravi

[R] Number of Cores limited to two in CRAN

, Ravi [[alternative HTML version deleted]] __ 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

Re: [R] How to use ifelse without invoking warnings

Thank you to Bert, Sarah, and John. I did consider suppressing warnings, but I felt that there must be a more principled approach. While John's solution is what I would prefer, I cannot help but wonder why `ifelse' was not constructed to avoid this behavior. Thanks & Best

[R] How to use ifelse without invoking warnings

rnings. Is there a better way to do this? Thank you, Ravi [[alternative HTML version deleted]] __ 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 ht

Re: [R] A Question about MLE in R

y do not see a good reason for not changing the default in optim. Best regards, Ravi [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

Re: [R] sqldf and number of records affected

V1 > 5 "))) sqldf("select * from main.con") sqldf() -Original Message- From: Gabor Grothendieck [mailto:ggrothendi...@gmail.com] Sent: Thursday, June 11, 2020 10:30 AM To: Ravi Jeyaraman Cc: r-help@r-project.org Subject: Re: [R] sqldf and number of records affected

Re: [R] sqldf and number of records affected

ldf() sqldf(c("pragma count_changes = 1", "update con set V1 = 0 where V1 > 5 ")) ans <- sqldf("select * from main.con") sqldf() -Original Message- From: Gabor Grothendieck [mailto:ggrothendi...@gmail.com] Sent: Thursday, June 11, 2020 9:12 AM To: Ravi Je

[R] sqldf and number of records affected

ss it's done on a database. Any ideas? Thanks Ravi -- This email has been checked for viruses by AVG. https://www.avg.com __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

Re: [R] how to add a calculated column into a data frame

How about something like this? df <- data.frame(ID=1:3, DTVAL=c("2009-03-21","2010-05-11","2020-05-05")) df <- df %>% mutate(YEAR = as.numeric(format(as.Date(DTVAL,'%Y-%m-%d'), '%Y'))) -Original Message- From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Charles Thuo Sent:

[R] Cumulative split of value in data frame column

Assuming, I have a data frame like this .. df <- data.frame(ID=1:3, FOO=c('A_B','A_B_C','A_B_C_D_E')) I want to do a 'cumulative split' of the values in column FOO based on the delimiter '_'. The end result should be like this .. ID FOO FOO_SPLIT1 FOO_SPLIT2 FOO_SPLIT

Re: [R] read_excel() ignore case of worksheet name?

Ista, With few tweaks this worked beautifully. Thank you so much. -Original Message- From: Ista Zahn [mailto:istaz...@gmail.com] Sent: Tuesday, May 26, 2020 11:38 PM To: Ravi Jeyaraman Cc: Erin Hodgess ; r-help@r-project.org Subject: Re: [R] read_excel() ignore case of worksheet name

Re: [R] read_excel() ignore case of worksheet name?

I’ve already tried that and doesn’t work From: Erin Hodgess [mailto:erinm.hodg...@gmail.com] Sent: Tuesday, May 26, 2020 10:55 PM To: Ravi Jeyaraman Cc: r-help@r-project.org Subject: Re: [R] read_excel() ignore case of worksheet name? Here’s a thought, please. Could you use the tolower

[R] read_excel() ignore case of worksheet name?

y thoughts? lapply(1:nrow(SIS), function(x) readxl::read_excel(SIS$FULL_FILEPATH[x], sheet='Tables', .name_repair = fixColNames)) Thanks in advance for your response. Cheers Ravi -- This email has been checked for viruses by AVG. https://www.avg.com ___

[R] Help with Parallel Processing

data to each child process? I am not doing that currently because I read somewhere that it's automatically done when using the multisession plan. Any idea what I'm doing wrong? Thanks Ravi if(!require('sqldf')) install.packages('sqldf') if(!require(

Re: [R] vectorizing the integrate function

find out how I can vectorize. I have obtained nice answers from Bert and I appreciate the help. I knew even previously that the vectorizing primarily improves the readability of the code. Thanks Bert for making it crystal clear that the looping is still at the interpreted level of R.Thanks,Ravi

Re: [R] vectorizing the integrate function

_along(a)){    res_i <- res[[i]]    for(j in seq_along(b)){      res_i_j <- res_i[[j]]      for(k in seq_along(m)){        r[i, j, k] <- res_i_j[[k]]$value        r.error[i, j, k] <- res_i_j[[k]]$abs.error      }    } } Hope this helps, Rui Barradas Às 21:14 de 10/08/19, ravi via R-help escreve

Re: [R] vectorizing the integrate function

function fv. That is, extract the fv$value as a vector. I can then perhaps try to improvise and extend it to arrays.Thanks,Ravi On Saturday, 10 August 2019, 20:03:22 CEST, Bert Gunter wrote: Ravi: First of all, you're calling Vectorize incorrectly. The first argument must be a fun

Re: [R] vectorizing the integrate function

)$value # How do I get this for a vector # Finally, if I want an array for all the values of a, b and m, how should I proceed? r <- array(0,c(3,6,3)) r.error <- array(0,c(3,6,3)) # I want the results of the vectorized integrate function in the above arrays. How do I extract these values? Tha

[R] vectorizing the integrate function

th dimensions of the lengths of a, b and c. I have tried several variants but am not having much luck. Will appreciate any help that I can get. Thanks,Ravi Sutradhara [[alternative HTML version deleted]] __ R-help@r-project.org mailing list --

Re: [R] Accessing C source files

lpful if you can explain what you mean by "downloading the source code".Thanks,Ravi On Monday, 5 August 2019, 15:32:00 CEST, Jeff Newmiller wrote: Seems like you are looking in the wrong place. Did you download the source code? [1] [1] https://cran.r-project.org/ On August

Re: [R] Accessing C source files

lpful if you can explain what you mean by "downloading the source code".Thanks,Ravi On Monday, 5 August 2019, 15:32:00 CEST, Jeff Newmiller wrote: Seems like you are looking in the wrong place. Did you download the source code? [1] [1] https://cran.r-project.org/ On August

[R] Accessing C source files

ry). I do not find any C file in the src subfolder  Where do I find the C source files? How do I get to look at them? I have a lot of follow-up questions on interfacing with the C files, but I would like to first know where I can find them. At least on how I can access them.Thank

Re: [R] plotting an isosurface on a 3d plot

ith the shading command results in a surface with a white patch in the middle. Can this be avoided?3. Just curious. When would one want to use homogenous coordinates in the definition of the vertices?4. If I want to paste separate images on the top and bottom planes, how would go about doing it

Re: [R] plotting an isosurface on a 3d plot

add an alpha value to the quad3d command? Once again, I would like to thank you for your fantastic help.Thanking you,Ravi On Thursday, 13 June 2019, 22:52:07 CEST, Duncan Murdoch wrote: On 13/06/2019 4:32 p.m., Duncan Murdoch wrote: > On 13/06/2019 12:47 p.m., ravi via R-help wr

[R] plotting an isosurface on a 3d plot

vectors. How can I plot this surface? In addition, I would like to add similiar surfaces to the same plot. I would appreciate any help that I can get.Thanking you,Ravi Sutradhara [[alternative HTML version deleted]] __ R-help@r-project.org mai

Re: [R] Fast matrix multiplication

Does Microsoft open R come with pre-compiled BLAS that is optimized for matrix computations? Thanks, Ravi -Original Message- From: Ista Zahn Sent: Monday, August 13, 2018 3:18 PM To: Ravi Varadhan Cc: r-help@r-project.org Subject: Re: [R] Fast matrix multiplication On Mon, Aug 13

Re: [R] Fast matrix multiplication

Hi Ista, Thank you for the response. I use Windows. Is there a pre-compiled version of openBLAS for windows that would make it easy for me to use it? Thanks, Ravi -Original Message- From: Ista Zahn Sent: Friday, August 10, 2018 12:20 PM To: Ravi Varadhan Cc: r-help@r-project.org

[R] Fast matrix multiplication

ng how to obtain more substantial speedup. Any suggestions would be greatly appreciated. Thanks, Ravi [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/l

Re: [R] adding overall constraint in optim()

%*% (cov.mat %*% wgt.vect)) LowerBounds<-c(0.2,0.05,0.1,0,0,0) UpperBounds<-c(0.6,0.3,0.6,0.15,0.1,0.2) spgSolution <- spg(wgt.vect, fn=opt.fun, lower=LowerBounds, upper=UpperBounds, project="projectLinear", projectArgs=list(A=matrix(1, 1, length(wgt.

[R] adding overall constraint in optim()

(p0)), b=1, meq=1)) Hope this is helpful, Ravi [[alternative HTML version deleted]] __ 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 ht

Re: [R] Subject: glm and stepAIC selects too many effects

More principled would be to use a lasso-type approach, which combines selection and estimation in one fell swoop! Ravi From: Ravi Varadhan Sent: Tuesday, June 6, 2017 10:16 AM To: r-help@r-project.org Subject: Subject: [R] glm and stepAIC selects too many

[R] Subject: glm and stepAIC selects too many effects

If AIC is giving you a model that is too large, then use BIC (log(n) as the penalty for adding a term in the model). This will yield a more parsimonious model. Now, if you ask me which is the better option, I have to refer you to the huge literature on model selection. Best, Ravi

Re: [R] Haplo.glm error: Failed to converge during EM-glm loop

verge when you set a small tolerance for convergence. You may fiddle with the control options for EM algorithm and see if you can get convergence. If not, you should contact the package authors. Hope this is helpful, Ravi Ravi Varadhan, Ph.D. (Biostatistics), Ph.D. (Environmental Engg) Associate Pr

[R] Linear Regressions with non-negativity constraint

Hi, Take a look at the package "ic.infer" by Ulrike Gromping. https://www.jstatsoft.org/article/view/v033i10 Best, Ravi Ravi Varadhan, Ph.D. (Biostatistics), Ph.D. (Environmental Engg) Associate Professor, Department of Oncology Division of Biostatistics & Bionformatics

Re: [R] help in maximum likelihood

Standard error = sqrt(diag(solve(opt$hessian))) Ravi From: Alaa Sindi [mailto:alaasi...@icloud.com] Sent: Wednesday, March 02, 2016 3:22 PM To: Ravi Varadhan Cc: r-help@r-project.org Subject: Re: help in maximum likelihood Thank you very much prof. Ravi, That was very helpful. Is there a way

Re: [R] help in maximum likelihood

6, 63, 59, 62, 60) # note: there is no need to have the choose(n, y) term in the likelihood fn <- function(p) sum( - (y*(p[1]+p[2]*x) - n*log(1+exp(p[1]+p[2]*x))) ) out <- nlm(fn, p = c(-50,20), hessian = TRUE) out eigen(out$hessian) Hope this is helpful, Ravi Ravi Varadhan, Ph.D.

[R] How to plot Geohash on a Map using R

0.00% ttnfs 0.04% ttnfy 0.01% Please let me know how this can be achieved. Thanks, Ravi [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

[R] Looking for a data set - Teratology data of Chen and Kodell (1989)

. Does anyone know of this or worked with this data set? If so, can someone share it with me? Thank you very much. Best, Ravi Ravi Varadhan, Ph.D. (Biostatistics), Ph.D. (Environmental Engg) Associate Professor, Department of Oncology Division of Biostatistics & Bionformatics Sidney Ki

Re: [R] Solve an equation including integral

te(fun, -Inf, t, nu=nu, exi=exi)$value return(res) } uniroot(myroot, c(-2, 2), nu=2, exi=0.5, alpha=.05) Ravi From: Ravi Varadhan Sent: Friday, January 08, 2016 11:29 AM To: r-help@r-project.org; 'sstol...@gmail.com' Subject: Re: Solve an equation including integral I think t

Re: [R] Solve an equation including integral

05) Hope this is helpful, Ravi Ravi Varadhan, Ph.D. (Biostatistics), Ph.D. (Environmental Engg) Associate Professor, Department of Oncology Division of Biostatistics & Bionformatics Sidney Kimmel Comprehensive Cancer Center Johns Hopkins University 550 N. Broadway, Suite -E Baltimore, MD

[R] grep/regexp

_percent_of_CD4" and "Fr_II_percent_of_CD4". How can I use grep() to do this? More generally, are there any good online resources with examples like this for the use of grep() and regexp() in R? I didn't find the help pages for these very user-friendly. Thank you very

Re: [R] Cautioning optim() users about "Nelder-Mead" default - (originally) Optim instability

mentations of Nelder-Mead. Best regards, Ravi From: ProfJCNash Sent: Sunday, November 15, 2015 12:21 PM To: Ravi Varadhan; 'r-help@r-project.org'; lorenzo.ise...@gmail.com Cc: b...@xs4all.nl; Gabor Grothendieck Subject: Re: Cautioning optim() user

[R] Cautioning optim() users about "Nelder-Mead" default - (originally) Optim instability

essage" in the optim help file that points this out to the users. Hope this is helpful, Ravi [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/li

Re: [R] About error: L-BFGS-B needs finite values of 'fn'

It seems like there is substantial finite-sample bias in the MLEs. Either that or there is some error in your procedure. See attached code. Ravi From: Ravi Varadhan Sent: Wednesday, November 11, 2015 2:33 PM To: 'denizozo...@gazi.edu.tr' ; r-help@r-project.org Cc: 'profj

Re: [R] About error: L-BFGS-B needs finite values of 'fn'

t can still happen). I would strongly suggest that when you are doing simulations, you should encapsulate the parameter estimation inside a `try' or `tryCatch' statement so that when there is an error, the simulation keeps going rather than crashing out. See the attached code. Best, Ravi

Re: [R] How to get variable name while doing series of regressions in an automated manner?

Thank you very much, Marc & Bert. Bert - I think you're correct. With Marc's solution, I am not able to get the response variable name in the call to lm(). But, your solution works well. Best regards, Ravi -Original Message- From: Bert Gunter [mailto:bgunter.4...@gm

[R] How to get variable name while doing series of regressions in an automated manner?

HS? In addition, is there a better paradigm for doing these type of series of regressions in an automatic fashion? Thank you very much, Ravi Ravi Varadhan, Ph.D. (Biostatistics), Ph.D. (Environmental Engg) Associate Professor, Department of Oncology Division of Biostatistics & Bi

Re: [R] Linear regression with a rounded response variable

Dear Peter, Charles, Gabor, Jim, Mark, Victor, Peter, and Harold, You have given me plenty of ammunition. Thank you very much for the useful answers. Gratefully, Ravi From: peter dalgaard Sent: Wednesday, October 21, 2015 8:11 PM To: Charles C. Berry Cc

[R] Linear regression with a rounded response variable

ered. Any suggestions? Thank you, Ravi Ravi Varadhan, Ph.D. (Biostatistics), Ph.D. (Environmental Engg) Associate Professor, Department of Oncology Division of Biostatistics & Bionformatics Sidney Kimmel Comprehensive Cancer Center Johns Hopkins University 550 N. Broadway, Suite -E Baltimore

Re: [R] Bug in auglag?

evaluates to `Inf' or some large value, when the parameter values are not in the constrained domain. constrOptim.nl() is a barrier method so it forces the initial value and the subsequent iterates to be feasible. Best, Ravi From: Rainer M Krug Sent: Tu

Re: [R] optimizing with non-linear constraints

I would recommend that you use auglag() rather than constrOptim.nl() in the package "alabama." It is a better algorithm, and it does not require feasible starting values. Best, Ravi -Original Message- From: Rainer M Krug [mailto:rai...@krugs.de] Sent: Thursday, October 01,

Re: [R] optimizing with non-linear constraints

nction(x, LAI) { h <- rep(NA, 2) h[1] <- LAI^x[2] / x[3] + x[1] h[2] <- 1 - x[1] - LAI^x[2] / x[3] h } Please take a look at the help page. If it is still not clear, you can contact me offline. Best, Ravi Ravi Varadhan, Ph.D. (Biostatistics), Ph.D. (Environmental Engg) Associate Profess

[R] what is the effective method to apply the below logic for ~1.2 million records in R

2 5 3 The above logic will take approximately 60 hours to generate the flag_1 column on a dataset consisting of ~1.2 million records. Is there any effective way in R to implement this logic in R ? Appreciate your help. Thanks, Ravi [[alternat

[R] lsqlin in pracma

not a deal breaker. Best, Ravi [[alternative HTML version deleted]] __ 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.o

Re: [R] plotting over a raster image with control over location and orientation

15)) wy<-(3-(-12))/(15-(-15)) xw<-convertWidth(unit(wx,"native"),"npc") yw<-convertHeight(unit(wy,"native"),"npc") vp2<-viewport(x = unit(x1, "npc"), y = unit(y1, "npc"),   width = unit(xw, "npc&quo

Re: [R] plotting over a raster image with control over location and orientation

[1] 1.02205882352941npc These converted values and the values for the limplt variable do not seem to be in agreement. How do I convert, for example, the point (10,-12) to "npc" units? I hope that I have explained my problem sufficiently well. If I can get help with an approach involving viewp

[R] plotting over a raster image with control over location and orientation

<-1:10 yy<-xx^2 plot(xx,yy,xlim=range(xx),ylim=range(yy),col="blue",type="b",xlab="x",ylab="square of x") # I would prefer to have the image more transparent just under the x and y labels and axis labelsThanks, Ravi [[alternative HTML version

Re: [R] : Ramanujan and the accuracy of floating point computations - using Rmpfr in R

What about numeric constants, like `163'? eval(Pre(exp(sqrt(163)*pi), 120))does not work. Thanks, Ravi From: David Winsemius Sent: Saturday, July 4, 2015 1:12 PM To: Duncan Murdoch Cc: r-help; John Nash; Ravi Varadhan Subject: Re: [R] : Rama

Re: [R] : Ramanujan and the accuracy of floating point computations - using Rmpfr in R

120), then all the constants, 163, pi, are automatically of 120 bits precision. Is this easy to do? Best, Ravi From: David Winsemius Sent: Friday, July 3, 2015 2:06 PM To: John Nash Cc: r-help; Ravi Varadhan Subject: Re: [R] : Ramanujan and the accuracy o

Re: [R] : Ramanujan and the accuracy of floating point computations - using Rmpfr in R

matters in Rmpfr. Even if it does, the answer you get is still wrong as you showed. Thanks, Ravi -Original Message- From: Richard M. Heiberger [mailto:r...@temple.edu] Sent: Thursday, July 02, 2015 6:30 PM To: Aditya Singh Cc: Ravi Varadhan; r-help Subject: Re: [R] : Ramanujan and the

[R] Ramanujan and the accuracy of floating point computations - using Rmpfr in R

ecision. Any thoughts as to what I am doing wrong? Thank you, Ravi [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read th

Re: [R] multiple parameter optimization with optim()

= "BFGS", hessian = TRUE) opt <- nlminb(startVal, fn, control=list(trace=1, rel.tol=1.e-06, iter.max=50)) # opt <- spg(startVal, fn) opt #list("coefficients" = opt$par, "LogLik" = -opt$val

Re: [R] Unable to use `eval(parse(text))' in nlme::lme

Yes, this is a very important point. Thank you, Bill. Best, Ravi From: William Dunlap Sent: Friday, February 13, 2015 4:56 PM To: Ravi Varadhan Cc: r-h...@stat.math.ethz.ch Subject: Re: [R] Unable to use `eval(parse(text))' in nlme::lme > ff <-

Re: [R] Unable to use `eval(parse(text))' in nlme::lme

) dd <- subset(labdata2, Transplant_type!=0 & time >0) lme(ff, random=~1|Patient, data=dd, correlation=corAR1(), na.action=na.omit) Best, Ravi Ravi Varadhan, Ph.D. (Biostatistics), Ph.D. (Environmental Engg) Associate Professor Department of Oncology Division of Biostatistics & Bionfor

[R] Unable to use `eval(parse(text))' in nlme::lme

yname') The same usage works well in lme4::lmer without any problems. It seems that there is a problem in how the formula object is evaluated in lme(). Is there an alternative way to do this? Thank you, Ravi [[alternative HTML version deleted]] _

Re: [R] Re-order levels of a categorical (factor) variable

Dear Bert, Martin, Bill, Thank you very much for the help. Bert & Martin's solutions solved my problem. Would it be useful to add an example like this to the help page on ?factor or ?levels? Thanks again, Ravi From: Martin Maechler Sent:

[R] Re-order levels of a categorical (factor) variable

used in whatever operation is done? I looked at relevel, reorder, etc., but they did not seem to be applicable to my problem. Thanks for any help. Best, Ravi [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To U

Re: [R] R 3.1.2 mle2() function on Windows 7 Error and multiple solutions

lems and it is quite difficult to know a priori which one would be the best for "my" problem. Thus, the benchmarking capability provided by optimx is a powerful feature. Ravi From: Ben Bolker Sent: Thursday, January 15, 2015 9:29 AM To: Rav

Re: [R] R 3.1.2 mle2() function on Windows 7 Error and multiple solutions

izers in optimx. If this solution is not what you are looking for, your problem may be poorly scaled. First, make sure that the likelihood is coded correctly. If it is correct, then you may need to improve the scaling of the problem. Hope this is helpful, Ravi [[alternativ

[R] Subject: short-sale constraint with nonpositive-definite matrix in portfolio optimization

of matrix projection is quite easy to do, and in fact, there are functions in R to do this (e.g., see posdefify() function in "sfsmisc" package). There is a recent review article on spectral projected gradient algorithm in J Statistical Software. http://www.jstatsoft.org/v60/i03 Hope th

Re: [R] porting an access database to sqlite

f starting from row 1, can I fetch rows from row number, 6e6+1?How do I know that I am near the tail end? There seems to be no simple method for finding the size of the table. Thanks,Ravi >>> ravi 11/25/14 12:42 PM >>> Hi,All my data is presently locked in a Microsoft access d

[R] porting an access database to sqlite

a detailed bit if code for doing this in a simple way. I would appreciate all help that I can get.Thanks,Ravi [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do rea

[R] R-tool - OS compatibility help

Dear Support, Please let us know Windows 2008R2 OS compatibility for "R" tool is available or not? If available, please let us know the details.

[R] Checking modeling assumptions in a binomial GLMM

binomial GLMM? Are there any standard checks that are commonly done? Are there any pedagogical examples or data sets where model assumptions have been examined for binomial GLMMs? Any suggestions/guidance is appreciated. Thank you, Ravi [[alternative HTML version deleted

Re: [R] Box-cox transformation

Thank you. It is very helpful. Ravi -Original Message- From: Joshua Wiley [mailto:jwiley.ps...@gmail.com] Sent: Monday, July 07, 2014 4:15 PM To: Ravi Varadhan Cc: r-help@r-project.org Subject: Re: [R] Box-cox transformation Dear Ravi, In my previous example, I used the residuals, so

[R] Box-cox transformation

="lambda", ylab="deviance") I am trying to create the profile likelihood for the Box-Cox parameter, but obviously I am not getting it right. I am not sure that ans$deviance is the right thing to do. I would appreciate any guidance. Thanks & Best, Ra

Re: [R] A combinatorial assignment problem

en maximum and minimum candidates for referees is 13 - 3 = 10. Of course, I have to increase the # iters to get a better solution, but for large K and R this may not converge at all. Best regards, Ravi From: Ravi Varadhan Sent: Wednesday, April 30, 2014 1:49 PM To: r-help@r-project.org Subject: A co

Re: [R] A combinatorial assignment problem

= N, R-Nrem, M*(K-i+1)) sel <- sample(1:R, size=size, replace=FALSE) end <- min((i+N-1),K) assignment[i:end, ] <- sel } assignment } sol <- assign(40,16,3) table(sol) Thanks, Ravi -Original Message- From: Bert Gunter [mailto:gunter.ber...@gene.com] Sent: T

[R] A combinatorial assignment problem

nment is not possible. I tried some obvious, greedy search methods but they are not guaranteed to work. Any hints or suggestions would be greatly appreciated. Best, Ravi [[alternative HTML version deleted]] __ R-help@r-project.org mailing list h

[R] A couple of issues with lattice::bwplot

type=rep(est.type, times=n.param*n.methods)) truth <- runif(16) # I would like to plot these points, one per each panel bwplot(method ~ Tx.effect|name, data=estimates, box.width=0, coef=0, pch=22) I would greatly appreciate any help. Thank you, Ravi [[alternativ

[R] Unable to Install a package from source in Windows

nks in advance, Ravi > install.packages("H:/Documents/computations/BB_2014.01-1.tar.gz", repos=NULL, > type="source") Installing package into '\\homer.win.ad.jhu.edu/users$/rvaradh1/Documents/R/win-library/3.0' (as 'lib' is unspecified) '\\homer.win.

Re: [R] convergence=0 in optim and nlminb is real?

will disappear. Best, Ravi [[alternative HTML version deleted]] __ R-help@r-project.org mailing list 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

Re: [R] [optim/bbmle] function returns NA

o notice that you do not use `tx' in your `g' function. There are likely a number of other issues as well. You may try the nmkb() function in the "dfoptim" package. It can handle box constraints and typically tends to perform a bit better than optim's Nelder-Mead for un

Re: [R] GARCH Optimization Problems

Hi Tony, Did you try the `auglag' algorithm in "alabma" package? Ravi [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting gui

[R] While using R CMD check: LaTex error: File `inconsolata.sty' not found

older. How can I eliminate this problem? Thanks, Ravi [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.htm

Re: [R] Vectorized code for generating the Kac (Clement) matrix

Thank you, Berend and Enrico, for looking into this. I did not think of Enrico's clever use of cbind() to form the subsetting indices. Best, Ravi From: Berend Hasselman [b...@xs4all.nl] Sent: Friday, April 26, 2013 10:08 AM To: Enrico Schumann Cc:

[R] Vectorized code for generating the Kac (Clement) matrix

n - i + 1 for (i in 2:(n+1)) Kacmat[i, i-1] <- i-1 The above code is fast, but I am curious about vectorized ways to do this. Thanks in advance. Best, Ravi Ravi Varadhan, Ph.D. Assistant Professor The Center on Aging and Health Division of Geriatric Medicine & Gerontology Johns Hopk

Re: [R] Crrstep help

Kathy, You need to make sure that you have installed the latest version of `cmprsk' package (version 2.2-4). crrstep will not work Otherwise. Ravi From: Ravi Varadhan Sent: Thursday, March 28, 2013 5:25 PM To: 'kathyhan...@gmail.com' Cc: r-help@r-project.org Subject: Re: [

Re: [R] Crrstep help

Hi Kathy, You should first contact the package maintainer (which is me!) before posting to r-help. Which version of the package are you using? Can you send me a minimally reproducible code? Ravi Ravi Varadhan, Ph.D. Assistant Professor The Center on Aging and Health Division of Geriatric

Re: [R] Constrained Optimization in R (alabama)

optimizers like alabama would still NOT be good approaches. It would be easiest to just go through the discrete values. Ravi From: Berend Hasselman [b...@xs4all.nl] Sent: Monday, February 11, 2013 12:45 AM To: Axel Urbiz Cc: R-help@r-project.org; Ravi

Re: [R] Constrained Optimization in R (alabama)

auglag(par = rep(1,4), fn = ff, heq = heq) Ravi From: Axel Urbiz [axel.ur...@gmail.com] Sent: Sunday, February 10, 2013 3:16 PM To: R-help@r-project.org; Ravi Varadhan Subject: Constrained Optimization in R (alabama) Dear List, I'm trying to solve this si

Re: [R] Constrained Optimization in R (alabama)

Your constraints are non-sensical. The only way for these constraints to be satisfied is for two of the parameters to be 0 and the other two to be 1. Please spend time to formulate your problem correctly, before imposing on others' time. Ravi From:

[R] Non linear programming: choose R that minimizes corr(y,x^R)

- function(p, x, y) sapply(p, function(p) cor(y, x^p)^2) plot(seq(0,5,length=1000), obj2(seq(0,5,length=1000),x,y), type="l") Best, Ravi Ravi Varadhan, Ph.D. Assistant Professor The Center on Aging and Health Division of Geriatric Medicine & Gerontology Johns Hopkins University

Re: [R] starting values in glm(..., family = binomial(link =log))

I did not get any warnings when I ran your data/model example. From: Fischer, Felix [mailto:felix.fisc...@charite.de] Sent: Wednesday, January 30, 2013 11:19 AM To: Ravi Varadhan Cc: r-help@r-project.org Subject: AW: [R] starting values in glm(..., family = binomial(link =log)) Thanks for your

[R] starting values in glm(..., family = binomial(link =log))

Try this: Age_log_model = glm(Arthrose ~ Alter, data=x, start=c(-1, 0), family=quasibinomial(link = log)) Ravi Ravi Varadhan, Ph.D. Assistant Professor The Center on Aging and Health Division of Geriatric Medicine & Gerontology Johns Hopkins University rvarad...@jhmi.edu<mailto

[R] texture mapping images on the faces of a cuboid with rgl

on the rwiki page referred to above). But, as can be seen from my questions, I don't understand this sufficiently to do this. But I hope that I have given sufficient information to explain my problem. Thanks, Ravi [[alternative HTML version deleted]] __

Re: [R] sqldf package: using variables in where condition

statement as list[1], list[2 }..tile list[15] > > Can someone help since, I am unable to reference the element location as > variable in the sqldf function? > > also, is there a better way? > would appreciate your help > > Rgds, Ravi > > > > -- >

[R] shell function in R

and is run from CMD console. How can I view the intermediate outputs while running NONMEM from R? With best regards, Ravi --- Ravi Singh, Ph.D. Post Doctoral Associate, Dept of Pharmaceutics Box 100494 JHMHC University of Florida Gainesville FL-32610 Ph: 352-273-7855 Fax: 35

Re: [R] model selection with spg and AIC (or, convert list to fitted model object)

, an important question is how do we know that the truth is supposedly on the boundary. If you can assume that the truth is in the interior, then you can use standard approaches: observed hessian or bootstrap to get standard errors. Best, Ravi From: Adam Zeilinger [mailto:ad...@ucr.edu] Sent

  1   2   3   4   5   6   7   8   9   >