similar to: avoiding (another) error

Dear all, is this a bug or do a miss something (this concerns version 0.63.1 on Linux)? > boxplot (c (1, 2, 3, 4, 5), ylim= c (-5, 5)) Warning: parameter "ylim" couldn't be set in high-level plot() function > boxplot (c (1, 2, 3, 4, 5), pars= list (ylim= c (-5, 5))) Warning: parameter "ylim" couldn't be set in high-level plot() function > and actually in

Dear all, well, today I just seem to run into these problems ... (I am still working on R 0.63.1 on a SuSE 5.2). I would like to have a data.frame that also contains character (and NOT factor) vectors. Two things that I tried (the second worked in R 0.62.2): > test.df <- data.frame (X= character (4)) > test.df$X [1] <- "a" Warning: invalid factor level, NAs generated or:

Dear all, I am working with R Version 0.62.2 (July 10, 1998) on a SuSE Linux 5.2. I have troubles in understanding the function image (). Can someone enlighten me? I'd like to draw a square image-plot (with every little part of the image a square as well) and add some (straight) lines to it. Also, I create an empty z matrix which I need to fill step by step. Thus I tried to find out how

Dear all, The following works with a simple result of lm but gives the error with a more complicated one (could this be due to singularities that occurr in the model?). > cooks.distance (X) Error: NAs in foreign function call (arg 6) > class (X) [1] "lm" But, actually, I would like to ask a more general question: Is there a place (I searched in the CRAN site) where I find

Dear all, I think this question has been raised very recently but I think I didn't see repsonses on the list nor did I find the thread in the archives. (If it is there I'd appreciate directions ...) I am about to install Linux on a Compaq notebook. One important thing I want to do is being able to run R. I have fiddled with different distributions of Linux (LST, SuSE, Debian) and can

Dear all, in a current project I have a pretty huge .Rdata. Thus I was working with R --vsize 100 --nsize 1000000. Today when I tried to restart R I get the following error message: Error: a read error occured Fatal error: unable to restore saved data (remove .RData or increase memory) I increased memory up to --vsize 180 and --nsize 2000000, but the error reoccurrs. Is there a way to know

Dear all, I would like to fit an lm in which a subset of the explanatory variables are linearly dependent. Thus I would like to include the restriction that all betas of these variables sum to 1. Is there a way to this in R? Or is this what happens automatically if I set singular.ok= T (which is the default, I believe)? Thanks for your help. Lorenz --

Dear all, I am trying to use ace on a 446 x 14 Matrix, with two categorical variables in the first two columns. If I try to run it I get: > aceiH <- ace (data.mat, resp, cat= 1:length(x.cat)) Error: NAs in foreign function call (arg 3) There is only one .Fortran call in ace () with the third argument being x = t(x). Thus I tried to debug ace, i.e. I went through ace until just before

Dear all, I would like to use the biological symbols for male and female as plotting symbols in a scatterplot (ideally filled and non-filled). R does not seem to have these symbols using pch= in plot() nor are they implemented via expression() or at least I did not find them. I found that the symbols are e.g. available in the wasysym and the marvosym package in LaTeX. I have coded two very rough

Dear all, I would like to use the biological symbols for male and female as plotting symbols in a scatterplot (ideally filled and non-filled). R does not seem to have these symbols using pch= in plot() nor are they implemented via expression() or at least I did not find them. I found that the symbols are e.g. available in the wasysym and the marvosym package in LaTeX. I have coded two very rough

Hi, Another question about graphics... Suppose I have: par(mfrow = c(4, 6)) which will set up a figure region for 24 plots. But, I actually only need 23 plots. In fact I wrote a little function, foo.plot, which will draw all the 23 plots at once with: foo.plot(x) where x is a data frame. Now, suppose I want to draw several set of these 23 plots, but want to put them all in one postscript

I am using lme to handle repeated measurements. So far i can follow the examples from the book from pinheiro and bates. Now i get the problem , that i have "nested" repeated measuremnts, and i cant find out how to do the grouping part of the lme formula. 1.The simple problem ist that i have different Samples , from which i make repeated measurements (each sample is measured 6 times) and

I have a data.frame with three columns + the row number Sample MidPoint Count [1,] 1 0.025 1 [2,] 1 0.075 3 [3,] 1 0.125 15 [4,] 1 0.175 12 [5,] 1 0.225 5 [6,] 1 0.275 8 When I export the data to a csv the column names are shifted because the row numbers are being exported as well. Sample MidPoint

I don't think, this has been answered: > I'm trying to run a 3-way within-subject anova in lme with 3 > fixed factors (Trust, Sex, and Freq), but get stuck with handling > the random effects. As I want to include all the possible random > effects in the model, it would be something more or less > equivalent to using aov > > > fit.aov <- aov(Beta ~ >

Dear all, I do not seem to grasp how contrasts are set for ordered factors. Perhaps someone can elighten me? When I work with ordered factors, I would often like to be able to reduce the used polynomial to a simpler one (where possible). Thus, I would like to explicetly code the polynomial but ideally, the intial model (thus, the full polynomial) would be identical to one with an ordered factor.

Hello I just installed openssh-3.7.1p1 on IRIX64 6.5 04100802 IP27. The compilation failed since the header file openbsd-compat/inet_ntoa.h was missing. I downloaded the package from different sources including ftp.openbsd.org, there was no difference. Since inet_ntoa.c has not been changed, copying the inet_ntoa.h from openssh-3.6p1 fixed the problem. The package should be fixed anyway.

Does someone out there have an example of R-code for a split-split plot ANOVA using aov or another function? The design is not balanced. I never set up one in R before and it would be nice to see an example before I tackle a very complex design I have to model. Thanks, Mike Mike Saunders Research Assistant Forest Ecosystem Research Program Department of Forest Ecosystem Sciences University of

I am trying to perform a model 3 ANOVA for a 2 factor (say factor A and factor B) anova in which factor A is fixed and factor B is random. Therefore, the error term for the test of factor A should be the A:B interaction term and the error terms for B and A:B should be the model residual (within) term. I have tried to work out how to specify such error strata using aov, however, I have had

I do not think anyone has answered this. > I'm trying to run a simple one-way ANCOVA with the lmer > function in R package lme4, but have encountered some > conceptual problem. The data file MyData.txt is like this: > > Group Subj Cov Resp > A 1 3.90 4.05 > A 2 4.05 4.25 > A 3 4.25 3.60 > A 4 3.60 4.20 > A 5 4.20 4.05 > A 6 4.05 3.85

Can someone verify for me if the for loop below is really calculating the nonzero min for each row of a matrix? I have a bug somewhere in the is section of code. My first guess is how I am find the the nonzero min of each row of my matrix. The overall idea is to make sure I am investing all of my money, i.e. new.set is a set of indicator variables for each stock for a particular portfolio, i.e.