Hi, I have a problem with R, after an update: this piece of code: cat("creating resolver data frame\n"); dfG<-cbind(dfG,2^(RGN$G)) dfR<-cbind(dfR,2^(RGN$R)) suddenly, creates values to inf. Is the syntax changed for the ^ symbol in the latest R downloadable version? Kind regards, wilfred -- Dr. Ir. Wilfred Van IJcken Labmanager Genomics Erasmus Center for Biomics, Erasmus MC P.O. Box 1738, 3000 DR Rotterdam Dr. Molewaterplein 50, 3015 GE Rotterdam The Netherlands Tel +31 (0) 10 40 88 454 Fax +31 (0) 10 40 89 468 www.erasmusmc.nl/biomics <http://www.erasmusmc.nl/biomics>
What are RGN$G and RGN$R? Consider the following: > .Machine$double.xmax [1] 1.797693e+308 > log(.Machine$double.xmax, base=2) [1] 1024 > 2^(log(.Machine$double.xmax, base=2)-1) [1] 8.988466e+307 > 2^log(.Machine$double.xmax, base=2) [1] Inf If you want something different from this, PLEASE do read the posting guide! "http://www.R-project.org/posting-guide.html". It may help you formulate your question to increase the chances of getting a more useful reply. Best Wishes, spencer graves WFJ van IJcken wrote:> Hi, > > I have a problem with R, after an update: > this piece of code: > cat("creating resolver data frame\n"); > dfG<-cbind(dfG,2^(RGN$G)) > dfR<-cbind(dfR,2^(RGN$R)) > > suddenly, creates values to inf. > Is the syntax changed for the ^ symbol in the latest R downloadable version? > > Kind regards, wilfred-- 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-7915
Dear Sir/Madam I would be pleased if anybody can help me. I'm using linear mixed model (lme) function.I'm doing some simulation in my research and need to be assigned variance components values during of my program. Specifically, when I use lme function, I can get some information by use summary() and I can assign some valuse like variance of fixed parameters and variance of random error term by using for example varFix and sigma.But I don't know how I can assign for variance of random effect. I know in SPLUS we have command var.ran, how about R ? Thanks alot. M.Torabi