search for: rhelp02

HI, I am using lmer() for a simple mixed effects model. The model is of the form logit(y)~ x + (1|z), where x is an indicator variable and z a multi-level factor. I would like an estimate of the response variable (either y or logit y) with an associated confidence interval for a given value of x. There does not appear to be a predict function written for lmer(). The output for the fixed

In trying to run John Fox's new Rcmdr package on Redhat Linux 9, I got this error message. The problem is in loading the tclck library. library(tcltk) produces the same error. It has been discussed before, e.g., http://finzi.psych.upenn.edu/R/Rhelp02/archive/11898.html but with respect to earlier versions of Linux and R. I tried updating to the latest RPMs of tcl and tk (8.3.5-89), but that didn't help. Otherwise I just have the standard RH9 installation. Has anyone gotten tcltk to work on RH9? Jon

Hi everyone, I'm building a matrix in R with a cycle for like this: pp_ris2=matrix(NA,6,6) for(i in 0:6){ for(j in 0:6){ if(i>j){ pp_ris2[i,j]=myfunction} else if(i==j){ print(c(i,j)) pp_ris2[i,j]=myfunction} }} but the result is: [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0.062 NaN NaN NaN NaN NaN [2,] 0.083 0.022 NaN NaN NaN NaN [3,]

...that some of the earliest threads here refer to the problem (now fixed) that the subset() documentation failed to note that the existing 'drop' argument would *not* (confusingly) drop unused levels. http://finzi.psych.upenn.edu/Rhelp10/2008-April/158566.html http://finzi.psych.upenn.edu/R/Rhelp02/archive/42976.html http://finzi.psych.upenn.edu/R/Rhelp02/archive/36961.html http://finzi.psych.upenn.edu/Rhelp10/2009-November/217878.html http://article.gmane.org/gmane.comp.lang.r.general/200395 This suggestion is milder and less wide-ranging than a global drop.unused.levels option, or than c...

Dear all, I am looking at the R-help entry below: http://finzi.psych.upenn.edu/R/Rhelp02/archive/26640.html I have a more complicatedt problem. I have a zoo time series frame with 100+ sequences. I want to cycle through them back and forth and compare them to the 1st column at any time. I need also a button to click when I need the viewed-selected sequence (that is being compared...

...ipt via an external file, this would be a tedious business! Is there some function already in R which outputs the bits in the binary representation directly? I see that Dabid Hinds asked a similar question on 17 Aug 2005: "Raw data type transformations" http://finzi.psych.upenn.edu/R/Rhelp02/archive/59900.html (without, apparently, getting any response -- at any rate within the following 3 months). With thanks for any suggestions, Ted. -------------------------------------------------------------------- E-Mail: (Ted Harding) <Ted.Harding at manchester.ac.uk> Fax-to-email: +44...

A common point made in discussion of contrasts, type I, II, III SS etc is that for sensible comparisons one should use contrasts that are 'orthogonal in the row-basis of the model matrix' (to quote from http://finzi.psych.upenn.edu/R/Rhelp02/archive/111550.html) Question: How would one check, in R, that this is so for a particular fitted linear model object? Steve Ellison ******************************************************************* This email and any attachments are confidential. Any use...{{dropped:8}}

Hi everyone, I have a very quick question: Is there a ready-made function in R or any R packages to find the prime factorization of an integer? -- View this message in context: http://r.789695.n4.nabble.com/Prime-Factorization-tp2548877p2548877.html Sent from the R help mailing list archive at Nabble.com.

Does anyone know which package include the computation of Cochran’s Q statistic in R? jlfmssm [[alternative HTML version deleted]]

Hi, I am aware that this question might be a basic one. I did browse the help and archives but I still haven't understood how to do the following: I run R locally and would like to read in data from a server for which there is a username and password. That is, how do I open a connection to a server with a password? Thanks, Hadassa -- Hadassa Brunschwig PhD Student Department of Statistics

Hi, How can I generate data from multivariate gamma distribution & multivariate beta distribution? I only found command for multivariate normal only. Many thanks in advance :) Regards Rofizah [[alternative HTML version deleted]]

hi,R users, Now I have some scatter figures, is there some method can plot equalprobable error ellipse ? -- TANG Jie Email: totangjie@gmail.com Tel: 0086-2154896104 Shanghai Typhoon Institute,China [[alternative HTML version deleted]]

Dear All, Can you help me, with the code below how do I obtain the fisher information from it. Is my q<-replicate(1000,x) the right way to do simulation. thank you. x<-rweibull(100,0.8,1.5) q<-replicate(1000,x) z<-function(p){ beta<-p[1] eta<-p[2] log1<-(n*log(beta)-n*beta*log(eta)+(beta-1)*sum(log(x))-sum((x/eta)^beta)) return(-log1) } zz<-optim(c(0.5,0.5),z) zz Chris

Dear Group: I want to use lattice with a formula such as y ~ x | v to plot a data frame in which v varies to indicate which "x" is really being plotted. I know how to make the x-axis scales vary with the panel but is it possible to let the x-axis label vary, i.e., to let the user specify a vector of x-axis labels? Thanks Frank -- Frank E Harrell Jr Professor and Chair

Dear all! I am using wireframe function from lattice package. Is it possible to remove box around the plot but to keep axis (x, y, z)? Thanks in advance! [[alternative HTML version deleted]]

...mine encountered some unexpected behavior regarding the postscript output from R. It's difficult for me to tell whether or not this is an R problem or a ghostview/gv/interpreter problem. Just to note, I think it's exactly the same situation reported here: http://finzi.psych.upenn.edu/R/Rhelp02/archive/25436.html The following code produces a working plot (no problems, as far as I can see): par(mfrow = c(2,2)) plot.new() plot(1:10, 1:10, xlab = "", ylab = "") plot(1:10, 1:10, xlab = "", ylab = "") plot(1:10, 1:10, xlab = "", ylab = &quo...

...rDevices utils datasets splines methods [8] base other attached packages: [1] tractor.session_1.0.0 tractor.base_1.0.3 tractor.nt_1.0.2 loaded via a namespace (and not attached): [1] tools_2.8.1 Regards, Jon (1) http://code.google.com/p/tractor/ (2) http://finzi.psych.upenn.edu/R/Rhelp02/archive/61872.html -- Jonathan D. Clayden, Ph.D. Research Fellow Radiology and Physics Unit UCL Institute of Child Health 30 Guilford Street LONDON WC1N 1EH United Kingdom t | +44 (0)20 7905 2708 f | +44 (0)20 7905 2358 w | www.homepages.ucl.ac.uk/~sejjjd2/

Hello UseRs, I need help on how to plot several raster images (such as those obtained from a kernel-smoothed intensity function) in a layout such as that obtained from the lattice package. I would like to obtain something such as obtained from using the "levelplot" or "xyplot" in lattice. I currently use: >par(mfrow=c(3,3) to set the workspace, but the resulting plots

...size = 2),), rectangles = list(col = c(2,3,4,5,6), size = 1, border = FALSE)) ) I could get things to work when the legend was on top using http://finzi.psych.upenn.edu/R/Rhelp02/archive/46654.html on something similar but could not get it to work in the last panel of a multipanel plot. Regards Duncan Mackay Department of Agronomy and Soil Science University of New England Armidale NSW 2351 Email: home mackay at northnet.com.au R version 2.12.1 (2010-12-16) Platform: i3...

Dear R-helpers, sorry if this is obvious, but I can't find this in the documentation. I am using Sweave, and have some code that does not actually work - but I want to include it anyway, including the error message that R produces. But on running Sweave() on my .Rnw file, it simply halts when it gets to the code chunk with the bad code. I can use \verbatim{} instead, but then I don't