similar to: nlme with a factor in R 2.4.0beta

Hi, the following R lines work fine in R 2.4.0, but not in R 2.4.1 or any devel versions of R 2.5.0 (see below for details). library(drc) # to load the dataset 'PestSci' library(nlme) ## Setting starting values sv <- c(0.43355869, 2.49963220, 0.05861799, 1.73290589, 0.38153146, 0.24316978) ## No error m1 <- nlme(SLOPE ~ c + (d-c)/(1+exp(b*(log(DOSE)-log(e)))), fixed =

Dear R-Help, I have used R2.6.0 until I recently installed also R2.7.2 (see details below) In R 2.6.0, the following script using the function 'var' (cor(stats)): x.test <- c(NA, NA, NA, NA) var(x.test, na.rm=T) gives the output: NA In R2.7.2 the output of the same script generates an error message and stops R: 'Error in var(x.test, na.rm = T)

Dear list members, is it intentional that: curve(cos, xlim = c(-5, 5)) plot(cos, xlim = c(-5, 5)) produce different plots? Shouldn't the 'xlim' argument in both cases set the 'from' and 'to' argument if they aren't supplied (at least that's what I understood reading the help page for 'curve' and 'plot.function')? Christian >

Dear friends - I am trying to understand fda and working with the 2009 book from Springer. I run the scripts directly from a new installed fda library and include sessionInfo() below. This is from the script fdarm-ch09.R - I seem to have got not fda_2.2.8 but 2.2.7 - is that the problem? Where do I get .8 as I just reinstalled? Best wishes Troels Ring, Nephrology Aalborg, Denmark >

Dear R-listers, I am using R for Windows Vista 64-bit. I have no experience working with maps, but now I have to plot a map of Denmark in which each Danish county is color-coded according to its incidence rate of a particular disease. My problem is that I don't know where to start. I have read the help files to packages like 'mapplot', and their examples indicate that I am on the

Dear R-listers, I am trying to assign colnames to a data frame within a loop, but I keep getting a "target of assignment expands to non-language object"-error. I need to split up a large dataset into about 20 smaller ones, and I would like to assign colnames within the loop, so I won't have to type the same thing 20 times over. I have concocted this really goofy example which

Dear list, Please consider the following call of sort > sort(c("a","f")) [1] "a" "f" > sort(c("f","a")) [1] "a" "f" > > sort(c("aa","ff")) [1] "ff" "aa" > sort(c("ff","aa")) [1] "ff" "aa" The last two results look strange

Hi, I am writing an Sweave document and am using 'xtable' to make frequency tables of diagnoses of people undergoing cholecystectomy. Some of these diagnoses contain Danish characters ("?", "?", and "?"), and these characters are all garbled in the Latex document after I run Sweave. The odd thing is, everything looks absolutely right in the R console, and if

Dear R-listers, I am building a prediction model starting with many, many variables. I use the 'stepAIC' procedure in the MASS package, and the model building process takes several hours to complete. At the very end, I am occasionally confronted with warnings like this one: Warning messages: 1: In fitter(X, Y, strata = Strata, offset = offset, weights = weights, : Loglik converged

I have created an R-package with datasets which I want my students to install (the package is not on CRAN). 1) I've put the package on the web in a directory called 'data' and I thought I could do: > install.packages("http://gbi.agrsci.dk/statistics/courses/2011-ISMLS-course/data/LiSciData_0.0-03.tar.gz",repos=NULL,type="source") Warning: invalid package

The example in the help file for pdf() plots the characters outside the plotting area in this for loop: for(i in c(32:255)) { x <- i y <- i points(x, y, pch=i) } The following loop seems to be working as intended for(i in c(32:255)) { x <- (i-31)%%16 y <- (i-31)%/%16 points(x, y, pch=i) } As an extra

And as follow up, fff<-function(x) gsub(",", ".", x) BX <- apply(apply(AX, 2, fff), 2, as.numeric) this seems to be easier. Cheers Petr S pozdravem | Best Regards RNDr. Petr PIKAL Vedouc? V?zkumu a v?voje | Research Manager PRECHEZA a.s. n?b?. Dr. Edvarda Bene?e 1170/24 | 750 02 P?erov | Czech Republic Tel: +420 581 252 256 | GSM: +420 724 008 364 petr.pikal at

Dear friends - I have a data.frame with "," instead of "." and found the discussion http://r.789695.n4.nabble.com/How-to-replace-all-commas-with-semicolon-in-a-string-tt4721187.html#a4721192 so copying the code of Ulrik(I hope:-)) I tried (making some data) AX <-

The following dump-core patches changes its format into ELF, adds PFN-GMFN table, HVM support, and adds experimental IA64 support. - ELF format Program header and note section are adopted. - HVM domain support To know the memory area to dump, XENMEM_set_memory_map is added. XENMEM_memory_map hypercall is for current domain, so new one is created. and hvm domain builder tell xen its

We can use the setfiles(8) command to relabel the guest filesystem, even though we don't have a policy loaded nor SELinux enabled in the appliance kernel. This also deprecates or removes the old and broken SELinux support. This patch isn't quite complete - I would like to add some tests to the new API. I'm posting here to garner early feedback. Rich.

v1 -> v2: - Add simple test of the setfiles API. - Use SELinux_relabel module in virt-v2v (instead of touch /.autorelabel). - Small fixes. Rich.

Hi everyone: Is there a function in R to calculate the confidence intervals for the predictions of a GLS(Generalized Least Square) model? The function "predict" gives confidence intervals for the predictions of other types of models (lm, glm, etc) but not gls. Any input will be much appreciated, Best, Camilo Camilo Mora, Ph.D. Department of Biology Dalhouisie University

Indexing mail attachments with Dovecot + Solr. This patch has been tested with these versions: * dovecot 2.0.9 * apache-solr 1.4.1 This is a patch for the fts-solr plugin (that indexes mail messages for Dovecot with Solr). In main stream, the plugin does not index attachments; With this patch, you can index mails and their attachments (pdf, docs, openoffice docs...) . You can get others

You need to give the model formula that gave your output. There are two sources of variation (at least), within and between locations; though it looks as though your analysis may have tried to account for this (but if so, the terms are not laid out in a way that makes for ready interpretation. The design is such (two locations) that you do not have much of a check that effects are consistent over

This question will probably get me in trouble on theoretical grounds, but I will pose it anyway. The situation: I recently ran a field study looking for differences in sugarbeet cultivar tolerance to a specific herbicide. The study was set up so that 37 cultivars were treated with 4 different applications of the herbicide (37*4 factorial). In doing so, we found that the interaction effect was