similar to: how do I read only specific columns using read.csv or other read function

Hello, I am puzzled by the behavior of hist() when generating multiple plots per page on the pdf device. In the following example two pdf files are generated. The first results in 4 plots on one pdf page as expected. However, the second, which swaps one of the plot() calls for hist(), results in a 4 page pdf with one plot per page. How might I get the histogram with 3 other scatter

Hi list, I have a question on how to identify the overlapped items from two vectors: >x<-c(0,1,2) >y<-c(1,2,2,3) and plot the number of the overlapped as well as non-overlapped in a diagram? Thanks much, Allen [[alternative HTML version deleted]]

I have a beginner question. After I finally get the data to a data.frame that I can work with I have the following a data frame that is fairly long: > length(r2007) [1] 17409 If I look at the first element: > r2007[1] $`100009` DayOfYear Quantity 1 66 1 2 128 1 3 137 1 4 193 1 Now how do I get the length of this list (actually it is

Dear R Gurus: What is the name of the person who has all of the Tcl/tk stuff, please? I know it's James W, but can't remember his last name. I wanted to look at some of his examples. Thanks in advance! Sincerely, Erin -- Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown mailto: erinm.hodgess at gmail.com

Hi: I have a matrix, named "mat," of numeric values that are in days since Jan 1, 1960. I want to convert them Date class values so that I can do things like seq( value , len=10, by="1 week") where "value" is any of the matrix elements. How do I do this? Please reply to: philipsmith at alumni.albany.edu Thanks! Phil Smith Duluth, GA

Hi R-People: I use 2 machines: a machine with a Windows XP operating system, and another with a Linux Ubuntu OS. I transport my code between these 2 machines. However, pathnames to data files always need to be "adjusted" to account for the OS that I'm working on. Here is my question: How do I write code to detect whether I'm using the XP or the Linux machine? If I knew

Hi All, I am dealing with a large data set which translates in to a sparse matrix, I want to load this data that is spread in approximately 17000+ files each defining a row and each file has variable number of records that come with its column number and the value that they store. I wish to load this data in memory from these files one by one. Is there anyway I can do this in R, before I start

I would like to revisit a problem that was discussed previously (see quoted discussion below). I am trying to do the same thing, using a string to indicate a column with the same name. I am making "foo" a string taken from a list of names. It matches the row where "item" = 5, and picks the corresponding "taxon" > foo <- list$taxon[match(5,list$item)] Let's

Hello all, I am trying to get the variance-covariance (VCOV) matrix of the parameter estimates produced from the nlminb minimizing function, using vcov.nlminb, but it seems to have been expunged from the MASS library. The hessian from nlminb is also producing NaNs, although the estimates seems to be right, so I can't VCOV that way either. I also tried using the vcov function after minimizing

This is the rsync version I am using. rsync version 2.5.5 protocol version 26 Copyright (C) 1996-2002 by Andrew Tridgell and others <http://rsync.samba.org/> Capabilities: 64-bit files, socketpairs, hard links, symlinks, batchfiles, IPv6, 64-bit system inums, 64-bit internal inums I have a file with a looooong name which I rsync from PC1 to PC2. (redhat-9) I am getting the

How to show panels for factor levels of conditioning variables which do have no values? E.g. there are panels for "Grand Rapids" when they have values: data( barley ) with( barley, dotplot(variety ~ yield | year * site, layout=c(6,2) ) ) There are no panels for "Grand Rapids" when there are no values for "Grand Rapids": my.barley <- subset( barley, ! ( site ==

Hello, I have generated 2 Poisson processes and want to plot them on a single graph in a step by step manner in order to be able to compare them. I tried plot and biplot but it does not help, I could connect two points by hand for point graph if they were 5 or 10 I have more than 200 such point to be connected and Poisson cluster makes it difficult for me to even read them properly. Can anyone

I have the following small dataset: > stardata NS HE EB CW RW PW 1 0 0.000 0.000 0.042 0.006 0 2 0 0.006 0.000 0.013 0.005 0 3 0 0.000 0.011 0.000 0.000 0 I have plotted the star diagrams as follows: stars(stardata, key.labels = dimnames(stardata)[[2]], labels = NULL, key.loc = NULL, draw.segments=TRUE, col.segments="gray", lty="blank") I am

Using the maptools function "map2poly" I have created a map/polylist object - a map of the Great Lakes. My goal is to passively plot additional data on this map. Unfortunately I am not able to change (distort) the relative scale of the axes for the map, and so the sample points do not line up correctly on the map. In the code below, notice that I have set the same axis ranges on each

Better would be to use interval censored data. Create your data set so that you have (time1, time2) pairs, each of which describes the interval of time over which the tag was lost. So an animal first captured at time 10 sans tag would be (0,10); with tag at 5 and without at 20 would be (5,20), and last seen with tag at 30 would be (30, NA). Then survit(Surv(time1, time2,

Good morning and Happy Friday. A Telecommunication/Cable client of ours in the Duluth, GA area has an immediate opening for a Mid-Level or Senior-Level Ruby on Rails Developer. Please see below for the job description: Apex Systems, the nation's 2nd largest IT Staffing organization, has an opportunity for a long-term contract Ruby on Rails Developer in Duluth, GA with a competitive

Hi R Users, I'm trying to re-order the "site names" ("Waseca", "Morris", ...). I'm using following code: libarry(lattice) barchart(yield ~ variety | site, data = barley, groups = year, layout = c(6,1), aspect=.7, ylab = "Barley Yield (bushels/acre)", scales = list(x = list(abbreviate = TRUE, rot=45,

Hello. I apologize in advance for the VERY lengthy e-mail. I endeavor to include enough detail. I have a question about survival curves I have been battling off and on for a few months. No one local seems to be able to help, so I turn here. The issue seems to either be how R calculates Kaplan-Meier Plots, or something with the underlying statistic itself that I am misunderstanding. Basically,

I can not get a laptop to running Win98 to mount samba shares on a Linux box, and I don't think it's the normal encrypted password problem. Here's the setup: Linux box has a class B ethernet connection to the internet and a class C connection to the laptop. +-------------------+ | Linux |--eth0--1.2.3.4---->internet | box

I know about the current.row, current.column, and panel.number functions that are useful within panel functions written for lattice. Are there easy ways to obtain the names of the conditioning variables (those appearing after |) and their values for the current panel? Thanks Frank ----- Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: