search for: dn2

...alled ## "PSPP". Actually, the executable is "psppire". ## ## The other thing I do is open the same data again in ## a numeric format, and compare the 2 combined data frames ## This is also a useful exercise because it helps you ## understand what a "factor" is in R. dn2 <- read.spss("ESS3e03_2.por", use.value.labels = F) dn2 <- as.data.frame(dn2) dn2$whichSurvey <- 2 dn3 <- read.spss("ESS2e03_1.por", use.value.labels = F) dn3 <- as.data.frame(dn3) dn3$whichSurvey <- 3 ## Might be smart to compare # dn2$HAPPY[1:50] # d2$HA...

...e regarding rtol and atol would be appreciated. p <- c(r1=.01, r2=.01, a12=0.5, a21=0.6, k1=1,k2=1) t <- 1:20 lvcomp2 <- function(y, t, p) { # Lotka - Volterra model of two competing populations dN1.dt <- p["r1"] * y[1] * (1-(y[1] + p["a12"]*y[2])/p["k1"]) dN2.dt <- p["r2"] * y[2] * (1-(y[2] + p["a21"]*y[1])/p["k2"]) list(c(dN1.dt,dN2.dt),c(???)) } outlv <- lsoda(c(.1,.1),t, lvcomp2, p,rtol=1e-4, atol=1e-6) outlv TIA, Henry Dr. M. Henry H. Stevens Postdoctoral Associate Department of Ecology, Evolution, & Natural...

>From my book on corpus linguistics with R: # (10) Imagine you have two vectors a and b such that a<-c("d", "d", "j", "f", "e", "g", "f", "f", "i", "g") b<-c("a", "g", "d", "f", "g", "a", "f", "a",

...I data. The data being visualized consists of a series of layers (mostly representing emissions from a particular tuple of agricultural {cultivar, cultivation technique}) over a surface (in the upper midwest US). I'm currently plotting https://github.com/downloads/TomRoche/ioapi-hack-R/compare.DN2.layers.pdf each layer, and including some information about the layer in a title and subtitle. The title is where I want it--above the image--but the subtitle is printing below the image; I'd prefer it between the title and the image. Can this be done? FWIW, the code that plots (in https://gi...

I have a North American BRI configured as National ISDN (NI-1) on an SBC (Nortel?) switch. SBC calls this "FastTrak ISDN", and it's surprisingly inexpensive -- less than two POTS lines. I've been trying to find an inexpensive PCI interface solution to connect this to Asterisk as voice lines. I prefer ISDN BRI over POTS lines because of the improved signalling, fast dialing,

Consider the following scrap of code: > x<- ts(1:50,start=c(1,11),freq=12) > y <- aggregate(x,nfreq=4) > c(y) [1] 6 15 24 33 42 51 60 69 78 87 96 105 114 123 132 141 > y Error in rep.int("", start.pad) : invalid number of copies in rep.int() > tsp(y) [1] 1.833333 5.583333 4.000000 So we can aggregate into quarters, but we cannot print it using

Consider the following scrap of code: > x<- ts(1:50,start=c(1,11),freq=12) > y <- aggregate(x,nfreq=4) > c(y) [1] 6 15 24 33 42 51 60 69 78 87 96 105 114 123 132 141 > y Error in rep.int("", start.pad) : invalid number of copies in rep.int() > tsp(y) [1] 1.833333 5.583333 4.000000 So we can aggregate into quarters, but we cannot print it using