search for: rhelp02a

something like "sprintf" in C? so I can do: print(sprintf("the correct result is %3.4f\n", myresult)); ------- Also, I am desperately looking for a "clear console screen" function in R... thanks a lot! [[alternative HTML version deleted]]

Is there a function in R for doing Greenhouse-Geiser correction in ANOVA models? Is it already available in the aov function? How do we use it? Best, Darren [[alternative HTML version deleted]]

Hello. Some months age an e-mail was posted in which a comparison between Glimmix and glm was discussed. I have not been able to find that e-mail on the R archive. Does anyone recall the date of the above e-mail? Thank you very much. ******************************************* Antonio Paredes USDA- Center for Veterinary Biologics Biometrics Unit 510 South 17th Street, Suite 104 Ames, IA 50010

...R 2.4.0, this is what happens: > gregexpr(" [a-z] [a-z] ", " a b c d e f ", perl=T) [[1]] [1] 1 7 attr(,"match.length") [1] 5 5 Looking at the archives, I came across these sites where the reverse issue has been discussed before: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/75843.html http://finzi.psych.upenn.edu/R/Rhelp02a/archive/76815.html http://finzi.psych.upenn.edu/R/Rhelp02a/archive/75846.html >From there, it seems as if the first result has been considered undesirable (apparently because it differs from Perl's output if not also for other reaso...

...ail before 2002 into one index - without affecting existing links). This function runs from within R and opens a web page with the results of the search. For example, RSiteSearch("multiple imputation") Jon -- RSiteSearch <- function(string, restrict=c("Rhelp02a","Rhelp01","Rhelp00","functions","docs"), format="normal", sortby="score", matchesPerPage=10) { string <- paste("http://finzi.psych.upenn.edu/cgi-bin/namazu.cgi?query=", gsu...

Hello- It has been several years since anyone has asked, so i am asking again- has anyone created a routine to estimate confidence intervals for predictions from nls models (ala Bates and Watts 1988)? Thanks - Alice Shelly [[alternative HTML version deleted]]

...ndexed by the search engine (Namazu). As a result, search should be a bit less slow and the lower items in the list should be less irrelevant. The down side of all this is that specific reference to other messages has been broken for the last two days. Please try not to refer to other message as Rhelp02a.new or tmp/Rhelp02a. They should be just Rhelp02a, and they will be. To fix this, I need to re-build the mailing list index. I plan to do this about 01:00 UTC (2 AM in London, 9 PM in Philadelphia) tomorrow night (June 23/24). Things may not work for a couple of hours around that time. The old...

The CRAN host in Tampa, FL (cran.hostingzero.net) isn't responding, and hasn't responded in quite a while -- at least problems were reported more than a year ago (Oct 2007) http://finzi.psych.upenn.edu/R/Rhelp02a/archive/116706.html although someone apparently succeeded in July 2008 http://finzi.psych.upenn.edu/R/Rhelp02a/archive/137451.html Perhaps it's just flaky, and not permanently missing, but it still causes problems (especially for new users who have no idea what's wrong, they just selec...

Hello: Could anyone point me to a nice example where someone has created methods for "[" on a user defined Class? I looked at the package Matrix but that was a little daunting. I'm looking for someone a little more introductory. I've tried to search the help section and the web but its difficult since "[" isn't searchable. Thanks in advance! Greg [[alternative

...lvable with uniroot to find "a": fn=function(a){ b=(0.7/a)-a (1/(a+b+1))-0.0025 } uniroot(fn,c(-500,500)) gives "Error in uniroot(fn, c(-500, 500)) : f() values at end points not of opposite sign" I read R-help posts and someone wrote a function: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/92407.html but it is not very precise. Is there any '"standard" function in R that can solve this? thanks. -- View this message in context: http://www.nabble.com/How-to-solve-difficult-equations--tf3643595.html#a10175603 Sent from the R help mailing list archive at Nabble.com...

...t SAS xport files that contain multiple files using package foreign's read.xport. I first attempted this back in 2005 and had problems. Some of files that were present in the SAS xport file weren't being created in R. I submitted my problem to the community: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/57864.html and Dr. Dalgaard confirmed the issue: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/57868.html >From what I've read, this problem was identified before 2005: http://tolstoy.newcastle.edu.au/R/help/03a/0527.html After updating to version 2.6.1 recently, I decided to tr...

...RSiteSearch("lme prediction intervals") produced 3 hits including 2 from James A Rogers over 3 years ago. In one, he said, "I am not aware of any published R function that gives you prediction intervals or tolerance intervals for lme models." (http://finzi.psych.upenn.edu/R/Rhelp02a/archive/42781.html) In the other, he provided sample code for prediction or tolerance intervals of lme variance components. (http://finzi.psych.upenn.edu/R/Rhelp02a/archive/44675.html)

...I get the following warning 'gamma' cannot be modified on this device in: par(args) with 2 empty devices. I do not understand why I get those warnings and why I get those empty devices. I looked up in RSiteSearch and I found the 2 following posts http://finzi.psych.upenn.edu/R/Rhelp02a/archive/38553.html In this post M. Schwarz explains that there is not a plot created that's why R generates the warning. In the case of tstpar, I generate the hist plots and should not be getting the warning. And in the following post Prof Ripley refers to an old R warning for pdf() devices W...

Hi to all, I found in the R-help archive how to calculate the p-value for a gee result: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/74150.html but there are two questions (I am afraid they are basic questions ...) 1. why is the result multiplicated with 2 2. how could I decide between lower.tail =TRUE and FALSE: example: 2*pnorm(c(1.8691945,0.5882351,2.4903091,1.9287802,2.3172983,2.2092593,2.2625959,1.6395695), lower...

Suppose I have the following data: y x id 44 0 104 48 58 104 48 55 204 47 105 204 41 275 206 18 67 209 ....... I fit the model >fit=lme(y~x+I(x^2),random=~1|id) Now I want to make a prediction plot: >time=seq(0,300,len=100) >plot(predict(fit,data.frame(x=time),level=0)) Very fine. It gives me the prediction curve based on the model. My further request is to make a confidence bands

hie all i am trying to carry out a categorical data analysis but my problem is that when in i use the chi squared test some of my expected values are less than 5. is there a test that can handle this situation. the data is not a 2*2 table. its more from the social sciences where you have from strongly agree to strongly disagree. i know i can collapse vthe tables but there is a loss of

Hello R-Group I found how to fill the data.frame -> http://finzi.psych.upenn.edu/R/Rhelp02a/archive/70843.html N1 <- rnorm(4) N2 <- rnorm(4) N3 <- rnorm(4) N4 <- rnorm(4) X1 <- LETTERS[1:4] ################### nams <- c(paste("N", 1:4, sep = ""), "X1") dat <- data.frame(lapply(nams, get)) names(dat) <- nams dat But I need also to cr...

I would like to add points to a wireframe but with a conditioning variable. I found a solution for this without a conditioning variable here, http://finzi.psych.upenn.edu/R/Rhelp02a/archive/65321.html. Does anyone know how to plot a wireframe conditioned on a variable and add the points conditioned on the same variable, similar to the solution at the link above?

...sep="") res <- charmatch(restrict, c("Rhelp", "doc", "function", "doc/fun")) if (res == 0) stop("wrong restriction specified") res <- switch(res, "Rhelp00/archive|Rhelp01/archive|Rhelp02a/archive", "finzi.psych.upenn.edu/R/doc", "finzi.psych.upenn.edu/R/library", "finzi.psych.upenn.edu/R/doc|finzi.psych.upenn.edu/R/library") res <- paste("restrict=", res, sep="") qstring <- pas...

Hi to all, maybe the last question was not clear enough. I did not found any hints how to decide whether it should use lower.tail or not. As it is an extra R-feature ( written in http://finzi.psych.upenn.edu/R/Rhelp02a/archive/66250.html ) I do not find anything about it in any statistical books of me. Regards Carmen