similar to: Thougt I understood factors but??

data: > tmp1 Date HrMn Temp Q.4 1 19450101 0900 -37.0 1 2 19450101 1000 -35.9 2 3 19450101 1100 -35.9 3 4 19450101 1200 -36.4 4 5 19450101 1300 -36.4 5 6 19450101 1400 -36.4 6 7 19450101 1500 -36.4 7 8 19450101 1600 -37.5 9 Accidentally, I did this (I meant to write Q.4 instead of Q here) > tmp1$Q[tmp1$Q!="1" & tmp1$Q!= "5"]<-NA I would

Dear all, I stumbled upon what appears to be a troublesome issue when sampling from an ARIMA model (from Rob Hyndman's excellent 'forecast' package) that contains a seasonal AR component. Here's how to reproduce the issue. (I'm using R 2.9.2 with forecast 2.19; see sessionInfo() below). First some data: > x <- c( 0.132475, 0.143119, 0.108104, 0.247291, 0.029510,

Dear all, I've recently purchased a new 64bit system with an intel i7 quadcore processor. As I understood (maybe wrongly) that to date the 32bit version of R is more stable than the 64bit, I installed the 32bit version and am happily using it ever since. Now I'm running a whole lot of models, which goes smoothly, and I thought out of curiosity to check how much processor I'm using. I

It looks like a timezone issue, and it's causing confusion to me at least. My original data: gmt <- c("19880101 0000", "19880101 0100", "19880101 0300", "19880101 0400", "19880101 0500", "19880101 0600") These were converted to local dates/times with akst<-strptime(gmt,format="%Y%m%d %H%M")-(3600*9) # because I want

Dear R list, I get a strange error in R: *** caught segfault *** address 0x6951c20, cause 'memory not mapped' Traceback: 1: .C("spline_eval", z$method, nu = as.integer(n), x = as.double(xout), y = double(n), z$n, z$x, z$y, z$b, z$c, z$d, PACKAGE = "stats") 2: spline(gam.data$x[, col.data], gam.smooths.all$fit[, m], xout = gam.results.global[m, ,

Hi folks, I am confused whether the following is a bug or it is fine Here is the explanation a <- zoo(c(NA,1:9),1:10) Now If I do rollapply(a,FUN=mean,width=3,align="right") I get > rollapply(a,FUN=mean,width=3,align="right") 3 4 5 6 7 8 9 10 NA NA NA NA NA NA NA NA But I shouldn't be getting NA right ? i.e for index 10 I should get (1/3)*(9+8+7)

Hi R users, I'm using the survey package to calculate summary statistics for a large health survey (the Demographic and Health Survey for Honduras, 2006), and when I try to calculate the variances for several variables, I get negative numbers. I thought it may be my data, so I ran the example on the help page: data(api) ## one-stage cluster sample dclus1<-svydesign(id=~dnum, weights=~pw,

Hello I tried the odfWeave package today, by running the formatting.odt and example.odt files that are included with the package. They both ran fine, but when I try to open them in my OpenOffice (OpenOffice 3.1.1 on Kubuntu 9.10) I get an error "Format error discovered in the file in sub-document content.xml at 1293,124(row,col)." I also tried to open them in MS Word 2003 (Windows

Dear All I have a problem with the height of the boxes in the key in the following. (The text is over 2 lines to accentuate the problem of no space between the rectangles.) Is there an easy way to put a space between the rectangles; size controls the width but there appears to be nothing for the height? xyplot(1~1, key = list(corner = c(0.8,0.8),

Hi R users, Im trying to install XML source package (in window platform) using the following command. *install.packages(pkgs = "C:\\Documents and Settings\\skang\\Desktop\\Softwares\\R\\Library\\XML_2.6-0.tar.gz", repos = NULL, type = "source")* and get the following result. gzip: stdin: unexpected end of file /cygdrive/c/Rtools/bin/tar: Unexpected EOF in archive

Dear all I am trying to fit Variable Coefficients Models on Unbalanced Panel Data. I managed to fit such models on balanced panel data (the example from the "plm" vignette), but I failed to do so on my real, unbalanced panel data. I can reproduce the error on a modified example from the vignette: > require(plm) > data("Hedonic") > Hed <- pvcm(mv ~ crim + zn + indus

The integral of any probability density from -Inf to Inf should equal 1, correct? I don't understand last result below. > integrate(function(x) dnorm(x, 0,1), -Inf, Inf) 1 with absolute error < 9.4e-05 > integrate(function(x) dnorm(x, 100,10), -Inf, Inf) 1 with absolute error < 0.00012 > integrate(function(x) dnorm(x, 500,50), -Inf, Inf) 8.410947e-11 with absolute error <

Hello All I have some problem for installing XML_2.6-0.tar . I am working in widows and R version is R-2.9.1 >*install.packages("XML")* After selecting a CRAN mirror ** *Error :-* Warning: unable to access index for repository http://cran.pt.r-project.org/bin/windows/contrib/2.9 Warning: unable to access index for repository

Dear R developers, I find the format of descriptions of the arguments in the read.table help text slightly inconsistent. For example, the logical arguments comes in seven different formats, more or less explicit about the consequences of a TRUE (or FALSE): 1. check.names logical. If TRUE then the names... 2. blank.lines.skip logical: if TRUE blank lines... 3. flush logical: if TRUE,

Hi all, I am just curious how to sort the groups in a categorical box plot chart bwplot here is the example: d<-data.frame(sample(rep(month.abb,20), 100), runif(100,1,10)); colnames(d) <- c("Month", "Value"); bwplot(d$Month ~ d$Value); as you can see Months are not sorted alphabetically. Does anybody know ho to sort those on the chart level? Thank you Jan

Hi, Doing 'R CMD check diveMove' is now throwing this message: Data codoc mismatches from documentation object 'sealLocs': Variables in data frame 'sealLocs' Code: id.time.class.lon.lat Docs: class id lat lon time with: R> sessionInfo() R version 2.10.1 (2009-12-14) x86_64-pc-linux-gnu locale: [1] LC_CTYPE=en_CA.UTF-8 LC_NUMERIC=C

Dear Lists: I recently ran quite annoyance problem while running R on Ubuntu 9.10. When running the program, the system suddenly exit from the R session with the following warnings: ############################ ############# OMP: Hint: This may cause performance degradation and correctness issues. Set environment variable KMP_DUPLICATE_LIB_OK=TRUE to ignore this problem and force the program

Hi, I summary() a variable with 409908 numeric observations. The variable is part of a data.frame. The problem is that the min and max returned by summary() do not equal the ones returned by min() and max(). Does anybody know why that is? > min(data$vc) [1] 15452 > max(data$vc) [1] 316148 > summary(data$vc) Min. 1st Qu. Median Mean 3rd Qu. Max. 15450 21670 40980

Dear list, since I have upgraded openOffice to version 3.2 I have some trouble to open very simple ODT files generated by odfweave: the file is apparently corrupted (but recovery is fine). I have observed this under windows and mac OS 10.4.11 with R 2.10.0, odfWeave_0.7.11, XML_2.6-0, lattice_0.17-26. However, I found nothing similar by RSiteSearching so that I'm suspecting a localization

Hi, I'm trying to get data from web page and modify it in R. I have a problem with encoding. I'm not able to get encoding right in htmlTreeParse command. See below > library(RCurl) > library(XML) > > site <- getURL("http://www.aarresaari.net/jobboard/jobs.html") > txt <- readLines(tc <- textConnection(site)); close(tc) > txt <- htmlTreeParse(txt,