similar to: adding a built-in drop.levels option for subset() in 2.12 ?

Hello. In this post: http://finzi.psych.upenn.edu/Rhelp10/2010-March/233815.html Uwe Ligges suggests using BRugs rather than R2WinBUGS under windows. He also notes that it is not in the main CRAN repository, but it is in "extras" which is a default repository under windows. I have OpenBUGS 3.1.0 (the latest that has a native Linux version which is no longer called linBUGS)

All - I have an uneven set of replicates and would like to sample from this set X number of times to generate a mean for each grouping variable. I was thinking the boot package would be the thing to use, but now I'm not so sure ... given the discussion here: http://finzi.psych.upenn.edu/Rhelp10/2010-June/243828.html Given these data (a small subset): > dput(x) structure(list(numSpp

Hello, It's several days I try to track this bug, and even cannot cook a reproducible example. Yet, it occurs consistently in a long-running task after a variable period of time. Here is an example: ... my long-running code [as I said, cannot give something simple that produces this bug in a reproducible manner] Error in match(x, table, nomatch = 0L) : formal argument

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

There is a thread on this topic already: http://finzi.psych.upenn.edu/Rhelp10/2010-August/250934.html I'm rather mystified by a similar problem and wondering whether I've overlooked something obvious. I'm running with latest versions of R and all packages, and latest version of JAGS running under Windows 7. Here's the problem. I have some source code. It's given below -

I have been trying to get the lme4 package installed on Mac OS X... with no success. The Mac OS binary is not available on any CRAN, and I also can¹t install the package from old source. Has anyone found a solution to this problem? I am happy to use nlme for now, but I tend to prefer to do my mixed model analyses with lmer(). Thanks, Nicholas Griffin Center for Genome Sciences Washington

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

Hi! I'm trying to implement in R in an easy way the negative reflection method described in Silverman (1986) on p.31, ie I have a non-negative dataset and would like to estimate the density by applying a reflection method where the reflected points have weight -1. I thought there should be a way to evaluate the wanted kernel at the required points, in other words to calculate formula (2.16):

Dear List, When I subset a data.frame, the levels are not re-adjusted (see example). Why is this? Am I missing out on some basic stuff here? Thanks Ulrik > m <- data.frame(gender = c("M", "M","F"), ht = c(172, 186.5, 165), wt = c(91,99, 74)) > dim(m) [1] 3 3 > levels(m$gender) [1] "F" "M" > s <- subset(m, m$gender ==

Hi, Back in June I posted the message below, but had no replies. I've made a little progress since then so this is to update anyone interested (!) and to ask for comments. Brief problem statement: Under Windows, some parts of R don't handle UNC paths beginning with backslashes. Specifically a) Sys.glob() fails to find some files breaking (e.g.) Rcmdr plugins

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

Hi R Users, I have two vectors, x and y, of equal length representing two types of data from two studies. I would like to test if they are similar enough to use them interchangeably. No assumptions about distributions can be made (initial tests clearly show that they are not normal). Here some result: Two-sample Kolmogorov-Smirnov test data: x and y D = 0.1091, p-value < 2.2e-16 alternative

Hi, a colleague of 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

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

Hi, Back in 2010 I had a problem with 'update.packages()', which I worked around by mapping a drive letter to a UNC path [described in <http://finzi.psych.upenn.edu/Rhelp10/2010-February/229820.html> but my current workaround is assign(".lib.loc", sub("\\\\\\\\Server02/stats", "L:", .libPaths(), ignore.case = TRUE), env=environment(.libPaths))].

A bit of a perennial I'm afraid. My wife has filled up her iPhone 6 with photos, and wants them moved onto my computer. I'm running CentOS 7 patched about 30 minutes ago. Needless to say the computer can't see the data on the iPhone, though it does recognise the phone as an iPhone. Any suggestions (well any that don't involve a steam roller, sledge hammer or GBH to the whole of

Hi! I have following data which is left truncated say at 10. I am trying to estimate the parameters of the Truncated Lognormal distribution to this data as given below. (I have referred to R code appearing in an earlier post - http://finzi.psych.upenn.edu/Rhelp10/2008-October/176136.html) library(MASS) x <- c(600.62,153.05,70.26,530.42,3440.29,97.45,174.51,168.47, 116.63,36.51, 219.77,

I am interested in studying the binary representation of numerics (doubles) in R, so am looking for possibilities of output of the internal binary representations. sprintf() with format "a" or "A" is halfway there: sprintf("%A",pi) # [1] "0X1.921FB54442D18P+1" but it is in hex. The following illustrate the sort of thing I want: 1.1001 0010 0001 1111

Hello, Is there a function in R with which we could find the greatest common divisor directly? Or is there a funciton, which could help us to find the remainder of a division directly? Best wishes! Thank you! Your truly, Qing -- View this message in context: http://r.789695.n4.nabble.com/greatest-common-divisor-tp2326262p2326262.html Sent from the R help mailing list archive at