similar to: R-alpha: Treatment of seed

The following commands crash: > plot(1:4,axes=F,log="x") > axis(1,c(1,2,3,4)) zsh: segmentation fault R Lines 35-37 of axis(): 35: else { 36: ind <- (usr[1] <= at & at <= usr[2]) 37: } should be replaced by else { if (log) ind <- (10^usr[1] <= at & at <= 10^usr[2])

Here are another three problems with logarithmic scales: 1) segments() does not work with logarithmic scales. I suggest to change lines 962-973 in "plot.c": for (i = 0; i < n; i++) { if (FINITE(xt(x0[i%nx0])) && FINITE(yt(y0[i%ny0])) && FINITE(xt(x1[i%nx1])) && FINITE(yt(y1[i%ny1]))) { GP->col = INTEGER(col)[i % ncol];

hello. I have a query about the Matrix package for R. I wrote some code a while ago using the Matrix package version 1.6.2 with an early version of R, to do some linear least squares for regression: xn [,1] [,2] [,3] [,4] [1,] 1 0.7429352 0.5519528 0.4100652 [2,] 1 0.7443713 0.5540886 0.4124477 [3,] 1 0.7447385 0.5546355 0.4130584 [4,] 1 0.7459597 0.5564558

This problem has been mentioned by Arne Kovac <maak@stats.bris.ac.uk> on May 9-10. The patch is really a replacement of 'axis' in RHOME/src/library/base/funs/ It contains the fixes which where proposed Arne Kovac and fixes (all?) the errors he/she (?) reported on May 9. "axis" <- function (which, at, labels = TRUE, ...) { if (which%%2 == 1) { axp <-

Dear list, Could somebody who is more experienced with the Tcl/Tk interface from R please help me clarify the issues I've put below with ### --> tags? Several things go wrong, and it's probably because of messy code, but I have a difficult time finding out what is the cause. Thanks very much, JeeBee. require(tcltk) || stop("Package tcltk is not available.") # Add path to

Gareth Did you ever get a solution to this? mary Mary Owen IT Manager Stride Treglown Limited Promenade House, The Promenade, Clifton Down, Bristol BS8 3NE T 0117 9157361 F 0117 974 5207 E maryowen@stridetreglown.com www.stridetreglown.co.uk

There probably is a way, but it involves some programming. You would need to clone a smooth constructor (e.g. for the "cr" class), and then modify it to add a linear constraint matrix C to the returned smooth object. If b are the smooth coefficients then C should? be the matrix such that s(0) = Cb (you can get this from the Predict.matrix method for the class). Then the constraint

Dear All I am trying to do a repeated measures analysis using lmer and have a number of issues. I have non-orthogonal, unbalanced data. Count data was obtained over 10 months for three treatments, which were arranged into 6 blocks. Treatment is not nested in Block but crossed, as I originally designed an orthogonal, balanced experiment but subsequently lost a treatment from 2 blocks. My

Thanks Simon, by cloning a smooth construct do you mean copying and modifying the smooth constructor code? Could you pleas elaborate on your answer? Which is the Predict.matrix method? 2018-01-12 3:20 GMT-06:00 Simon Wood <simon.wood at bath.edu>: > There probably is a way, but it involves some programming. You would need to > clone a smooth constructor (e.g. for the "cr"

Brian, Thanks for doing this. As the primary author, could I have a desk copy, please? Bill Venables. Bill Venables, CMIS, CSIRO Marine Laboratories, PO Box 120, Cleveland, Qld. 4163 AUSTRALIA Phone: +61 7 3826 7251 Fax: +61 7 3826 7304 Mobile: +61 419 634 642 <mailto: Bill.Venables at csiro.au> http://www.cmis.csiro.au/bill.venables/ -----Original Message----- From: Brian

Dear all How can I obtain the name of the argument passed in a function? Here is a simplistic example of what I would like to obtain: myfunction= function(name) { print(paste("The parameter name was",unknownFunction(name)) } myfunction(myobject) [1] "The parameter name was myobject" Thanks Francisco

Hi, If any one has time I need some help understanding the P-values given in the lmer output. Using AIC for model selection I find my minimal model is FOLLOW~MOVERSTATUS+DISTANCE however it appears DISTANCE is not significant at 95% confidence, see output quoted below. However, removing DISTANCE gives a higher AIC=433.5, therefore I will keep it in, but am confused as to what is adds to the

Hello, I have stored R objects as hexadecimals in a mysql database, I then usually transform them to binary and then save it into a file. E.g. hex <- getBlob #gets blob from database as hexadecimal binary <- transformToBinary(hex) #moves from hex to binary I would usually save them into a file as follows: writeBin(object=binary,con="fileName.txt") Then eventually if

Hi all, There are quite a few of us working on LLVM in the Bristol area. I invite you all to the first Bristol LLVM social: Where: The Hope & Anchor 38 Jacobs Wells Road Bristol BS8 1DR http://www.hopeandanchor.net When: Thursday 20 September from 7pm Look forward to meeting you all there. Jeremy -- Tel:

Hi, I'm running a GLMM on binomial choice data. The outputs I receive are sensible except for the degrees of freedom, which come out much larger than expected. Can anyone advise please? Exptl design: Response = Choice Fixed Factors = Position, Treatment and Sex Random Factor = ID nested within Treatment and Sex Covariate = Delay The model: glmmPQL(FreeChoice ~ Position * Treatment + Sex

Dear list, I have been using multivariate multiple regression (MMR) in the form lm(Y~X) where Y and X are matrices of response and predictor variables. I know that summary(mlm.object) would give the usual lm statistics for each response variable separately and that anova.mlm(mlm.object) will give the analysis of variance table of the mlm object. However, anova.mlm (also manova(mlm.object))

Hi, I noticed that the R reference manual seems to be missing the license permissions on the copyright page. Below is a patch to "refman.top" which adds the same permissions as the @permission{} macro in "R-defs.texi". I'm filing this bug as I have been preparing to publish the R Reference Manual (in several volumes). I plan to donate 5 EUR to the R-foundation for each

Hi all, I am very new to R and this might be a simple question but I have looked everywhere you suggest before writing to you. I am trying to change font type from san-serif to a serif (Times New Romans) on all labels and axis of my boxplot. I have used this function in other plots before, e.g.: plot(residuals~lnlifespan, data=mydata, pch=psymb, font=6, xlab="ln reproductive

Hi, Below I've flagged some sentences in the R documentation (from cvs) which are garbled in some way. If anyone can clarify what these sentences should say I'd be grateful. The problems were picked up as part of the proof-reading of the R manuals by Ciaran O'Riordan. After looking at the text neither of us could figure out what the documentation was meant to say. -- Brian Gough

Hi, I need to generate two normal variables with covariance matrix: 0.25, 0.20 0.20, 0.25 but I have no idea how to do that. Can anyone help me? Thanks, Joaquim ---------------------------------------------------------------------------- ---------- Joaquim J. S. Ramalho University of Bristol Department of Economics 8 Woodland Road Bristol BS8 1TN United Kingdom