similar to: R-help Digest, Vol 128, Issue 5

Dear R-list, It is not uncommon for laboratory equipment (e.g. spectrophotometers) to have a linear response in a certain interval and then go into saturation. I wonder if there is an R-function that models this; for instance by estimating the breakpoint and fitting a line below the breakpoint and a constant above. Best regards Thomas Poulsen -- Thomas Poulsen Research Scientist. PhD,

Re: Heteroscedasticity and mgcv. (Collin Lynch) The GAMLSS package can model heterogeneity in the scale parameter (e.g. standard deviastion) [and also heterogeity in skewness and kurtosis parameters].of the response variable distribution. For parametric models a generalized likelihood ratio test can be used to test whether the heterogeity is needed. Alternatively a generalized Akaike

Has anyone had trouble installing the package "ROpenOffice" from omegahat: http://www.omegahat.org/ROpenOffice/ I downloaded the .tar.gz file to my local hard drive and tried to install it with the following code, but got an error message: > install.packages("/Users/earlbrown/Desktop/ROpenOffice_0.4-0.tar.gz", repos = NULL, type = "source") Installing package(s)

Dear All, I have the following problem when reading files (a lot of them) in the spreadsheetML format into R. The spreadsheetML format is an xml format to allow easy import of multisheet data in Excel. As far as I can see, a direct import into R (using the XML package) is not feasible. I use the software ssconvert (included in Gnumeric) and call it from R. It converts the spreadsheetML into xls

Hi, I had some excel files that I needed to read into R, there were alot so I didn't want to do them by hand in gnumeric. I tried all the recommendations in the data import/export manuals, including RODBC, and for some reason they all failed on these files. Then I stumbled on ssconvert, a script that wraps all gnumerics converters, see

Hello, for my data preparation and administration (data, labels, etc.) I use OpenOffice.org ODS spreadsheet files with several sheets in one file. However, I find it inconvenient to export every single sheet to a csv file whenever I apply changes to the labelling or so, and I haven't found a plugin or an application which does this batch in OOorg. Is there new development on a direct import

Dear List, Please bear with a poor newbee, who might be doing everything backwards (I was brought up in pure math). I want to make a simple multi-linear regression on a set of data. I did some expreiments, and if X is a 4 by 2 array and Y is a 4 by 1 array, I can do a linear regression by lm(y~x). Now I have a tab-delimited text file with 10 rows of 300 measurements and an other file with

Dear list, I would like to call a Kerberos-authenticated web-service from within R. Curl can do it: $ curl --negotiate -u : "http://my.web.service/" so I would expect that RCurl also has the capability, but I have not been able to find the correct options to set. listCurlOptions() does not return anything with negotiate, and searching the source of RCurl, the only thing I found was

I would like to open OpenOffice (LibreOffice) .ods files in R. I've tried the ROpenOffice package from omegahat, but unfortunately, the read.ods() function attempts to use the values of the first column in a worksheet as row names, and thus does not allow duplicates in there (which, even more unfortunately, occur in my files). Also, the function does not allow to forward any other

Hi I need to use rep() to get a vector out, but I have spotted something very strange. See the reproducible example below. N <- 79 seg <- 5 segN <- N / seg # = 15.8 d1 <- seg - ( segN - floor(segN) ) * seg d1 # = 1 rep(2, d1) # = numeric(0), strange - why doesn't it print one "2"? rep(2, 1) # 2, ok rep(2, d1 / 1,1) # 2, this

Hi list! For a package providing a GUI intended in particular at R newcomers (RcmdrPlugin.temis[1]), I would like to use the Rstem package [2] by Duncan Temple Lang, which is hosted on Omegahat only. This means the package cannot be installed automatically as a dependency, which I can work around by installing the package from my code when the user needs it. But Omegahat does not provide binaries

Is there a way of executing SAS codes in R environment?   Peter Maclean Department of Economics UDSM [[alternative HTML version deleted]]

I have tried to compile the working Fortran-90 routines (see attachment) to interface them with R. The compiler automatically caled by R CMD command line prints a bunch of error messages (as follows). I think I'm better off translating the Fortran implementation into R ... which is what I hoped I could avoid. Maura host705:Desktop mauede$ R CMD SHLIB ford.f90 gfortran -arch i386

I am estimating Value at Risk using PerfomanceAnalytics package. The?variables are stored in a data frame. I formated the data variables using zoo() and as.xtx() but it is not working. The working example is below. ##########################################################? reguire(zoo) require(PerformanceAnalytics) reguire(xts) ? year<- c(1991-12-30, 1992-12-30, 1993-12-30, 1994-12-30) R1

I have csv data that extend beyond the area I want for an existing map. I want using the boundaries of the polygon shape file as a cookie cutter so that I can overlay the csv data on map without including anything outside the map boundaries and create a dbf file or shapefile of the clipped data .  The reproducible example: ############################################### library(RColorBrewer)

Hi Mackay and anybody (a) Is it possible to select randomly (let say five grids) and plot? (b) Is it possible to plot five nearest grid in one figure? The original question and improved codes: #I am ploting gridded time series data. I would like the actual lat #and lon value appear on the graph-if possible inside the graph as #numbers. If there is also more elegant ways to plot the graphs I

Hi, Can R read open office.org Calc files Thanks -- Shane [[alternative HTML version deleted]]

I want a list of the number of times some factor levels appear together, similar to the following SQL statement: SELECT A, B, COUNT(C) FROM TBL GROUP BY A, B How do I do that with a data.frame in R? Thanks, Jesper Frickmann Statistician, Quality Control Novozymes North America Inc. 77 Perry Chapel Church Road Franklinton, NC 27525 USA Tel. +1 919 494 3266 [[alternative HTML version deleted]]

Yes it is. Example from my program: workdir <- "S:\\Stat\\C03202\\R-output\\" ... sink(file=paste(workdir, assayname, "-1.txt", sep="")) ... win.metafile(filename = paste(workdir, assayname, "-2%02d.wmf", sep="")) You get the idea! Kind regards, Jesper Frickmann Statistician, Quality Control Novozymes North America Inc. Tel. +1 919 494

I am trying to estimate a gamma function using real data and I am getting the following error messages. When I set a lower limit; the error message is "L-BFGS-B needs finite values of fn" ? For other method the error message is: Error in optim(x = c(0.105286666666667, 0.3472275, 2.057625, 0.329675,? : ? non-finite finite-difference value [1] The codes works fine for simulated data