[R] Any package to perform HLM (PROC GENMOD) like logistic regression in R?

2005-10-21 Thread William M. Grove
I know there are very nice facilities in Pinhiero and Bates for doing HLM-type modeling for continuous dependent variables. But I would like to be able to do repeated measures logistic regression, or LR on clustered observations (at least with exchangeable correlation structures, and

Re: [R] suppress messages on loading its

2005-10-21 Thread Uwe Ligges
Prof Brian Ripley wrote: Those messages mean it has not been installed correctly -- they should occur during installation only. They are written to the message connection, which capture.output does not mess with (but sink can). On Thu, 20 Oct 2005, Omar Lakkis wrote: On loading the

Re: [R] Any package to perform HLM (PROC GENMOD) like logistic regression in R?

2005-10-21 Thread Prof Brian Ripley
On Fri, 21 Oct 2005, William M. Grove wrote: I know there are very nice facilities in Pinhiero and Bates for doing HLM-type modeling for continuous dependent variables. But I would like to be able to do repeated measures logistic regression, or LR on clustered observations (at least with

Re: [R] Subsetting a list

2005-10-21 Thread Martin Maechler
PaCo == Patrick Connolly [EMAIL PROTECTED] on Fri, 21 Oct 2005 13:26:08 +1300 writes: PaCo On Wed, 19-Oct-2005 at 05:09PM +0200, Martin Maechler wrote: PaCo | Lists can have 'dim' attributes and hence be treated as arrays; PaCo | Note that this is pretty rarely used and not too

Re: [R] problème d'import de fichier

2005-10-21 Thread Yves Magliulo
hello!! si tu veux plus de réponse, you should speak in english here. mais ok je vais essayer de te répondre vu ton erreur, tu devrais regarder la ligne 15 qui doit ne pas avoir le même nombre de séparateur que tes autres lignes. de façon générale, quand tu utilises read.table, assure toi que,

Re: [R] R package that Depends on methods

2005-10-21 Thread Martin Maechler
Thomas == Thomas Schönhoff [EMAIL PROTECTED] on Fri, 21 Oct 2005 00:27:33 +0200 writes: Thomas Hello, Thomas 2005/10/20, Daya Atapattu [EMAIL PROTECTED]: I am trying to build a package that depends on methods package. I am getting an error protection stack overflow

Re: [R] Filter design in R?

2005-10-21 Thread Martin Maechler
tom == tom wright [EMAIL PROTECTED] on Thu, 20 Oct 2005 09:26:05 -0400 writes: tom On Wed, 2005-19-10 at 17:49 -0400, Israel Christie wrote: Dr. Williams, I ran across your inquiry on one of the R-help mailing lists regarding digital filter design and implementation. I

[R] combinging plots

2005-10-21 Thread v . schlecht
If I have one set of points Line 1: {(0,0.87),(0.1,0.88),(0.2,0.89)} and another set of points {(0,0.75),(0.1,0.76),(0.2,0.77)} I can easily produce two separate diagrams: x-c(0,0.1,0.2) y1-c(0.87,0.88,0.89) y2-c(0.75,0.76,0.77) p1-plot(x,y1,type=b) p2-plot(x,y2,type=b) But what can I do if I

Re: [R] forrest plot

2005-10-21 Thread Michela Ballardini
Dear Thomas, thank you very much for the informations, I want just that one! But how can I do? Can you tell me where can I try the commands to do that graphic? Thank you a lot Michela ** Dr.ssa Michela Ballardini Unità di Biostatistica e Sperimentazioni

[R] partial rank correlation coefficient

2005-10-21 Thread Federico Calboli
Hi All, a colleague asked me if R has a function producing a Partial Rank Correlation Coefficient, sensu Blower and Dowlatabadi 1994 [1]. I personally would not have a clue, and I could not find something like that on the search page... although I would not be surprised if it's there under a

Re: [R] combinging plots

2005-10-21 Thread Petr Pikal
Hi from ?plot help page you can find see also ?points ?lines If you have common x as you suggest by your example, matplot will make the desired result for you matplot(x, cbind(y1,y2)) see ?matplot for details HTH Petr On 21 Oct 2005 at 11:59, [EMAIL PROTECTED] wrote: Date sent:

[R] Graph layout software

2005-10-21 Thread Sean Davis
This is a bit off-topic, but I thought I would ask. What are folks typically using for graph layout when drawing graph structures. I have used Rgraphviz for relatively small graphs (100 nodes or fewer) with nice results. However, quick google searches turn up a number of graph layout packages

Re: [R] Subsetting a list

2005-10-21 Thread Prof Brian Ripley
On Fri, 21 Oct 2005, Martin Maechler wrote: PaCo == Patrick Connolly [EMAIL PROTECTED] on Fri, 21 Oct 2005 13:26:08 +1300 writes: PaCo On Wed, 19-Oct-2005 at 05:09PM +0200, Martin Maechler wrote: PaCo | Lists can have 'dim' attributes and hence be treated as arrays; PaCo | Note

[R] [R-pkgs] new package bundle: onion

2005-10-21 Thread Robin Hankin
Dear List I have just uploaded a new package, onion, to CRAN. It provides some functionality for manipulating and visualizing quaternions and octonions. All comments welcome! enjoy rksh -- Robin Hankin Uncertainty Analyst National Oceanography Centre, Southampton European Way,

[R] Finding the neighbors of the point

2005-10-21 Thread Leaf Sun
Dear all, I got point data of trees. I was wondering if anybody has experience in searching the neighbors within a specified distance efficiently. XY Z 99 34 65 98 35 29 98 34 28 99 33 33 98 32 23 99 33 21 99 33

[R] please recommend

2005-10-21 Thread Weiwei Shi
Dear listers: These days I am analyzing some microarray data and I am wondering if there are some good books and references on the following topics: R, data mining and microarray analysis. I appolagize if this is a bit off topic.. Thanks for the recommendation! -- Weiwei Shi, Ph.D Did you

[R] change maxiter for nls

2005-10-21 Thread Xiaodong Jin
I typed the following commands but it still use maxiter=50 after the 2nd command: nls.control(maxiter = 1000) nls(..) Thanks! - [[alternative HTML version deleted]] __

Re: [R] [BioC] please recommend

2005-10-21 Thread James MacDonald
Here is an excellent choice: http://www.amazon.com/exec/obidos/tg/detail/-/0387251464/qid=1129910139/sr=8-1/ref=pd_bbs_1/002-2566394-3995236?v=glances=booksn=507846 Best, Jim James W. MacDonald Affymetrix and cDNA Microarray Core University of Michigan Cancer Center 1500 E. Medical Center

Re: [R] change maxiter for nls

2005-10-21 Thread Peter Dalgaard
Xiaodong Jin [EMAIL PROTECTED] writes: I typed the following commands but it still use maxiter=50 after the 2nd command: nls.control(maxiter = 1000) nls(..) Thanks! You need nlc - nls.control(maxiter = 1000) nls(..., control=nlc) And yes, I suppose both help pages could be

Re: [R] [BioC] please recommend

2005-10-21 Thread Berton Gunter
and this: http://www.amazon.com/exec/obidos/tg/detail/-/1584883278/qid=1129910385/sr=2 -1/ref=pd_bbs_b_2_1/002-5914105-6054455?v=glances=books -- Bert Gunter Genentech Non-Clinical Statistics South San Francisco, CA The business of the statistician is to catalyze the scientific learning

[R] using markstat

2005-10-21 Thread Sara Mouro
Dear all, I am using the spatstat package, and in particular the algorithm markstat(X,fun,N,R,.) I need that output to detect which is the nearest neighbour of each point in that point pattern. Therefore, I would use X=my ppp object already defined, N=2 But what function should I use? It

[R] read data from pdf file

2005-10-21 Thread Marco Venanzi
Hi, I'm trying to read data from a PDF file.Is it possible to do it with R? Thanks, Marco [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] order of approximation used in glmmPQL

2005-10-21 Thread Edmond Ng
Dear all, Could anyone let me know what order of approximation (first, second or higher order) is used in glmmPQL in the MASS library please? I have read the section on this wrapper command in Venables and Ripley but am not sure if the relevant information is there. Many thanks in advance.

[R] make three plot to one plot

2005-10-21 Thread Jan Sabee
Dear all, I want to make three plot below to only one plot together with legend, how can I do that? I have tried with matplot function but I did not succeed. Thanks for your help. Sincerelly, Jan Sabee test.five.x - c(0.02,0.05,0.07,0.09,0.10,0.12,0.13,0.14,0.16,0.17,0.20,0.21,0.34,0.40)

Re: [R] read data from pdf file

2005-10-21 Thread Ted Harding
On 21-Oct-05 Marco Venanzi wrote: Hi, I'm trying to read data from a PDF file.Is it possible to do it with R? Thanks, Marco Basically, No. But you may be lucky with copypaste using the mouse, from the display generated in Acrobat Reader to a text file. The basic procedure here is 1. Click

[R] locator

2005-10-21 Thread tom wright
I'm trying to use the locator function on a drawing area with multiple graphs par(mfrow=c(1,2)) Is it possible to identify which graph has been clicked? thanks tom __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] read data from pdf file

2005-10-21 Thread Thomas Schönhoff
2005/10/21, Ted Harding [EMAIL PROTECTED]: On 21-Oct-05 Marco Venanzi wrote: Hi, I'm trying to read data from a PDF file.Is it possible to do it with R? Thanks, Marco Basically, No. But you may be lucky with copypaste using the mouse, from the display generated in Acrobat Reader to a

Re: [R] locator

2005-10-21 Thread Duncan Murdoch
On 10/21/2005 11:23 AM, tom wright wrote: I'm trying to use the locator function on a drawing area with multiple graphs par(mfrow=c(1,2)) Is it possible to identify which graph has been clicked? locator() will return coordinates based on the active graph (typically the last one you drew), so

Re: [R] read data from pdf file

2005-10-21 Thread Thomas Schönhoff
Hello again, 2005/10/21, Thomas Schönhoff [EMAIL PROTECTED]: 2005/10/21, Ted Harding [EMAIL PROTECTED]: On 21-Oct-05 Marco Venanzi wrote: Hi, I'm trying to read data from a PDF file.Is it possible to do it with R? Thanks, Marco Hmm, if this doesn't work you should have a look to

Re: [R] make three plot to one plot

2005-10-21 Thread Marc Schwartz (via MN)
On Fri, 2005-10-21 at 21:13 +0200, Jan Sabee wrote: Dear all, I want to make three plot below to only one plot together with legend, how can I do that? I have tried with matplot function but I did not succeed. Thanks for your help. Sincerelly, Jan Sabee test.five.x -

Re: [R] locator

2005-10-21 Thread tom wright
Perfect, thankyou On Fri, 2005-21-10 at 15:46 -0400, Duncan Murdoch wrote: On 10/21/2005 11:23 AM, tom wright wrote: I'm trying to use the locator function on a drawing area with multiple graphs par(mfrow=c(1,2)) Is it possible to identify which graph has been clicked? locator() will

Re: [R] Finding the neighbors of the point

2005-10-21 Thread Roger Bivand
On Fri, 21 Oct 2005, Leaf Sun wrote: Dear all, I got point data of trees. I was wondering if anybody has experience in searching the neighbors within a specified distance efficiently. XY Z 9934 65 9835 29 9834 28 9933 33 9832

Re: [R] make three plot to one plot

2005-10-21 Thread Gabor Grothendieck
# we can create a zoo object library(zoo) z - merge(five = zoo(test.five.y, test.five.x), six = zoo(test.six.y, test.six.x), seven = zoo(test.seven.y, test.seven.x)) # and then plot it all in one go using na.approx to fill in generated NAs plot(na.approx(z), plot.type = single,

[R] Help with SOM clustering package

2005-10-21 Thread Nayeem Quayum
Hello everybody, I am trying to use SOM package to cluster some gene expression data and there fore requesting you some help. I would really appreciate if you cold help me out. Q1. I clustered the data using SOM function and plotted it using plot function. But I can't figure out how I can list

[R] peculiar matrices

2005-10-21 Thread Ben Bolker
As far as I can tell from reading The Fine Documentation (R Language Definition and Intro to R), matrices are supposed to be of homogeneous types. Yet giving matrix() an inhomogeneous list seems to work, although it produces a peculiar object: v = list(1:3,4,5,a) m = matrix(v,nrow=2) m

Re: [R] read data from pdf file

2005-10-21 Thread Thomas Schönhoff
Hi, 2005/10/21, Thomas Schönhoff [EMAIL PROTECTED]: Hello again, 2005/10/21, Thomas Schönhoff [EMAIL PROTECTED]: 2005/10/21, Ted Harding [EMAIL PROTECTED]: On 21-Oct-05 Marco Venanzi wrote: Hi, I'm trying to read data from a PDF file.Is it possible to do it with R? Thanks, Marco

Re: [R] make three plot to one plot

2005-10-21 Thread Jan Sabee
Thanks to Marc and Gabor. Have a nice weekend. Best, Jan Sabee On 10/21/05, Jan Sabee [EMAIL PROTECTED] wrote: Dear all, I want to make three plot below to only one plot together with legend, how can I do that? I have tried with matplot function but I did not succeed. Thanks for your help.

Re: [R] peculiar matrices

2005-10-21 Thread Marc Schwartz (via MN)
On Fri, 2005-10-21 at 21:32 +, Ben Bolker wrote: As far as I can tell from reading The Fine Documentation (R Language Definition and Intro to R), matrices are supposed to be of homogeneous types. Yet giving matrix() an inhomogeneous list seems to work, although it produces a peculiar

Re: [R] peculiar matrices

2005-10-21 Thread Berton Gunter
matrix coerces it's arguments to a single mode first. Try mode(m) -- iit's a list (with a dim attribute). If you print each entry of m separately, you'll find they're all lists. As the docs say, if either of nrow or ncol aren't given, it tries to infer what they should be from the data and other

[R] R Programmer/Analyst

2005-10-21 Thread Earl F. Glynn
From http://www.stowers-institute.org/ScientistsSought/ScientistsSought.asp#positions Programmer/Analyst The Stowers Institute for Medical Research has an opening for a Programmer/Analyst to support scientific data analysis and assist with computational biology tasks. Responsibilities

Re: [R] read data from pdf file

2005-10-21 Thread Dave Roberts
In linux (and possibly other *nixes) you can view the file with xpdf and simply cut and paste it into another window (I use vi) and it's converted to ASCII text on the fly. For large documents you might have to scroll quite a bit to convert the whole document, but this process has saved my

Re: [R] Any package to perform HLM (PROC GENMOD) like logistic

2005-10-21 Thread Alan Cobo-Lewis
I think you're looking for lmer( ... family=binomial ) in package lme4 by Bates Sarkar. William M. Grove [EMAIL PROTECTED] wrote: I know there are very nice facilities in Pinhiero and Bates for doing HLM-type modeling for continuous dependent variables. But I would like to be able to do

[R] Generalised rbind/cbind

2005-10-21 Thread hadley wickham
Dear list, Is there a generalised form of rbind/cbind for combining matrices/arrays into higher-D structures? ie. if I have: a - matrix(2,2,2) b - matrix(3,2,2) how can I get array(rep(c(3,2), each=4), c(2,2,2)) ? It seems like this would be the job of a generalised abind function:

Re: [R] Generalised rbind/cbind

2005-10-21 Thread Marc Schwartz
On Fri, 2005-10-21 at 17:42 -0500, hadley wickham wrote: Dear list, Is there a generalised form of rbind/cbind for combining matrices/arrays into higher-D structures? ie. if I have: a - matrix(2,2,2) b - matrix(3,2,2) how can I get array(rep(c(3,2), each=4), c(2,2,2)) ? It

Re: [R] random coefficient multinomial logistic regression (was: request)

2005-10-21 Thread Spencer Graves
I haven't seen a reply to this post, so I just experimented with a few terms in R Site Search. RSiteSearch(multinomial hierarchical model) produced 8 hits, the first of which was http://finzi.psych.upenn.edu/R/library/bayesm/html/rhierMnlRwMixture.html;. hope this helps.

Re: [R] npmc package

2005-10-21 Thread Kjetil Holuerson
Martin Maechler wrote: Carlos == Carlos Mauricio Cardeal Mendes [EMAIL PROTECTED] on Wed, 19 Oct 2005 15:11:32 -0300 writes: Carlos So, is there another package to substitute those Carlos functions described on ORPHANED npmc package ? May be not. But nobody stops you from

Re: [R] npmc package

2005-10-21 Thread Kjetil Holuerson
Martin Maechler wrote: Carlos == Carlos Mauricio Cardeal Mendes [EMAIL PROTECTED] on Wed, 19 Oct 2005 15:11:32 -0300 writes: Carlos So, is there another package to substitute those Carlos functions described on ORPHANED npmc package ? May be not. But nobody stops you from

Re: [R] peculiar matrices

2005-10-21 Thread Gabor Grothendieck
The reference manual of 2.2.0 says in section 2.2 that Matrices and arrays are simply vectors with the attribute dim and optionally dimnames. Now earlier in section 2.1 it discusses vectors and I think that that is where the confusing part lies. Section 2.1 starts out saying that Vectors can be

[R] Sorting Numeric and Character Data

2005-10-21 Thread Thomas Hopper
Hello, I have what seems like an easy question to answer, but I'm struggling with it. I have a set of categorical data that I am reading in, looking something like: category result A .234 B .123 C .564 D -.452 E .112 F -.106 I'd like to plot this twice on two separate dot charts, once with

Re: [R] Sorting Numeric and Character Data

2005-10-21 Thread Gabor Grothendieck
Use read.table(myfile, header = TRUE, as.is = TRUE) where as.is=TRUE causes read.table not to convert character data to factors. On 10/21/05, Thomas Hopper [EMAIL PROTECTED] wrote: Hello, I have what seems like an easy question to answer, but I'm struggling with it. I have a set of