similar to: CART analysis

I all I am just trying to get the idea of a loop if statement for another purpose. I was going to label a new variable based on date ie first week 1, second week 2 etc. My code in bold is wrong but seems logical to me, could someone help. x <- rep(c(2660156,2663703,2658165,2659303,2661531,2660914),c(2,2,2,2,1,1)) y <- rep(c(6476767,6475013,6475487,6479659,6477004,6476388),c(2,2,2,2,1,1))

Hi all I hope you might help me with some aspects of producing a graph in lattice. There are three things I have struggling with and that is: 1. to separate the horizontal box rows from each other; 2. to change the colour of the horizontal and vertical strips to white; and 3. to place the axes labels on the left y axes and on the bottom x axes. I would really appreciate some help. I have put the

Hi R users I apologise in advance for this question as I suspect it is simple and perhaps others have had this problem. I am struggling to sort out how to fix the x axes so that the labels don't overlap. I have put the following example together to show my problem. library(lattice) titre <- as.factor(rep(c(10999,20999,30999,40999,50999,60999,

Hi all I know this must be an easy one so sorry for the trouble. I would like to select a list of variables within a factor The following example is given in help for subset: subset(airquality, Temp > 80 select = c(Ozone, Temp)) So how do I select all temperatures of 90 and 80 ie Temp = c(80,90) I would appreciate your help. Kind regards andy Andrew McFadden MVS BVSc Incursion

Hi R users I am trying to create a spatial join between two datasets. The first data set is large and contains descriptive data including x and y co-ordinates. The second dataset is small and has been selected spatially. The only data contained within the second dataset is the x and y coordinates only i.e. no descriptive data. The aim of a join made between the two datasets is to select

I all I am trying to combine columns from two dataframes to make a completely new dataframe consisting of one column of dates (ie a combination of dates from the two dataframes). >From the following dataframes a b 1 2008-07-27 1 2 2008-10-01 2 3 2008-08-15 3 4 2008-08-14 4 5 2008-08-14 5 6 2008-09-20 6 c d 1 2008-07-27 1 2 2008-10-01 2 I would like to get: z

Hi all The code of trying to write relates to selecting properties (given by x and y co-ordinates) spatially (distance "X" from "infected" properties identified by date) over a certain time period. i.e. what properties are within 3 km from properties infected on "2008-01-01" over the last 14 days. Is any one able to give me some clues on how to write code to solve

Hi all I am trying to do a simple recode which I am stumbling on. I figure there must be any easy way but haven't come across it. Given data of A","B","C","D","E","A" it would be nice to recode this into say three categories ie A and B becomes "Treat1", C becomes "Treat 2" and E becomes "Treat 3". I tried

Hi R users The scatterplot (used package R-Cmdr ) obtained has the legend off the plot in the upper left I would like to decide where to put the legend - is that option still available and how do I do it? code: Dataset <- read.table("J:/Pigmatingsb.txt", header=TRUE, sep="\t", na.strings="NA", dec=".", strip.white=TRUE)

Hi all I have created a lattice graph but the text appears too big. Just wondering if there was a way I could change the font size. tp1.nDai=xyplot(value ~ yr | Name, data = nDai,type = 'o') update(tp1.nDai, aspect = c(0.8), layout = c(4, 4)) I would appreciate your help. Kind regards andy Andrew McFadden MVS BVSc Incursion Investigator Investigation & Diagnostic Centres -

Hi R users > I am trying to add a regression line to a graph for "c" for factor 2 > only. Any suggestions? > library(lattice) a=(1:10) b=c(1:5,16:20) c=as.factor(c(rep(1,5),rep(2,5))) d=data.frame(a,b,c) xyplot(a~b, pch=c(6,8),data = tavg, groups=d$c, reg.line=lm, smooth=FALSE, type=c("p","g"),xlab="a",ylab="b") > I would

Hi R users I trying to some aggregate statistics on a very large dataset. The null values are causing a problem. I would like to calculate aggregate values for groups. I am just no sure how to deal with the "" I would ideally like them to ignored ie if there is only 1 value over several columns than that value would be the summary statistic. I have put the following example together to

Hi R users I am trying to find unmatched data from two dataframes. I would like to find unmatched data based on several factors. For the following data: dat1 <- data.frame(x = paste("A", 1:6, sep=""), y = c("andy","bob","ciaran","dan", "eion", "fred")) dat1 dat2

Hi all I was wondering if it is possible to get rid of the horizontal strips and produce each barchart with a left y axes and lower x axes only. Also can you specify an exact size of graph ie 88mm wide with a font size of 'x'. library(lattice) library(latticeExtra) n=as.factor(c(1:5,1:5))

. Monday, July 5, 2004 15:50:06 (-08:00hrs UTC) Hello asterisk-users, From the following post: On Mon, 5 Jul 2004, Joe Baptista wrote: > On Mon, 5 Jul 2004, Randy Bush wrote: >> i did not criticize the protocol. remember, my question started >> with >> >> >> i am looking at iax to

Full_Name: Andrew Thornton Version: Webpage OS: Linux Submission from: (NULL) (81.154.9.63) The Writing R Extensions documentation at: http://cran.r-project.org/doc/manuals/R-exts.html lists two different filenames for making the example's Rout.save a test result. Section 1.1.3 says: "If tests has a subdirectory Examples containing a file pkg-Ex.Rout.save, this is compared to the

Try something like with(df, predict(smooth.spline(x = altitude, y = atm_values), deriv = 1)) Cheers, Andrew -- Andrew Robinson Chief Executive Officer, CEBRA and Professor of Biosecurity, School/s of BioSciences and Mathematics & Statistics University of Melbourne, VIC 3010 Australia Tel: (+61) 0403 138 955 Email: apro at unimelb.edu.au Website: https://researchers.ms.unimelb.edu.au/~apro

Greetings R Developers, I've been testing the alpha release of R-3.0.0 and I noticed that the plotting functionality in package tkrplot was not working correctly. Further diagnosis found the issue to be related to the command '.Tcl("image create Rplot plotname")' from package "tkrplot". The following example is taken from here:

Hi, I'm thinking about adding support for PDU (Power Distribution Unit) devices from Digipower [1] and Aviosys [2] to NUT via the SNMP interface. The specific devices I have in mind are the following: Digipower AMz-1623/SW-08-1 Aviosys IP Power 925x (both are 8-port PDUs which are SNMP capable via an ethernet interface) Does anybody have other DigiPower or Aviosys devices and if so

Hi all I am trying to calculate the instantaneous hazard (rather than cumulative hazard) for an event "bird death". Some sample data is as follows: year<-c(rep(2003,10),rep(2004,10)) time<-c(5,10,5,15,20,10,15,15,25,25,5,15,15,5,25,10,25,15,35,25) event<-c(1,0,1,1,0,1,1,0,1,0,1,0,0,1,0,0,0,0,1,0) Ideally I would like to know if the hazard is greater at specific time periods