search for: rhelp10

...here is some documentation from r-help over the years). Note 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 w...

Hi all, I am trying to fit a logistic regression for a bivariate response using five independent variables in a stepwise procedure. My outputs look okay but does any one know (or is there any literature on) how the confidence intervals are calculated for the reported odds ratios..? Thanks! [[alternative HTML version deleted]]

Dear all, I am confused with the output of survfit.coxph. Someone said that the survival given by summary(survfit.coxph) is the baseline survival S_0, but some said that is the survival S=S_0^exp{beta*x}. Which one is correct? By the way, if I use "newdata=" in the survfit, does that mean the survival is estimated by the value of covariates in the new data frame? Thank you very much!

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

Dears! I have a system that generates links on a https webserver. I have to be logged on with a username and password; Then follow at least fife links to get into the link where I can download the file. The question is: How can I accomplish that with R? Caveman

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

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

Hi, My questions concerns the ROCR package and I hope somebody here on the list can help - or point me to some better place. When evaluating a model's performane, like this: pred1 <- predict(model, ..., type="response") pred2 <- prediction(pred1, binary_classifier_vector) perf <- performance(pred, "sens", "spec") (Where "prediction" and

Hi list, I want to draw a bar plot with color indicating one grouping and different shading on top of the color indicating another grouping.? How should I proceed? Thanks! ...Tao

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):

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 = c(7, 8, 6, 4, 5, 4, 10, 12, 7, 13, 4, 11, 4, 3, 3, 10, 8, 16, 11, 6, 3, 5, 6, 13, 15, 2, 6, 11, 11, 10, 6, 9, 11, 14, 10, 7), TrSeasYr = c("1", "1", "1", "1",...

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

...ne gets such an error message? And why? Many thanks, I am desperate :-( I got this error on R 2.11.1 on Mac OS X 10.6.4, and on R 2.10.1 on Windows XP SP3 (but it does not matter, since I cannot cook a reproducible example). Philippe P.S.: seems related to this: http://finzi.psych.upenn.edu/Rhelp10/2008-June/165101.html -- ..............................................<?}))><........ ) ) ) ) ) ( ( ( ( ( Prof. Philippe Grosjean ) ) ) ) ) ( ( ( ( ( Numerical Ecology of Aquatic Systems ) ) ) ) ) Mons University, Belgium ( ( ( ( ( .............................................

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

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&...

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

...ith Jewell "Keith Jewell" <k.jewell at campden.co.uk> wrote in message news:... > 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))]. > > More recently a colleague had problems using the 'FactoMineR...

Dear All, Consider the following trivial code snippet rm(list=ls()) name_vec <- c("color1", "color2") pdf("test_color.pdf") plot(seq(5), seq(5), main=paste(name_vec[1]," and ",name_vec[2], sep="")) dev.off() What I would like to achieve is rather simple to explain, but it is giving me a headache: how can I have two colors in main? Let us

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, 231.67,110.62,173.11,202.85,689.27,23.82,234.15,625.89,192.88,29.53,109.99, 80.53, 46.85,82.69,231.53,75.38,217.52,85.29,182.18) ltlnorm <- function(x, meanlog, sdlo...

I'm trying to identify and remove rows in a data frame that are duplicated only on particular columns within it (i.e. not on all columns). The "unique" function looks for uniqueness across all columns of a data frame. Identifying unique rows based only on specific columns of interest returns only those columns, not all of the columns in the original frame. I tried this, and then