Re: [R] GLMM

2004-08-26 Thread Douglas Bates
Bossarte, Robert wrote: I am trying to use the LME package to run a multilevel logistic model using the following code: --- Model1 = GLMM(WEAP ~ TSRAT2 , random = ~1 | GROUP ,

Re: [R] Why terms are dropping out of an lm() model

2004-08-26 Thread John Pitney
John Pitney wrote: Hi all! I'm fairly new to R and not too experienced with regression. Because of one or both of those traits, I'm not seeing why some terms are being dropped from my model when doing a regression using lm(). I am trying to do a regression on some experimental data d,

Re: [R] Plotting groupedData objects

2004-08-26 Thread Deepayan Sarkar
Quoting David Scott [EMAIL PROTECTED]: On Thu, 26 Aug 2004, Deepayan Sarkar wrote: Quoting David Scott [EMAIL PROTECTED]: I am trying to create a plot similar to Figure 3.2 in Bates and Pinheiro. I have repeated measurements on about 80 subjects from 2 treatment groups.

RE: [R] Surprise when mapping matrix to image

2004-08-26 Thread Deepayan Sarkar
Quoting Glynn, Earl [EMAIL PROTECTED]: Prof Ripley: Thank you for your prompt reply. It's pure convention: see below. Did you try reading the help for image? You don't seem to understand it if you actually did. It seems you are looking for image(t(x)[ncol(x):1, ]) I

[R] Error TukeyHSD

2004-08-26 Thread T. Murlidharan Nair
I am running the following code on the coagulation data and I am getting an error. Please let me know if I am missing anything from my code. coag- matrix( scan(//Samba3/nair/R/blood.dat, sep=,), 24, 3, byrow=TRUE) colnames(coag) - c(time,diet,order) coag - as.data.frame(coag) oneway.test(time ~

[R] coplot and par

2004-08-26 Thread Mihai Nica
R 1.9.1 on Win2000 or Win98SE. I am using coplot as follows: coplot(AVG~LRPI| REGION) the output seems normal but I get: Warning message: calling par(new=) with no plot This is the only explanation that I have for being unable to use par() with coplot for changing the way the xlab and ylab

Re: [R] Surprise when mapping matrix to image

2004-08-26 Thread A.J. Rossini
I think I'd have to respectfully disagree with both Brian and Deepayan, as to whether it should be obvious. It is reasonable (principle of least suprise) to expect orientation of the plot to match the print order of the matrix. I would have expected Brian's one-liner to be in the help page,

[R] text() with text, variables and math HOWTO?

2004-08-26 Thread Johannes Graumann
Hello, One more question from the 'abusing R for blotting - particularly anally' department: How can I in the expression below make the '%~~%' show up as the aprrox-sign I want it to be? Thanks for any hint, Joh text( 500,1.5, cex=0.75, substitute(

[R] multiple regression with dummy variables

2004-08-26 Thread franc Li
Hi I'm a beginner in R, please help... I'm trying to do regression analysis on a categorical data to find b's of each Xi. Essentially, f(x) is the output, and two independent variables Xd (days of week (1,2,...,7) and Xwk (week of month (1,2,...,6)including partial wks). So the data looks

Re: [R] Surprise when mapping matrix to image

2004-08-26 Thread Jim Brennan
I was also surprised by the image orientation this summer and used the easy fix of matrix manipulation. There is however another issue and that is when you start flipping around the matrix, orders etc. and for the case you want to have sensibly labeled axes, you may have to use the axis commands

[R] expressions/approximate text

2004-08-26 Thread Erin Hodgess
Please forgive the general posting. To get the x approximately = y, please try text(6,2,expression(x %~~% y)) HTH. Sincerely, Erin Hodgess mailto: [EMAIL PROTECTED] __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

Re: [R] Label using equivalent of \mathbb{R}

2004-08-26 Thread Paul Murrell
Hi Deepayan Sarkar wrote: On Thursday 26 August 2004 11:41, Martin Maechler wrote: Simon == Simon Cullen [EMAIL PROTECTED] on Thu, 26 Aug 2004 15:30:06 +0100 writes: Simon On Thu, 26 Aug 2004 15:24:33 +0200, Trenkler, Simon Dietrich [EMAIL PROTECTED] Simon wrote: [Dietrich

[R] FIML in lme

2004-08-26 Thread F Z
Hi I was asked if lme can use FIML (Full Information Maximum Likelihood) instead of REML or ML but I don't know the answer. Does anybody know if this is implemented in R? Thanks Francisco __ [EMAIL PROTECTED] mailing list

RE: [R] Surprise when mapping matrix to image

2004-08-26 Thread Henrik Bengtsson
I've got some simple code example at http://www.maths.lth.se/help/R/image/, or just try: source(http://www.maths.lth.se/help/R/image/image.R;) Best wishes Henrik Bengtsson -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Jim Brennan Sent: Friday,

[R] read.spss warning: unrecognized record type??

2004-08-26 Thread Michael Jerosch-Herold
When using read.spss (library: 'foreign') I get the following warning message: Warning message: E:/R4win/mesamri.sav: Unrecognized record type 7, subtype 13 encountered in system file. I don't see anything wrong with record #7 in the database I am trying to read in, but I suspect that the

Re: [R] Plotting groupedData objects

2004-08-26 Thread David Scott
On Thu, 26 Aug 2004, Deepayan Sarkar wrote: Quoting David Scott [EMAIL PROTECTED]: On Thu, 26 Aug 2004, Deepayan Sarkar wrote: Quoting David Scott [EMAIL PROTECTED]: I am trying to create a plot similar to Figure 3.2 in Bates and Pinheiro. I have repeated

Re: [R] introduction slides for beginners

2004-08-26 Thread Deepayan Sarkar
Quoting Deepayan Sarkar [EMAIL PROTECTED]: On Thu, 26 Aug 2004, Armin Roehrl wrote: Hi all, does anybody have slides for a 1 to 3 hour crash course into R, that I would be allowed to recycle to show a few people the light? I have some stuff I used for 5 x 1 hour presentations at

RE: [R] Surprise when mapping matrix to image

2004-08-26 Thread Richard A. O'Keefe
Deepayan Sarkar [EMAIL PROTECTED] wrote: You seem to be thinking that Prof Ripley's solution had something to do with image(). Glynn, Earl [EMAIL PROTECTED] could certainly be excused for thinking so, because what Prof Ripley wrote included this: Did you try reading the

RE: [R] predict.mvr error message

2004-08-26 Thread Liaw, Andy
What version of R, what version of pls.pcr, and on what OS? Have you checked whether your versions of software are up to date? I get: n - 1350 p - 180 y - rnorm(n) x - matrix(sample(0:1, n*p, replace=TRUE), n, p) fit - mvr(x, y, method=SIMPLS, validat=none, ncomp=2) xt -

RE: [R] predict.mvr error message

2004-08-26 Thread Stewart T Chang
Well, I figured it out on my own. It turns out, even though my test set of data appeared to be a matrix, it wasn't. Adding as.matrix to the front of the name of the data set solved the problem. (Maybe there's some alternative way of doing it, such as yours?) Thanks for your suggestion.

RE: [R] Surprise when mapping matrix to image

2004-08-26 Thread Liaw, Andy
From: Richard A. O'Keefe [snip] Given the wide range of ways that the x, y, z arguments can be passed to image(), it would actually make sense to have some kind of flip and/or mirror operations specified via an argument to image(). The source code of image is available (image.default) so it

RE: [R] Surprise when mapping matrix to image

2004-08-26 Thread Deepayan Sarkar
Quoting Richard A. O'Keefe [EMAIL PROTECTED]: Deepayan Sarkar [EMAIL PROTECTED] wrote: You seem to be thinking that Prof Ripley's solution had something to do with image(). Glynn, Earl [EMAIL PROTECTED] could certainly be excused for thinking so, because what Prof Ripley wrote included

RE: [R] Surprise when mapping matrix to image

2004-08-26 Thread Richard A. O'Keefe
Deepayan Sarkar [EMAIL PROTECTED] wrote: I agree with you that the issues are related, but your (edited) quoting is misleading. My quotes were selective but not otherwise edited. I interpret Prof Ripley's comment as implying that reading ?image would tell the user

Re: [R] Surprise when mapping matrix to image

2004-08-26 Thread A.J. Rossini
Richard A. O'Keefe [EMAIL PROTECTED] writes: Once you know what image() is _intended_ to do with a matrix, you can figure out the transposing and reversing you need for any other view. It's knowing you need to that's the problem. This is a much better statement of my point, and precisely why

[R] Keyboard input into functions

2004-08-26 Thread Neil Leonard
Hi, Does anybody know if it is possible to have keyboard input into functions? Kind of like 'scanf()' in C. I want to write a function that asks me what to put in certain variables. Thanks, Neil __ [EMAIL PROTECTED] mailing list

RE: [R] Keyboard input into functions

2004-08-26 Thread Manoj - Hachibushu Capital
?readline HTH -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Neil Leonard Sent: Friday, August 27, 2004 2:04 PM To: [EMAIL PROTECTED] Subject: [R] Keyboard input into functions Hi, Does anybody know if it is possible to have keyboard input into

Re: [R] Difficulties in starting up with its package

2004-08-26 Thread Gabor Grothendieck
In looking at this again I don't think my x2 corresponds to yours. It should be: x2 - zoo(matrix(1:30,ncol=3),d)[as.numeric(format(d,%w)) %in% 1:5,] Gabor Grothendieck ggrothendieck at myway.com writes: : : Ajay Shah ajayshah at mayin.org writes: : : : : : Folks, : : : : I'm trying to

Re: [R] Order with scientific notation

2004-08-26 Thread Prof Brian Ripley
Your column of P.values is not numeric. If it were, it would be left justified and not printed with trailing zeroes suppressed. As in (X - data.frame(P.values=c(0.0379, 0.068, 0.0025, 8e-04))) P.values 1 0.0379 2 0.0680 3 0.0025 4 0.0008 So I guess it is factor or character, as in

Re: [R] Plotting groupedData objects

2004-08-26 Thread Deepayan Sarkar
Quoting David Scott [EMAIL PROTECTED]: I am trying to create a plot similar to Figure 3.2 in Bates and Pinheiro. I have repeated measurements on about 80 subjects from 2 treatment groups. I would like to have the panels for the two treatment groups in separate groups and within those

[R] A bug in ternaryplot (vcd)

2004-08-26 Thread Jong Hee Park
Does anyone know how to make multiple ternaryplots? My first figure in the following command appears in the bottom row and then the next figure goes to the next screen. par(mfrow=c(2,1)) ternaryplot(list[,1:3], grid=F, pch=.); ternaryplot(list[,4:6], grid=F, pch=.) I can't put these

Re: [R] postings and archives

2004-08-26 Thread Martin Maechler
Louize == Louize Hill [EMAIL PROTECTED] on Thu, 26 Aug 2004 01:34:15 +0100 writes: Louize PS - how long do postings take to get in the Louize archives - I remember seeing something about scale Louize values (pure numeric?) this week, I erased it and Louize now I

Re: [R] Beginners Question: Make nlm work

2004-08-26 Thread Arne Henningsen
Hi, did you try nls (nonlinear least squares)? ?nlm says: For nonlinear regression, 'nls' may be better. Another option would be optim. Best wishes, Arne On Thursday 26 August 2004 01:07, Johannes Graumann wrote: Hello, I'm new to this and am trying to teach myself some R by plotting

[R] Problems with par() and labels with boxplot

2004-08-26 Thread michael watson (IAH-C)
Quite a simple one really! When I run boxplot(), the labels on the X axis are horizontal, and I want them vertical. So I did: par(las=3) boxplot(...) And my labels just aren't there anymore Any help??? __ [EMAIL PROTECTED] mailing list

[R] o bty with different axes

2004-08-26 Thread bavorak
Is it possible to create plot bty=o with different scale of left- and right x-axis? In my documents, by typesetting, bty=u looks disruptively. At 2nd: how to frame the whole box (including title and axes labels) and resize to landscape-oriented rectangle? Tomas Bayer

Re: [R] Beginners Question: Make nlm work

2004-08-26 Thread Peter Dalgaard
Arne Henningsen [EMAIL PROTECTED] writes: Hi, did you try nls (nonlinear least squares)? ?nlm says: For nonlinear regression, 'nls' may be better. Another option would be optim. Actually, just using nlm according to specs might help: setwd(~/Biology/R_versuch)

Re: [R] Problems with par() and labels with boxplot

2004-08-26 Thread Uwe Ligges
michael watson (IAH-C) wrote: Quite a simple one really! When I run boxplot(), the labels on the X axis are horizontal, and I want them vertical. So I did: par(las=3) boxplot(...) And my labels just aren't there anymore Any help??? I guess they get clipped and you have to enlarge the margins,

RE: [R] Order with scientific notations

2004-08-26 Thread Ted Harding
On 26-Aug-04 Kevin Wang wrote: Hi, If a value like 8e-04 is in a data frame, is the following behaviour normal? final.df Chr P.values 11 0.0379 260.068 32 0.0025 4 138e-04 5 14 0.0244 63 0.0279 74 0.1561 8

Re: [R] o bty with different axes

2004-08-26 Thread Uwe Ligges
[EMAIL PROTECTED] wrote: Is it possible to create plot bty=o with different scale of left- and right x-axis? See the mailing list archives, this has been discussed several times during the last years. Why is this related to bty=o (which is the default)??? In my documents, by typesetting, bty=u

Re: [R] Problems with par() and labels with boxplot

2004-08-26 Thread Peter Dalgaard
michael watson (IAH-C) [EMAIL PROTECTED] writes: Quite a simple one really! When I run boxplot(), the labels on the X axis are horizontal, and I want them vertical. So I did: par(las=3) boxplot(...) And my labels just aren't there anymore Any help??? You may need to make

RE: [R] Problems with par() and labels with boxplot

2004-08-26 Thread michael watson (IAH-C)
Hello I thought it was such a simple problem it didn't warrant too much detail! Sorry :-) OS is Suse Linux 8.2, R version is 1.9.0, device is x11. Example reproducible code: data(InsectSprays) subset - InsectSprays[1:24,] # this creates normal horizontal text boxplot(count ~ spray, data =

RE: [R] Problems with par() and labels with boxplot

2004-08-26 Thread michael watson (IAH-C)
Yes, I thought of that, but when one is creating lots of images automatically, one doesn't always know how long one's labels are in advance. So I guess I need to check how long my labels are and truncate them if they are too large, or adjust the margins if they are just slightly bigger than

Re: [R] Problems with par() and labels with boxplot

2004-08-26 Thread Uwe Ligges
michael watson (IAH-C) wrote: Hello I thought it was such a simple problem it didn't warrant too much detail! Sorry :-) OS is Suse Linux 8.2, R version is 1.9.0, device is x11. Example reproducible code: data(InsectSprays) subset - InsectSprays[1:24,] # this creates normal horizontal text

RE: [R] Problems with par() and labels with boxplot

2004-08-26 Thread michael watson (IAH-C)
Nope. Even after par(las=3, mar=c(15,4,4,2)), it doesn't draw the labels, and at those margins the plot os all skwished up. I guess if I keep enlarging the margins then it maybe will draw my labels, but what use is that? What I would like it to do is simply try and squeeze the labels into the

[R] Column names of a data frame

2004-08-26 Thread michael watson (IAH-C)
I'm coming out with all the classic problems today ;-) Can anyone tell me how to change the column names of an already existing data frame? I've read the docs for ?data.frame and ?as.data.frame but can't figure it out. I want to make a new data.frame from some of the columns of an existing

RE: [R] Problems with par() and labels with boxplot

2004-08-26 Thread Prof Brian Ripley
On Thu, 26 Aug 2004, michael watson (IAH-C) wrote: The different behaviour between Windows and Linux is a bit odd though (see my previous post) No, it's not. It is the graphics device, not the OS which matters. How a graphics display system responds to a request to partially clip text is not

Re: [R] Column names of a data frame

2004-08-26 Thread Prof Brian Ripley
?names data frames have names, rather column names (they are lists with extra structure). On Thu, 26 Aug 2004, michael watson (IAH-C) wrote: I'm coming out with all the classic problems today ;-) Can anyone tell me how to change the column names of an already existing data frame? I've

Re: [R] Column names of a data frame

2004-08-26 Thread Sean Davis
Michael, See ?colnames. Sean On Aug 26, 2004, at 5:31 AM, michael watson (IAH-C) wrote: I'm coming out with all the classic problems today ;-) Can anyone tell me how to change the column names of an already existing data frame? I've read the docs for ?data.frame and ?as.data.frame but can't

Re: [R] integrate function

2004-08-26 Thread Jari Oksanen
On Wed, 2004-08-25 at 23:44, Peter Dalgaard wrote: Ronaldo Reis Jr. [EMAIL PROTECTED] writes: Is possible to integrate this diferential equation: dN/dt = Nr(1-(N/K)) in R using the integrate() function? No. However, you could use N = K/(1 + exp(log((K-N0)/K) -r*t)), where

Re: [R] Column names of a data frame

2004-08-26 Thread Arne Henningsen
On Thursday 26 August 2004 11:31, michael watson (IAH-C) wrote: I'm coming out with all the classic problems today ;-) Can anyone tell me how to change the column names of an already existing data frame? I've read the docs for ?data.frame and ?as.data.frame but can't figure it out. This is

Re: [R] Column names of a data frame

2004-08-26 Thread Marta Rufino
Hello, Can anyone tell me how to change the column names of an already existing data frame? I've read the docs for ?data.frame and ?as.data.frame but can't figure it out. names(object.name)=c(col.name1, col.name2, col.name3, etc.) # note that the number of names as to the the same as the

Re: [R] Beginners Question: Make nlm work

2004-08-26 Thread Peter Dalgaard
Peter Dalgaard [EMAIL PROTECTED] writes: Arne Henningsen [EMAIL PROTECTED] writes: Hi, did you try nls (nonlinear least squares)? ?nlm says: For nonlinear regression, 'nls' may be better. Another option would be optim. Actually, just using nlm according to specs might help:

[R] Re:how to compute a mode of a distribution

2004-08-26 Thread John Wilkinson
On 08/24/04 13:50, Paolo Tommasini wrote: Hi my name is Paolo Tommasini does anyone know how to compute a mode ( most frequent element ) for a distribution ? try median(density(my.vector)) John John Wilknson [EMAIL PROTECTED] __ [EMAIL PROTECTED]

[R] Image recognition in R

2004-08-26 Thread Adrian Baddeley
Angel Lopez writes: I have some images of bugs (insects) with many bugs in each image. I want to count the number of bugs and to have an estimate of the area of each one. See also the book A. Lawson D. Denison Spatial cluster modelling Chapman Hall where chapters 4-6 describe

RE: [R] image recognition in R

2004-08-26 Thread Stephen Eglen
Philippe wrote: You should use ImageJ to binarize the picture and then extract features (measurements) for each blob (individual). Save results as a text file and import it in R (read.table). There you have plenty of packages to analyze and classify your data. Things to try range from

Re: [R] image recognition in R

2004-08-26 Thread Barry Rowlingson
Angel wrote: I have some images of bugs (insects) with many bugs in each image. I want to count the number of bugs and to have an estimate of the area of each one. I've tried searching for an R package to do so with no success. Is this a task that I should pursue doing in R or

[R] Label using equivalent of \mathbb{R}

2004-08-26 Thread Simon Cullen
Hi, I'm trying to label the horizontal axis of a plot with a symbol that is the equivalent of \mathbb{R} in LaTeX. I've had a look through the help pages for plotmath and for Hershey and haven't found the symbol. Could someone give me a pointer, please? Using R 1.9.1 on Win32. -- SC

[R] library(car) Anova() and Error-term in aov()

2004-08-26 Thread Paul Lemmens
Dear all, Type III SS time again. This case trying to reproduce some SPSS (type III) data in R for a repeated measures anova with a betwSS factor included. As I understand this list etc, if I want type III then I can do library(car) Anova(lm.obj, type=III) But for the repeated measures anova, I

RE: [R] image recognition in R

2004-08-26 Thread Philippe Grosjean
Hello, Yes, UTHSCSA Image Tools is a good package too, but it is less powerful than ImageJ and it is not updated since quite a long time. There is an active development of ImageJ, plenty of plugins, and an easier interface to make its own plugins and scripts than in image tools. So, unless

Re: [R] Re:how to compute a mode of a distribution

2004-08-26 Thread Kjetil Brinchmann Halvorsen
John Wilkinson wrote: On 08/24/04 13:50, Paolo Tommasini wrote: Hi my name is Paolo Tommasini does anyone know how to compute a mode ( most frequent element ) for a distribution ? try median(density(my.vector)) That cannot be right. Try: test - rnorm(1000) d - density(test) median(d)

[R] PLS and categorical variables

2004-08-26 Thread David Galvan
Dear all, I am trying to analyze a Y matrix of fish abundanceof different fish species vs. an X matrix of environmental variables using PLS regression,implemneted inthe pls.pcr package.Some of my environmental variables are continuous and some are categorical. Is it possible tomix continuous

RE: [R] Label using equivalent of \mathbb{R}

2004-08-26 Thread Trenkler, Dietrich
-Original Message- From: Simon Cullen Sent: Thursday, August 26, 2004 1:26 PM To: [EMAIL PROTECTED] Subject: [R] Label using equivalent of \mathbb{R} Hi, I'm trying to label the horizontal axis of a plot with a symbol that is the equivalent of \mathbb{R} in LaTeX.

RE: [R] integrate function

2004-08-26 Thread Wayne Jones
I suggest that you check out the R package odesolve for solving this. However, from a more ecological standpoint (which I am guessing from the fact that it is logistical growth) you might be interested in SOLVER. You can find it at http://www.stams.strath.ac.uk/ecodyn/ Or indeed you might be

[R] introduction slides for beginners

2004-08-26 Thread Armin Roehrl
Hi all, does anybody have slides for a 1 to 3 hour crash course into R, that I would be allowed to recycle to show a few people the light? On the R-website http://cran.r-project.org/other-docs.html I found several additional books, but no slides. It might be nice to add a page with slides,

RE: [R] library(car) Anova() and Error-term in aov()

2004-08-26 Thread John Fox
Dear Peter and Paul, As Paul discovered, Anova() doesn't handle aovlist objects. As a general matter, one should be careful with type-III tests, since it's easy to test hypotheses that aren't sensible (e.g., tests ostensibly of main effects that aren't reasonably interpretable as tests of main

RE: [R] introduction slides for beginners

2004-08-26 Thread John Fox
Dear Armin, I have a variety of materials at http://socserv.socsci.mcmaster.ca/jfox/Courses/S-course/index.html, though the slides are minimal; most of the material is in the script files. The introductory lecture (of seven lectures) takes me two hours. I hope that this helps, John

RE: [R] Problems with par() and labels with boxplot

2004-08-26 Thread Deepayan Sarkar
On Thu, 26 Aug 2004, michael watson (IAH-C) wrote: Yes, I thought of that, but when one is creating lots of images automatically, one doesn't always know how long one's labels are in advance. bwplot in the lattice package will try to do that for you, which IMO is one of the big perks of

[R] EM norm package (NA/NaN/Inf in foreign function call (arg 2))

2004-08-26 Thread [EMAIL PROTECTED]
Greetings! I am bootstrapping and I am using EM in the norm package to fill in missing data for a financial time series with each step of the loop. For the most part EM works fine for me, but the following error message is guaranteed before I hit the 200th scenario: Iterations of EM:

RE: [R] Help using Hmisc / Latex

2004-08-26 Thread Frank E Harrell Jr
Sorry about non-indenting of quoted text - I haven't figured out a good way to answer mail when receiving messages in digest mode -FH I'm afraid you need to modify your approach. You're trying to pass latex an lm object, which

Re: [R] Label using equivalent of \mathbb{R}

2004-08-26 Thread Simon Cullen
On Thu, 26 Aug 2004 15:24:33 +0200, Trenkler, Dietrich [EMAIL PROTECTED] wrote: [Dietrich Trenkler] plot(rnorm(10),xlab=expression(bold(x)),ylab=expression(bold(y))) Not quite what I am looking for, I'm afraid. \mathbb gives blackboard fonts - the capitals with two vertical parallel

Re: [R] introduction slides for beginners

2004-08-26 Thread Deepayan Sarkar
On Thu, 26 Aug 2004, Armin Roehrl wrote: Hi all, does anybody have slides for a 1 to 3 hour crash course into R, that I would be allowed to recycle to show a few people the light? I have some stuff I used for 5 x 1 hour presentations at http://www.stat.wisc.edu/~deepayan/SIBS/slides/

Re: [R] Help using Hmisc / Latex

2004-08-26 Thread Gabor Grothendieck
Frank E Harrell Jr f.harrell at vanderbilt.edu writes: Sorry about non-indenting of quoted text - I haven't figured out a good way to answer mail when receiving messages in digest mode -FH What I do is post through http://news.gmane.org/gmane.comp.lang.r.general

[R] glmm in R

2004-08-26 Thread Diego Rubolini
Dear all, I'm new to R and to the list, and I have a problem which I'm unable to solve. Consider the following simple simulated data frame:

Re: [R] Label using equivalent of \mathbb{R}

2004-08-26 Thread Simon Cullen
On Thu, 26 Aug 2004 16:40:09 +0200, Trenkler, Dietrich [EMAIL PROTECTED] wrote: [Dietrich Trenkler] Oh yes, I see. I was too overhasty. Maybe using psfrag http://www.ctan.org/tex-archive/macros/latex/contrib/psfrag/ is a way out. Thanks! That is exactly what I

[R] How to change pch in plot.factor ?

2004-08-26 Thread Ernesto Jardim
Hi, I want to change the symbol in plot.factor but using pch is not working. It keeps plotting the horizontal lines ... How can I do it ? Thanks EJ __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

Re: [R] Re:how to compute a mode of a distribution

2004-08-26 Thread John Wilkinson
you are quite right ,I forgot to insert $x, I should have written median(density(my.vector)$x) that gives a median result but it is still not the proper result. density(my.vector)$y) gives the probabilites y of each x value The correst answer would be the value of x corresponding to

[R] GLMM

2004-08-26 Thread Bossarte, Robert
I am trying to use the LME package to run a multilevel logistic model using the following code: --- Model1 = GLMM(WEAP ~ TSRAT2 , random = ~1 | GROUP , family = binomial,

Re: [R] Adding labels to variables

2004-08-26 Thread Mike Prager
At 04:36 PM 08/25/2004, Peter Dalgaard wrote: There's a good chance that R-2.x will introduce usePackage() and eventually remove library() as a tool for loading packages. That would be a welcome instance of consistency. -- Michael Prager, Ph.D.[EMAIL PROTECTED] NOAA Beaufort

[R] plot.new?

2004-08-26 Thread Filippo Biscarini
Hi all, I am just beginning to use R and encountered a problem in trying to draw a regression line using the command abline(coef(var) and similar ones but got the error plot.new has not been called yet. Does anybody know how can I call plot.new? Is this an additional package to be found

Re: [R] o bty with different axes

2004-08-26 Thread bavorak
I don't see why using different scalings is related to bty settings??? Really, after shifting the box type from u to o I received an error message that it can't plot with different scaling on the left and right ordinate. How to fix it? Tomas Bayer __

Re: [R] o bty with different axes

2004-08-26 Thread Uwe Ligges
[EMAIL PROTECTED] wrote: I don't see why using different scalings is related to bty settings??? Really, after shifting the box type from u to o I received an error message that it can't plot with different scaling on the left and right ordinate. How to fix it? Tomas Bayer Please specify an

Re: [R] plot.new?

2004-08-26 Thread Jan Goebel
From help(abline): Examples: data(cars) z - lm(dist ~ speed, data = cars) plot(cars) abline(z) You first have to create a plot to which you can add a line. jan On Thu, 26 Aug 2004, Filippo Biscarini wrote: Hi all, I am just beginning to use R and encountered a

Re: [R] How to change pch in plot.factor ?

2004-08-26 Thread Uwe Ligges
Ernesto Jardim wrote: Hi, I want to change the symbol in plot.factor but using pch is not working. It keeps plotting the horizontal lines ... Have you read ?plot.factor. What kind of a symbol do you expect in a barplot or boxplot (instead of lines)? It would not work for boxplot and barplot

Re: [R] plot.new?

2004-08-26 Thread Uwe Ligges
Filippo Biscarini wrote: Hi all, I am just beginning to use R and encountered a problem in trying to draw a regression line using the command abline(coef(var) and similar ones but got the error plot.new has not been called yet. Does anybody know how can I call plot.new? Is this an additional

Re: [R] plot.new?

2004-08-26 Thread Marta Rufino
Hello, Probably the only probem is that you have to do the plot first, than you put the line in. so, first: plot(x,y) # then you do the abline command. abline(coef..., col=2) Marta __ [EMAIL PROTECTED] mailing list

[R] newbie question: how to read a file into a matrix

2004-08-26 Thread Gaby Aguilera
I am trying to read a file with 8 columns and about 200 rows into a matrix. I tried something like: x-matrix(data=data.txt, nrow=200, ncol=8) but this produces a matrix with data.txt in it. Is there a way to read the data in? Any pointers in the right direction would be appreciated. Thank you!

Re: [R] How to change pch in plot.factor ?

2004-08-26 Thread Ernesto Jardim
On Thu, 2004-08-26 at 16:42, Uwe Ligges wrote: Ernesto Jardim wrote: Hi, I want to change the symbol in plot.factor but using pch is not working. It keeps plotting the horizontal lines ... Have you read ?plot.factor. What kind of a symbol do you expect in a barplot or boxplot

Re: [R] Modalwert

2004-08-26 Thread Thomas Petzoldt
Sonja Dornieden wrote: Hai - kann mir jemand sagen, wie ich den Modalwert in R berechne?! IRgendwie finde ich den Befehl nicht greetz und herzlichen Dank Sonja Hi Sonja, this question seems to appear in regular intervals and you find something about it in the R-Help archives, e.g.:

Re: [R] newbie question: how to read a file into a matrix

2004-08-26 Thread Sundar Dorai-Raj
Gaby Aguilera wrote: I am trying to read a file with 8 columns and about 200 rows into a matrix. I tried something like: x-matrix(data=data.txt, nrow=200, ncol=8) but this produces a matrix with data.txt in it. Is there a way to read the data in? Any pointers in the right direction would be

Re: [R] newbie question: how to read a file into a matrix

2004-08-26 Thread Ernesto Jardim
On Thu, 2004-08-26 at 16:47, Gaby Aguilera wrote: I am trying to read a file with 8 columns and about 200 rows into a matrix. I tried something like: x-matrix(data=data.txt, nrow=200, ncol=8) but this produces a matrix with data.txt in it. Is there a way to read the data in? Any pointers

RE: [R] newbie question: how to read a file into a matrix

2004-08-26 Thread Rau, Roland
Hi, -Original Message- From: [EMAIL PROTECTED] [SMTP:[EMAIL PROTECTED] On Behalf Of Gaby Aguilera Sent: Thursday, August 26, 2004 5:47 PM To: [EMAIL PROTECTED] Subject:[R] newbie question: how to read a file into a matrix matrix

Re: [R] newbie question: how to read a file into a matrix

2004-08-26 Thread Arne Henningsen
Please read the R Data Import/Export manual. It should be on your harddisk and it is also available from http://cran.r-project.org/doc/manuals/R-data.pdf Arne On Thursday 26 August 2004 17:47, Gaby Aguilera wrote: I am trying to read a file with 8 columns and about 200 rows into a matrix.

Re: [R] Label using equivalent of \mathbb{R}

2004-08-26 Thread Martin Maechler
Simon == Simon Cullen [EMAIL PROTECTED] on Thu, 26 Aug 2004 15:30:06 +0100 writes: Simon On Thu, 26 Aug 2004 15:24:33 +0200, Trenkler, Simon Dietrich [EMAIL PROTECTED] Simon wrote: [Dietrich Trenkler]

[R] gls: Newton-Raphson or EM?

2004-08-26 Thread Brendan A. Klick
Hello, Does anyone know whether the gls function in the nlme library uses the Newton-Raphson or EM algorithm to find the restricted log-likelihood or maximum log-likelihood estimates? Brendan Klick [EMAIL PROTECTED] [[alternative HTML version deleted]]

[R] plot.lm basic question

2004-08-26 Thread christian_mora
Hi all, Does anyone know how can I change the color of the 'captions' and 'point labels' in the plots obtained from plot.lm? Thanks for any hint Christian __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] gls: Newton-Raphson or EM?

2004-08-26 Thread Dimitris Rizopoulos
Hi Brendan, according to @Book{pinheiro.bates:00, author= {J. Pinheiro and D. Bates}, title = {Mixed-Effects Models in S and S-PLUS}, year = {2000}, address = {New York}, publisher = {Springer-Verlag} } Section 2.2.8, the optimization procedure for lme (so I suspect

Re: [R] Label using equivalent of \mathbb{R}

2004-08-26 Thread Deepayan Sarkar
On Thursday 26 August 2004 11:41, Martin Maechler wrote: Simon == Simon Cullen [EMAIL PROTECTED] on Thu, 26 Aug 2004 15:30:06 +0100 writes: Simon On Thu, 26 Aug 2004 15:24:33 +0200, Trenkler, Simon Dietrich [EMAIL PROTECTED] Simon wrote: [Dietrich Trenkler]

[R] Replying to digested R-help messages -- get MIME digests!

2004-08-26 Thread Martin Maechler
Hi Frank, Frank == Frank E Harrell [EMAIL PROTECTED] on Thu, 26 Aug 2004 09:08:39 -0500 writes: Frank Sorry about non-indenting of quoted text - I haven't Frank figured out a good way to answer mail when receiving Frank messages in digest mode -FH use a mail client that can

[R] Why terms are dropping out of an lm() model

2004-08-26 Thread John Pitney
Hi all! I'm fairly new to R and not too experienced with regression. Because of one or both of those traits, I'm not seeing why some terms are being dropped from my model when doing a regression using lm(). I am trying to do a regression on some experimental data d, which has two numeric

Re: [R] Why terms are dropping out of an lm() model

2004-08-26 Thread Sundar Dorai-Raj
John Pitney wrote: Hi all! I'm fairly new to R and not too experienced with regression. Because of one or both of those traits, I'm not seeing why some terms are being dropped from my model when doing a regression using lm(). I am trying to do a regression on some experimental data d, which has

RE: [R] EM norm package (NA/NaN/Inf in foreign function call (ar

2004-08-26 Thread Ted Harding
On 26-Aug-04 [EMAIL PROTECTED] wrote: The following code should replicate the error by downloading the dataset from the internet (it is not too big): library(norm) df-download.file(http://www.tariqkhan.org/R/DataFromExcel.csv;, C:/Program Files/R/d.csv) mat-as.matrix(read.table(C:/Program

  1   2   >