similar to: gremlin in rep()

********************************************************** ** Dr. David Lucy ** ** Centre for Forensic Statistics and Legal Reasoning ** ** Department of Mathematics and Statistics ** ** The University of Edinburgh ** ** James Clerk Maxwell Building ** ** King's Buildings

here is a distillation of a problem encountered in transfering some working code from R-0.63 to R-1.1.0 a1 <- 1:10 b1 <- tapply(a1,a1, sum) c1 <- ts(b1) c1 Error in if (NCOL(x) == 1) { : missing value where logical needed note that the error is returned as the value of calling ts() and is not automatically displayed problem seems to be that is.array(b1) returns TRUE yet dim(b1)

The following works as expected plot(1:5, panel.first=grid(2,2)) and if my.panel<-function() grid(2,2) then plot(1:5, panel.first=my.panel() ) is also OK but, do.call("plot", list(x=1:5, panel.first=grid(2,2)) do.call("plot", list(x=1:5, panel.first=my.panel)) do.plot("plot", list(x=1:5, panel.first=my.panel() ))

Hi I have installed Festival on the same box as asterisk and followed the instructions to integrate it with asterisk. Festival seems to work fine on its own performing text to speech from the command line or via a file. Asterisk answers the call but there is no speech. I can see no errors in the Festival log file The asterisk console shows --Executing Answer("SIP/81801-c091",

Hello, We have recently added a 64-bit 2 x Dual-Core Athlon server running Red Hat Enterprise Linux AS release 3 to allow for processing of large data sets (>4GB) in R. To integrate this server into our Linux cluster, I have been trying (unsuccessfully) to use Rmpi/lam-mpi to parallel process some of our scripts. While I have successfully compiled R 2.4.0 and lam-mpi 7.1.2 on this server using

Hi all, May be there is a plot which shows 3 variables in a triangle, their sum being constant. I have forgotten its name, and a search in the engine using "plot" did not help. Does anyone know of such a procedure in R or S? Thank you --christian Dr.sc.math.Christian W. Hoffmann Mathematics and Statistical Computing Landscape Dynamics and Spatial Development Swiss Federal Research

Use of the graphical parameter 'sub' in a high level graphics function is inconsistent, and incompatible with the blue book > plot(1:10, main="test", sub="demo") Warning: parameter "sub" couldn't be set in high-level plot() function Warning: parameter "sub" couldn't be set in high-level plot() function Warning: parameter "sub"

I have come across three problems in the past few days, in spell of heavy R ( version R0.50-a4/Sun Solaris2.5.1 ) use. 1. I was using lwd=2 to get thicker lines on plots for printing, but although the 'lwd' parameter works fine with x11(), the thickened lines do not print with print.plot, or by using postscript() directly. 2. Try the following, plot(1:10, -(1:10))

Setup: Sparcstation 4 with Solaris 2.5.1 and Sun compilers After downloading R-0.16.1 and runnning configure, make fails as follows cd main; make cc -o R.binary arithmetic.o array.o attrib.o bind.o builtin.o character.o coerce.o colors.o context.o cov.o cum.o dataframe.o debug.o deparse.o dotcode.o dstruct.o duplicate.o edit.o envir.o errors.o eval.o format.o gram.o

Dear All, I am trying to produce log-log plots in R and I was wondering if any of you have a 'template' for generating these with 'nice' labels and log-log grids? I know I can set up axes individually and use the intervals I want, however, I will be producing a large number of these plots and would not like to do this manually for each of them + I am very new to R and at the

A trick I use is to give tab-delimited text files created by excel or R the .xls extension: this may apply, and it may be better, to CSV: R couldn't see any difference whichever the extension will be, and windows will be foolished by the .xls extension an will open the file with Excel that will handel it CORRECTLY. If you need to have the data always alligned in Excel and R this trick will

I'm running Windows NT4. Rcmd does not seen to read the command line. Rterm and other things work well. For example "Rcmd check --help" fails with the command interpreter trying to find a program or file called "check." What have I missed? -- Bob Wheeler --- (Reply to: bwheeler at echip.com) ECHIP, Inc.

Hi, Say I have this vector x <- sort(rnorm(n)). I'd like to know if there is a nice way to find out the index of the nearest x to say 1.5? Thanks YG -------------- next part -------------- An HTML attachment was scrubbed... URL: https://stat.ethz.ch/pipermail/r-help/attachments/20000625/be875885/attachment.html

Hi, Is there any way to change the y axis range of values in a plot.gam()? I need that two different GAM plots to be of the same scale. Also, it is possible to change the labels? I tried with "ylab" and "ylim" and did not work Thanks in advance Ricardo Lopes Ricardo Lopes ............................................. Instituto do Mar Departamento de Zoologia

Hello. I have searched the archives but have not found anything. I need to solve a constrained optimisation problem for a nonlinear function (“maximum entropy formalism”). Specifically, Optimise: -1*SUM(p_ilog(p_i)) for a vector p_i of probabilities, conditional on a series of constraints of the form: SUM(T_i*p_i)=k_i for given values of T_i and k_i (these are constraints on

Hi, is there a way in R to generate a polynomial spline with multiple predictor variables? I have one response and two predictors and I'm trying to fit a spline model for this... Please cc me on the reply.. Thanks, nirmal

Hi, I've got the R library "mgcv" for GAM written by Simon Wood which works well in many instances. However, over the years I got attached to the S+ implementation of GAM which allows loess smoothing in more than 1 dimension as well as spline smoothing. Has anyone ported the S+ GAM library to R? Regards, Doug Beare. Fisheries Research Services, Marine Laboratory, Victoria Road,

All, Sorry for bombarding you with GAM related questions, but... I know a partial residual option in plot.gam() is on Simon Wood's todo list, but since I'm in the midst of a project and not yet having acquired sufficient R knowledge to code something usable myself I'll have to put my trust in you. Anybody got some code lying around for doing this? Or if someone can supply me with

hello, I am looking for a step() function for GAM's. In the book Statistical Computing by Crawley and a removal of predictors has been done "by hand" model <- gam(y ~s(x1) +s(x2) + s(x3)) summary(model) model2 <- gam(y ~s(x2) + s(x3)) # removal of the unsignificant variable #then comparing these two models if an significant increase occurs. anova(model, model2,

Dear R Users, I would like to make a lesast squares regression similar to that what is done by the command "lm". But additionally, I would like to impose some restrictions: 1) The sum of all regression coefficients should be equal to 1. 2) Each coefficient should assume a value between 0 and 1. (inequality restrictions) Which command is the best to use in order to solve this problem