similar to: windows() options

I am wanting to create a model matrix and keep the NAs. stratmat <- model.matrix(myformula,mydata) Is there any way to do this? model.matrix doesn't have na.action as a parameter. Elsewhere I have made use of na.keep <- function(x){x}. Many thanks, Rachel Cunliffe -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hi there, I am wanting to create 8 side-by-side histograms which have been rotated 90 degrees clockwise from how they usually sit.. all with the same scales. Is someone able to help me out? Thanks so much, Rachel Cunliffe -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

Hi there Probably a very simple solution to this problem. I have a character vector eg c("name1","name2","name3") and I want to pull out these columns from a data.frame, converting each of these columns into factors also. Many thanks Rachel -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hi again Why does this occur? > test <- apply(data,2,factor) # where data is a data.frame > is.factor(test[,1]) [1] FALSE Many thanks, Rachel -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the

Hello, When I am writing in sqldf or RSQLite I lose the decimals in my matrix. The only way I can get decimals is by multiplying by 1.00000, etc. I have tried manipulating the options, but it is only effective once I multiply by 1.0000. I appreciate any suggestions! Thanks! Example: z <- sqldf ("select ST, SUM(AGEP*PWGTP)*1.000000000/SUM(PWGTP)*1.00000000000000 as wgtage from ss09

I've tried Zorrin, Ubuntu and Centos. Zorrin and Ubuntu work fine, but Centos doesn't load, or run. I asked for advice, didn't get any and seem to have been added to all your Centos mailing lists. Please remove me from your mailing lists. Sincerely, Bob Cunliffe -------------- next part -------------- An HTML attachment was scrubbed... URL:

Ok, I finally got the slideshow code to a state worth showing it off. The site is a very rough cut of a site I''m building for my wife''s photography, so ignore the unfinished design for now :) http://rachel.kathihill.com/ To see the ajax version, go to: http://rachel.kathihill.com/?ajax=1 To randomize the order the images show: http://rachel.kathihill.com/?random=1 To change

I just tried downloading the lattice package from CRAN but got a syntax error when loading the library (see below). > {a <- CRAN.packages() + install.packages(select.list(a[,1],,TRUE), .lib.loc[1], available=a)} trying URL `http://cran.r-project.org/bin/windows/contrib/PACKAGES' Content type `text/plain; charset=iso-8859-1' length 7505 bytes opened URL downloaded 7505 bytes trying

Hi All, I want to take a simple random sample from a large dataset, gly, but I'm getting an error message. Any help? dim(gly) [1] 112371 37 > s1 <- sample(gly,100) Error in `[.data.frame`(x, .Internal(sample(length(x), size, replace, : cannot take a sample larger than the population when 'replace = FALSE' Thanks, Rachel [[alternative HTML version

Quick question please....A user logs into windowsXP and tries to create a folder/document and the ownership on the new file/folder defaults to nobody:nobody. I have the user set up in samba on the IRIX machine. All other users have no problem. Anyone have any suggestions? Thanks Rachel

Hi All, I'm trying to automate a data summary using summary or describe from the HMisc package. I want to stratify my data set by patient_type. I was hoping to do something like: Describe(myDataFrame ~ patient_type) I can create data subsets and run the describe function one at a time, but there's got to be a better way. Any suggestions? Rachel [[alternative HTML

Hi, I am trying to create a function for a goodness-of-fit test for the Pareto Distribution for some loss data that I have. So far I have the following: function(X=OTOL) { n <- length(X)-1 #calculated the number of values (extra as 0 included) i <- 2:640 #values of i j <- 1:639 #values of i-1 Y <- (n-j+1)*((X[i])-(X[j])) #First part of GoF model Y } Where OTOL is the ordered loss

Hi All, I have a factor variable with 52 levels -with long, annoying names. I want to keep only rows with some variables. I can do this using this code: test1 <- subset(nih2009,ic_name %in% c('NATIONAL EYE INSTITUTE','Veterans Affairs')) dim(test1) [1] 2396 38 But this doesn't work: t1 <- subset(nih2009, ic_name %in% c(27,51)) dim(t1) [1] 0 38 I know

I have attempted to follow posting guidelines but I have failed to find out what I am doing wrong here. I am trying to use lines.survfit to plot a second curve onto a survival curve produced by plot.survfit. In my case this is to be a progression free survival curve superimposed upon an overall survival curve, but I will illustrate my problem using the example given in the help for

[global] workgroup = IIG netbios name = TUX server string = Samba Server interfaces = 192.168.0.10/24 update encrypted = Yes logon path = \\TUX\profile\%u logon drive = H: domain logons = Yes os level = 65 preferred master = Yes domain master = Yes wins server = 127.0.0.1 wins support = Yes

I''m a rails newb and have been Googling about this, but I''m still stumped. Not showing everything here, but basically it should be a pretty common setup so I''m sure others know what I''m doing wrong. - A meter can belong to many meter_groups - A meter_group can have many meters. - A user can ''subscribe'' to viewing a meter_group (Subscription)

Hi, I would like to view matrices I am working with in a clean, easy to read, separate window. A friend showed me how to do something like I want with edit(). I can view the matrix in the 'R Data Editor': For a sample matrix: > mat=matrix(1:15,ncol=3) > mat [,1] [,2] [,3] [1,] 1 6 11 [2,] 2 7 12 [3,] 3 8 13 [4,] 4 9 14 [5,] 5 10 15

Hi! Whenever I save a graphic in ps/eps format generated by filled.contour or image, an undesired grid is added to it (not visible on the X11 screen). For example: postscript("volcano.eps") filled.contour(volcano,col=gray(seq(0,1,,50)),levels=seq(min(volcano),max(volcano),,50)) dev.off() Any ideia how to eliminate this grid? Thanks, Rachel [[alternative HTML version deleted]]

I posted several months about the problem with adding tick marks to curves using lines.survfit. This occurs when lines.survfit is used to add a curve to survival curves plotted with plot.survfit. The help for this function implies that mark.time=TRUE thus: plot(pfsfit,conf.int=FALSE,xscale=365.25,yscale=100,xlab="Years",ylab="% surviving",lty=2,mark=3)

Hello, I am trying to deal with a database table that contains a date in the form of the number of seconds since epoch. Data is inserted in this format by an existing script that I can''t change, so a schema change isn''t realistic. I would like to allow humans to view and edit this data in the form of an actual human-readable date. I figured a facade column would be the