Re: [R] R version conflict.

2003-03-04 Thread Peter Dalgaard BSA
Bai Yan [EMAIL PROTECTED] writes:

 My mistake :(
 Seems that there's another R installed by root. So when I use .libPaths()
 there're two of them:
 
  .libPaths()
 [1] /disk/hopper/projects/class/cse514/R/library
 [2] /usr/lib/R/library
 
 How could the second path be enabled then? The R_LIBS have already been
 set to as the [1] only.

If that happens, you're running the installed R, not the one you
compiled. You need to run it as /where/you/put/it/bin/R (or bin/R if
you're in the toplevel build dir).

-- 
   O__   Peter Dalgaard Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics 2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark  Ph: (+45) 35327918
~~ - ([EMAIL PROTECTED]) FAX: (+45) 35327907

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] writing several command line in R console

2003-03-04 Thread Petr Pikal
Hallo

On 4 Mar 2003 at 2:36, Vincent Stoliaroff wrote:

 
 Hi R lovers
 
 I would like to know how to step to the next line in the R console
 editor without breaking the continuity of my code more clearly : if
 for example I write a function, so far i have to write the all code
 inside on the same line wich may become obscure as the function is
 more and more complex. I would like to do like in the example of the
 manuels:
 
 twosam - function(y1, y2) {
 n1  - length(y1); n2  - length(y2)
 yb1 - mean(y1);   yb2 - mean(y2)s1  - var(y1);s2  -
 var(y2) s - ((n1-1)*s1 + (n2-1)*s2)/(n1+n2-2) tst - (yb1 -
 yb2)/sqrt(s2*(1/n1 + 1/n2))tst  }
 
 
 all I can do is something like that:
 
 twosam - function(y1, y2) {n1  - length(y1); n2  - length(y2)
 
 +yb1 - mean(y1);   yb2 - mean(y2)s1  - var(y1);s2  -
 +var(y2) s - ((n1-1)*s1 + (n2-1)*s2)/(n1+n2-2)
 +   tst - (yb1 - yb2)/sqrt(s2*(1/n1 + 1/n2))tst  }
 
 with the sign + in front of each line
 What does this sign mean? and how could I solve my problems
 Thanks

Maybe looking to the manuals will be more precise than my 
answer.

+ is a continuity sign and means you did not finished your imput 
and you shall continue typing your command.

But why you do not use any text ditor for writing functions and 
than copy/paste to R command window?




 
 __
 [EMAIL PROTECTED] mailing list
 http://www.stat.math.ethz.ch/mailman/listinfo/r-help

CheersPetr Pikal
[EMAIL PROTECTED]
[EMAIL PROTECTED]

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


RE: [R] Embed R in other applications

2003-03-04 Thread Pfaff, Bernhard
see the contributed pacakge R2HTML; in case you have access to a databank:
RODBC package and the like might also be an option in combination with
*.php and/or *.asp-files for web-publishing.
I have implemented both and it runs succesfully without any hickups.
Updating the data, running R-program as batch via scheduler and outputting
directly or writing it first into a databank for further processing.

HTH,
Bernhard

-Original Message-
From: Bai Yan [mailto:[EMAIL PROTECTED]
Sent: 04 March 2003 06:01
To: [EMAIL PROTECTED]
Subject: [R] Embed R in other applications


Hello all,

I want to create webpages which can display the R results. The first step
is essentially getting input from client, generating graphics and
returning a page containing those images.

What should I do for this? Are there some sources on the web? Any hints
will be very helpful.

Thanks very much :)

Yan

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


--
If you have received this e-mail in error or wish to read our e-mail 
disclaimer statement and monitoring policy, please refer to 
http://www.drkw.com/disc/email/ or contact the sender.
--

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] How to estimate or approximate a 3-D surface

2003-03-04 Thread Petr Pikal
Hi

On 3 Mar 2003 at 22:03, Feng Zhang wrote:

 Hey, R-listers
 
 Now I am going to estimate or approximate a surface in
 3-D space given a large enough number of (x,y,z) data sets.
 
 So for these 3-D data points, is it possible to get a surface
 function, like z=f(x,y) to represent this underlying surface?

I am not sure but interp from akima package shall give you an 
interpolation of your data.
 
 Thanks for your time and point.
 
 Fred
 
 __
 [EMAIL PROTECTED] mailing list
 http://www.stat.math.ethz.ch/mailman/listinfo/r-help

Cheers

Petr Pikal
[EMAIL PROTECTED]
[EMAIL PROTECTED]

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Recall: [R] Embed R in other applications

2003-03-04 Thread Pfaff, Bernhard
Pfaff, Bernhard would like to recall the message, [R] Embed R in other
applications.


--
If you have received this e-mail in error or wish to read our e-mail 
disclaimer statement and monitoring policy, please refer to 
http://www.drkw.com/disc/email/ or contact the sender.
--

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] R in your pocket on a Sharp Zaurus

2003-03-04 Thread Stuart Leask
R v.1.6.2 is now available for the Sharp SL-5500 Zaurus linux pda. Thanks to
the effort of Simon Pickering at Bath University, NA handling is now correct
and graphics function fine. It can run on the unexpanded 'out of the box'
machine. The FPU emulation can slow some operations a lot, but it can do
lowess smoothing on over 10,000 data points in just a few minutes, as good
as the 486SX25 I was happy with not that many years ago. And R (+Zaurus) in
the hand/jacket pocket is worth any number of laptops in the bush.

You'll need R (see http://students.bath.ac.uk/enpsgp/Zaurus/#R), plus a few
libraries (from http://ipkgfind.handhelds.org/). To add graphics, you can
either run native X-windows - faster - or run X within qtopia (see
http://students.bath.ac.uk/enpsgp/Zaurus/#X) for convenience.

Stuart

PS.
At present the help system doesn't function. We're working on this - we're
not sure why, but note that in the current compilation, the CONTENTS folders
end up empty. All suggestions gratefully received...

Dr Stuart Leask MA MRCPsych, Clinical Lecturer in Psychiatry
University of Nottingham Dept of Psychiatry, Duncan Macmillan House
Porchester Road, Nottingham. NG3 6AA. UK
http://www.nottingham.ac.uk/psychiatry/staff/sjl.html

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] tex/pdf?

2003-03-04 Thread Kurt Sys
Hello,

I'm having some problems, now with building libraries/packages. Actually, with making 
pdf- or latex-files. So, everyting is ok, I can generate latex-files ('.tex') using 'R 
CMD Rdconv -t tex', but apparently, it are not valid latex-files (no 'documentclass', 
no 'body' etc). How come? Do I have to add these things manually?


Kurt.

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


RE: [R] log axis assignment

2003-03-04 Thread Morten Sickel
From: Jeremy Z Butler [mailto:[EMAIL PROTECTED]
 Surely I should be able to set the 
 ylog option to T using another object (logaxis).

for (n in colnames(raw))
{
if(n==pH){logaxis-T} else {logaxis-F}
plot(full.age,raw[,n],type=n,ylog=logaxis)
...
}

What (probably) will work is:

if(n==pH){logaxis-Y} else {logaxis-}
plot(full.age,raw[,n],type=n,log=logaxis)

The plot-argument log takes which axis that shuld be logaritmic. 


Morten

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] tex/pdf?

2003-03-04 Thread ripley
On Tue, 4 Mar 2003, Kurt Sys wrote:

 I'm having some problems, now with building libraries/packages. 
Actually, with making pdf- or latex-files. So, everyting is ok, 
I can generate latex-files ('.tex') using 'R CMD Rdconv -t tex', 
but apparently, it are not valid latex-files (no 'documentclass', 
no 'body' etc). How come? Do I have to add these things manually?

No, you use help(offline=TRUE) or R CMD Rd2dvi (which is documented in 
R-exts).


Please wrap your lines in emails 

-- 
Brian D. Ripley,  [EMAIL PROTECTED]
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel:  +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UKFax:  +44 1865 272595

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] How to estimate or approximate a 3-D surface

2003-03-04 Thread Simon Wood
Have you tried looking at packages mgcv or gss? 

 Now I am going to estimate or approximate a surface in
 3-D space given a large enough number of (x,y,z) data sets.
 
 So for these 3-D data points, is it possible to get a surface
 function, like z=f(x,y) to represent this underlying surface?

Simon
_
 Simon Wood [EMAIL PROTECTED]www.stats.gla.ac.uk/~simon/
  Department of Statistics, University of Glasgow, Glasgow, G12 8QQ
   Direct telephone: (0)141 330 4530  Fax: (0)141 330 4814

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] How to estimate or approximate a 3-D surface

2003-03-04 Thread Stuart Leask
Locfit has worked well for me.
http://cm.bell-labs.com/cm/ms/departments/sia/project/locfit/

Stuart

Dr Stuart Leask MA MRCPsych, Clinical Lecturer in Psychiatry
University of Nottingham Dept of Psychiatry, Duncan Macmillan House
Porchester Road, Nottingham. NG3 6AA. UK
http://www.nottingham.ac.uk/psychiatry/staff/sjl.html

- Original Message - 
From: Simon Wood [EMAIL PROTECTED]
To: Feng Zhang [EMAIL PROTECTED]
Cc: R-Help [EMAIL PROTECTED]
Sent: Tuesday, March 04, 2003 9:55 AM
Subject: Re: [R] How to estimate or approximate a 3-D surface


 Have you tried looking at packages mgcv or gss? 
 
  Now I am going to estimate or approximate a surface in
  3-D space given a large enough number of (x,y,z) data sets.
  
  So for these 3-D data points, is it possible to get a surface
  function, like z=f(x,y) to represent this underlying surface?
 
 Simon
 _
  Simon Wood [EMAIL PROTECTED]www.stats.gla.ac.uk/~simon/
   Department of Statistics, University of Glasgow, Glasgow, G12 8QQ
Direct telephone: (0)141 330 4530  Fax: (0)141 330 4814
 
 __
 [EMAIL PROTECTED] mailing list
 http://www.stat.math.ethz.ch/mailman/listinfo/r-help

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] writing several command line in R console

2003-03-04 Thread Peter Dalgaard BSA
Petr Pikal [EMAIL PROTECTED] writes:


  I would like to know how to step to the next line in the R console
  editor without breaking the continuity of my code more clearly : if
  for example I write a function, so far i have to write the all code
  inside on the same line wich may become obscure as the function is
  more and more complex. I would like to do like in the example of the
  manuels:
  
  twosam - function(y1, y2) {
  n1  - length(y1); n2  - length(y2)
  yb1 - mean(y1);   yb2 - mean(y2)s1  - var(y1);s2  -
  var(y2) s - ((n1-1)*s1 + (n2-1)*s2)/(n1+n2-2) tst - (yb1 -
  yb2)/sqrt(s2*(1/n1 + 1/n2))tst  }

...
 
 + is a continuity sign and means you did not finished your imput 
 and you shall continue typing your command.
 
 But why you do not use any text ditor for writing functions and 
 than copy/paste to R command window?

I'm not sure if and how it works on Windows, but on Unix one of the
better kept secrets of the readline library is that you can embed
newlines by typing ctr-V ctr-J, e.g.

 x - matrix(c(1,2,3,
4,5,6,
7,8,9), 3)

The nice thing is that it allows you to recall the entire command and
edit it (I don't think it survives being saved to the history file
though). If you have ever had to correct a typo in an expression that
has been split over 8 lines, you'll know what I mean:

up,up,up,up,up,up,up,up,RET,
up,up,up,up,up,up,up,up,RET,
up,up,up,up,up,up,up,up,RET,
up,up,up,up,up,up,up,up,fix typo,RET,
up,up,up,up,up,up,up,up,RET,
up,up,up,up,up,up,up,up,RET,
up,up,up,up,up,up,up,up,RET,
up,up,up,up,up,up,up,up,RET

-- 
   O__   Peter Dalgaard Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics 2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark  Ph: (+45) 35327918
~~ - ([EMAIL PROTECTED]) FAX: (+45) 35327907

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] How to estimate or approximate a 3-D surface

2003-03-04 Thread ripley
On Tue, 4 Mar 2003, Stuart Leask wrote:

 Locfit has worked well for me.
 http://cm.bell-labs.com/cm/ms/departments/sia/project/locfit/

Unfortunately it has not been updated for current versions of R: it also 
writes diagnostic output to stderr which is not shown on GUI versions of R
(especially on Windows).  There is a version in the Devel area on CRAN.

Package akima is the most obvious choice for interpolation.  For
smoothing, there are all the various spatial packages (and package fields)
as well as mgcv and gss already mentioned.

-- 
Brian D. Ripley,  [EMAIL PROTECTED]
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel:  +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UKFax:  +44 1865 272595

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] Parallel programming using R

2003-03-04 Thread lun li
Dear All,

For my regression problem(thousand dependents by thousand predictors by 
hundred months observations), I prefer to  parallel computation in R.  Is it 
possible to parallel programming using R as using C and Fortran?

Cheers,

Lun

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] tex/pdf?

2003-03-04 Thread Kurt Sys
Mail from [EMAIL PROTECTED], sent on Tuesday March 4
2003 at 09:42 (GMT+):

   On Tue, 4 Mar 2003, Kurt Sys wrote:
   
I'm having some problems, now with building
   libraries/packages.  Actually, with making pdf- or
   latex-files. So, everyting is ok, I can generate
   latex-files ('.tex') using 'R CMD Rdconv -t tex', but
   apparently, it are not valid latex-files (no
   'documentclass', no 'body' etc). How come? Do I have to
   add these things manually?
   
   No, you use help(offline=TRUE) or R CMD Rd2dvi (which is
   documented in R-exts).
   

Ok... tnx. Was a stupid question. Should have found that myself...

Kurt.

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] Instability in package gld in R?

2003-03-04 Thread Steve Su
Dear All,

Is there anyone using package gld in R?

I have managed to load this package into Splus.

I found that in some cases the use of this package will cause Splus to crash and R to 
hang, I am not sure why but this is just an example:

qgl(c(0,1),45.23587,0.2363868,-1.5,0)

Which means the functon ranges from -Inf to Inf, however if you put:

pgl(150,45.23587,0.2363868,-1.5,0)

Then Splus will crash and R will hang.

There are other examples in dgl that also causes R to hang and Splus to quit. 

I think this may have to do with using numerical methods to approximate the 
distribution function of gld. Does anyone have any ideas about how to fix this problem 
or correct this difficulty?

Steve. 







**

 Steve Su ([EMAIL PROTECTED])
 PhD student.  

 School of Accountancy 
 School of Mathematical Sciences 
 Queensland University of Technology  
 
 Postal Address: Steve Su, School of Accountancy, QUT, PO Box 2434, Brisbane, 
 Queensland, Australia, 4000.   


 Phone:  +61 7 3864 2017
 Fax:+61 7 3864 1812  
 Mobile: 0421  840  586  
 .   
   _--_|\  
  /  QUT   
  \_.--._/   
v   
  
** 
  

[[alternate HTML version deleted]]

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] error using try() and coxph()

2003-03-04 Thread vito muggeo
Dear all,
I'm experiencing some problems in using the following code to perform
simulations with Cox models (using coxph() in survival package):

for(i in 1:1000){
dd-#simulate the dataframe
o0-try(coxph(Surv(start,stop,status)~x,data=dd))
o1-try(coxph(Surv(start,stop,status)~x+x:stop,data=dd))
o2-try(coxph(Surv(start,stop,status)~x+x:poly(stop,2),data=dd))

}

The error (both in Rgui and in batch mode) causes R to crash immediately
with a DOS message like Error in survival.dll
Using coxph() without try() simulation stops because sometimes the model may
be not fitted, I get the explicit error message. Therefore I suspect that
the problem depends on using try(coxph()). Also note that using
try(arima()) for a similar problem works well.
versions: survival 2.9-4 and R 1.5.1 on WinMe
I know I should update R, but I wonder whether the problem depends on the
old version

Could anyone suggest me anything?

cheers,
vito

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] Unit root test for panel data

2003-03-04 Thread Cenktan Ozyildirim

Dear Sir/Madam,
I need to apply IPS, Fischer or any other more contemporary panel data unit
root test. I would be glad if you could inform me whether there are
r-modules which can perform these tests. I thank you very much in advance
for your precious time and interest.
Yours Sincerely,
Cenktan zyldrm

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] How to modify XML documents and save changes

2003-03-04 Thread Duncan Temple Lang
Again, Stephen is right on the mark with his explanation of modifying
the XML objects and needing to put the modified values back into the
containing structure. Trees  graphs are slightly cumbersome in the S
language since it is not a reference-based language.  There are some
tricks one can use that are in the XML package for indexing them with
a cursor, but these are (currently) only for use when constructing a
tree from scratch within S.

 A little cumbersome, but doable. One other option is to use xmlEventParse
 and write a handler that would add the element after the one you're
 interested in.  I hope there is a better way, but haven't seen it yet. :-(

Here there is some good news. xmlEventParse is extremely low-level
relative to xmlTreeParse and typically used for efficiency in
minimizing the amount of storage used when parsing the contents of an
XML document. The idea is that an S function gets invoked when an XML
tag is opened, and another when it is closed.

xmlTreeParse() can also be used in this event-driven style
programming.  When the XML parser encounters the close of a tag, it
looks for a suitable S function to invoke specified in the handlers
argument of xmlTreeParse.  If there is an element in that list of
functions whose name matches the tag name, it is invoked.  Thus if you
programmatically want to augment/modify the contents of an XML node,
you can supply a function that operates on that node.  The function
can return the resulting updated node and it will be inserted into the
overall tree as one would expect.

So, suppose one wants to _always_ add a node 
  Norm/Norm
to each tagname/tagname node, you can use the following

myFun - function(node,...) {
  node - append.XMLNode(node, xmlNode(Norm))
  node
}

and then use this in the command

xmlTreeParse(file.xml, handlers = list(tagname = myFun), asTree = TRUE)

One can also use XSL to do this, either directly from the command
line or from S via the Sxslt package.

If one only wants to modify particular tagname elements, it is
probably easiest to read the entire tree and modify it directly.

For people familiar with other XML parsers, the default xmlTreeParse
behaves like DOM and xmlEventParse behaves likes SAX.  With handlers,
xmlTreeParse() provides a hybrid parser for XML somewhere between DOM
and SAX.

 D.


Stephen C. Upton wrote:
 Steffen,
 
 As with most R objects, you're basically putting a copy of the R object into
 the new object. Any operation or function you apply to that object does not
 affect the original. Same goes for append.xmlNode - you're appending to the
 original and getting back another structure that is the original plus the
 new node. Finally, saveXML works on an object of XMLInternalDocument and doc
 (the object returned by xmlTreeParse) is a XMLDocument object.
 
 Here's one suggestion:
 1. Read in the doc object as you've done, but manipulate the structure as
 you would any other R object, e.g.,
 QTListNode [[1]] -
 append.xmlNode(QTListNode[[1]],xmlNode(name=Norm,attrs=NULL))
 2. you then need to modify that node within the root with this new, modified
 node, e.g. (assume I've assigned root - xmlRoot(doc)),
 root[[whateverindextagnameis]] - QTListNode[[1]]
 (I use the index here rather than the name, since it's unique - if you use a
 name,e.g., root[[tagname]], that just adds a list element to root with
 name tagname)
 3. to write out, use write with toString.XMLNode,
 write(toString.XMLNode(root, file=out.xml)
 
 Since this is probably a little cumbersome, suggest writing a function to do
 the finding, appending, and replacement.
 
 A little cumbersome, but doable. One other option is to use xmlEventParse
 and write a handler that would add the element after the one you're
 interested in.  I hope there is a better way, but haven't seen it yet. :-(
 
 HTH
 steve
 
 
 Steffen Durinck wrote:
 
  Dear,
 
  I want to read XML documents, add child nodes to some elements and store
  everything back as an XML document.
 
  I've tryed the following:
 
  doc - xmlTreeParse(file.xml)
  QTListNode-xmlElementsByTagName(xmlRoot(doc)[[1]],tagname)
  append.xmlNode(QTListNode[[1]],newXMLNode(name =Norm, attrs = NULL))
  saveXML(doc, file = out.xml, compression = 0, indent=T)
 
  This doesn't seem to work.
  Can anyone help?
 
  Thanks,
  Steffen
 
  __
  [EMAIL PROTECTED] mailing list
  http://www.stat.math.ethz.ch/mailman/listinfo/r-help
 
 __
 [EMAIL PROTECTED] mailing list
 http://www.stat.math.ethz.ch/mailman/listinfo/r-help

-- 
___

Duncan Temple Lang[EMAIL PROTECTED]
Bell Labs, Lucent Technologiesoffice: (908)582-3217
700 Mountain Avenue, Room 2C-259  fax:(908)582-3340
Murray Hill, NJ  07974-2070   
 http://cm.bell-labs.com/stat/duncan

__
[EMAIL PROTECTED] mailing list

[R] linear model with arma errors

2003-03-04 Thread Gregory BENMENZER
Dear all,

I'm looking for how can I estimate a linear model with ar(ma) errors :

y(t)=a*X(t)+e(t) with
P(B)e(t)=Q(B)u(t)

where u is a white noise and P, Q are some polynomes.

Could you help me ?

Grégory Benmenzer

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] Wavelets correlation test

2003-03-04 Thread Brandon Whitcher

 Hello, I use wavethresh packages to perform wavelet analysis. In
 particular, I would like to compare 2 signals (vectors) after a wavelet
 decomposition. I would like to use cor.test function, but this function
 acts on the entire vector values. I plan to perform a cor.test on each
 level of the wavelet decomposition, say N. So I will have at the end of
 a first step N results of cor.test.

 How can I deal with this N results to have an answer globaly ?

A lot of this code is already implemented in an alternative wavelet
analysis package for R (waveslim).  There, I have placed code to compute
multiscale covariance/correlation estimates that may also be a function of
lag.  This would provide a correlation coefficient (or autocorrelation
sequence) per scale.  Approximate confidence intervals for these estimates
are also available, there might be a reference in wave.covariance or
wave.correlation -- I can't remember.  The papers of interest are:

- Whitcher, Guttorp, Percival (2000).  Wavelet analysis of covariance with
application to atmospheric time series, JGR 105 (D11), 14,941-14,962.

- Serroukh and Walden (2000).  Wavelet scale analysis of bivariate time
series I: Motivation and estimation, Journal of Nonparametric Statistics,
13 (1), 1-36.

- Gencay, Selcuk, Whitcher (2001).  An Introduction to Wavelets and Other
Filtering Methods in Finance and Economics, Academic Press.

I think it would be difficult to make a statement for the whole series.
Some scales may agree and some scales may not.  What is the scientific
question you are trying to answer?

cheers...

Brandon

--
 Senior Researcher in Imaging GlaxoSmithKline
 Research Statistics Unit  New Frontiers Science Park (South)
  Third Avenue, Harlow   CM19 5AW
 phone:  +44 (0)127 963 1285   United Kingdom
 fax:+44 (0)127 964 4004

[EMAIL PROTECTED]
--

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] linear model with arma errors

2003-03-04 Thread Spencer Graves
That's discussed in Pinhiero and Bates, Mixed Effects Models in S and 
S-Plus.  I don't have the book in my hand now, so I can't tell you 
exactly where to find it.  However, I'm pretty confident that it can be 
done in lme.

Spencer Graves

Gregory BENMENZER wrote:
Dear all,

I'm looking for how can I estimate a linear model with ar(ma) errors :

y(t)=a*X(t)+e(t) with
P(B)e(t)=Q(B)u(t)
where u is a white noise and P, Q are some polynomes.

Could you help me ?

Grégory Benmenzer

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help
__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] suggestion for addition to R-lang.texi

2003-03-04 Thread Peter Dalgaard BSA
Stephen Eglen [EMAIL PROTECTED] writes:

 Small suggestion for the documentation:
 
 Should the = operator for assignments, described on
 http://developer.r-project.org/equalAssign.html, now be included in
 the operator table in doc/manual/R-lang.texi?  If so, one possibility
 could be after lines 1096/7:
 
 @item @code{-}
 @tab Left assignment, binary
 @item @code{=}
 @tab Left assignment, binary
 
 Also, where is the best place to send suggestions of this nature? I
 didn't want to send it to r-bugs, since it is not a bug.

It could easily be considered one... There is a Documentation category
in the bug repository, so why not?

BTW, I wonder: why did your mail come with sender r-help-bounces
instead of the usual r-help-admin and confuse my mail filter?

-- 
   O__   Peter Dalgaard Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics 2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark  Ph: (+45) 35327918
~~ - ([EMAIL PROTECTED]) FAX: (+45) 35327907

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] r-help-bounces and r-devel-bounces e-mail?

2003-03-04 Thread Marc Schwartz
Hi all,

Just curious on the sudden change in the e-mail address coming from
r-help and r-devel today. They are indicated as being
[EMAIL PROTECTED] and
[EMAIL PROTECTED]

Any clues?  I just noticed that my e-mail filter rule suddenly stopped
putting these messages in their respective folders this morning (my
time).

Regards,

Marc

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] r-help-bounces and r-devel-bounces e-mail?

2003-03-04 Thread Jeff Gentry
 Any clues?  I just noticed that my e-mail filter rule suddenly stopped
 putting these messages in their respective folders this morning (my
 time).

I believe they're upgrading the Mailman software that runs the list.

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] suggestion for addition to R-lang.texi

2003-03-04 Thread A.J. Rossini
Peter Dalgaard BSA [EMAIL PROTECTED] writes:

 BTW, I wonder: why did your mail come with sender r-help-bounces
 instead of the usual r-help-admin and confuse my mail filter?

Mailman upgrade?

best,
-tony

-- 
A.J. RossiniRsrch. Asst. Prof. of Biostatistics
U. of Washington Biostatistics  [EMAIL PROTECTED]   
FHCRC/SCHARP/HIV Vaccine Trials Net [EMAIL PROTECTED]
-- http://software.biostat.washington.edu/ 
FHCRC:Tu: 206-667-7025 (fax=4812)|Voicemail is pretty sketchy/use Email
UW:   Th: 206-543-1044 (fax=3286)|Change last 4 digits of phone to FAX
(CHANGE: monday/wednesday/friday locations are completely unpredictable.)

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


RE: [R] r-help-bounces and r-devel-bounces e-mail?

2003-03-04 Thread Marc Schwartz
-Original Message-
From: Charles Annis, P.E. 
[mailto:[EMAIL PROTECTED] 
Sent: Tuesday, March 04, 2003 11:06 AM
To: 'Jeff Gentry'; 'Marc Schwartz'
Cc: [EMAIL PROTECTED]
Subject: RE: [R] r-help-bounces and r-devel-bounces e-mail?


You need to add [EMAIL PROTECTED] to your 
filtering rules to augment your current rule, 
[EMAIL PROTECTED]

Charles Annis, P.E.


Easy enough to do of course. Just curious if this was a temporary (as
Jeff note suggests) or a permanent change.  The 'bounces' suggests a
problem, which raised my suspicions.

From Peter Dalgaard's reply earlier, it would seem to be unexpected.
I wasn't sure if I had perhaps missed an announcement somewhere.

Thanks,

Marc

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] r-help-bounces and r-devel-bounces e-mail?

2003-03-04 Thread Martin Maechler
 JeffG == Jeff Gentry [EMAIL PROTECTED]
 on Tue, 4 Mar 2003 11:53:13 -0500 (EST) writes:

 Any clues?  I just noticed that my e-mail filter rule suddenly stopped
 putting these messages in their respective folders this morning (my
 time).

JeffG I believe they're upgrading the Mailman software that runs the list.

Exactly.  They := { (our head sys.admin.), me }

The new mailman (2.1.1 instead of 2.0.13) has quite a few new
features; particularly it says it allows unsubscribing without a
password --- something which will hopefully save quite a bit of
my time {from answering e-mails of people who are not patient
enough to read a page on the web to its end}.


I'm sorry we overlooked the change from 
Sender: list[EMAIL PROTECTED]
to
Sender: list[EMAIL PROTECTED]

BTW: I think there are still some new problems that we hadn't
 seen on the testing machine on which we tried everything.
But now is past the end of the working day around here...  
Please send more feedback directly to me, aka  [EMAIL PROTECTED]

Regards,
Martin Maechler [EMAIL PROTECTED] http://stat.ethz.ch/~maechler/
Seminar fuer Statistik, ETH-Zentrum  LEO C16Leonhardstr. 27
ETH (Federal Inst. Technology)  8092 Zurich SWITZERLAND
phone: x-41-1-632-3408  fax: ...-1228   

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] implementing ecological models in R

2003-03-04 Thread Thomas Petzoldt
Hi all,

we found, that R is not only *the* superior data analysis, graphics and 
statistics package, but it is also a general tool for implementing, 
running and teaching ecological models. Last week we held our second 
course Modelling in Limnology using R. Some material (code and docs) 
was produced for this purpose and may be interesting to someone else, so 
I put the tutorial Konstruktion ökologischer Modelle mit der Open 
Source Software R on my homepage.

The tutorial covers
* a short introduction,
* some simple regression models,
* the modelling of populations with ODEs (thanks to Woodrow Sezter!),
* some basic lake modelling stuff and
* some basic individual-based techniques (population dynamics, movement, 
Conway's game of life, an individual-based predator-prey model).

All examples are given in full source code. The tutorial is in German 
language, but possibly someone may be interested to work together with 
me on an English version in the future. It is work in progress, so 
please excuse typos and mistakes (some PDF errors are already known).
On the other hand, I would be very glad about remarks and constructive 
criticism.

The tutorial can be found on

http://www.tu-dresden.de/fghhihb/petzoldt/zeugs.html

together with an absolutely preliminary version of an accompanying 
library (the name simecol may be changed soon).

Thank you for your feedback!

Many thanks to all the people who are making the R-project possible!

Thomas

--
Thomas PetzoldtTel. +49-351-463 3 4954
Dresden University of Technology   Fax  +49-351-463 3 7108
Institute for Hydrobiology [EMAIL PROTECTED]
01062 Dresden
http://www.tu-dresden.de/fghhihb/petzoldt/

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


RE: [R] r-help-bounces and r-devel-bounces e-mail?

2003-03-04 Thread Marc Schwartz
-Original Message-
From: [EMAIL PROTECTED] 
[mailto:[EMAIL PROTECTED] On Behalf Of Martin
Maechler
Sent: Tuesday, March 04, 2003 11:42 AM
To: [EMAIL PROTECTED]
Subject: Re: [R] r-help-bounces and r-devel-bounces e-mail?


 JeffG == Jeff Gentry [EMAIL PROTECTED]
 on Tue, 4 Mar 2003 11:53:13 -0500 (EST) writes:

 Any clues?  I just noticed that my e-mail filter rule 
suddenly stopped
 putting these messages in their respective folders this 
morning (my
 time).

JeffG I believe they're upgrading the Mailman software 
that runs the list.

Exactly.  They := { (our head sys.admin.), me }

The new mailman (2.1.1 instead of 2.0.13) has quite a few new 
features; particularly it says it allows unsubscribing without 
a password --- something which will hopefully save quite a bit 
of my time {from answering e-mails of people who are not 
patient enough to read a page on the web to its end}.


I'm sorry we overlooked the change from 
Sender: list[EMAIL PROTECTED]
to
Sender: list[EMAIL PROTECTED]

BTW: I think there are still some new problems that we hadn't
 seen on the testing machine on which we tried everything. 
But now is past the end of the working day around here...  
Please send more feedback directly to me, aka  [EMAIL PROTECTED]

Regards,
Martin Maechler [EMAIL PROTECTED]


Martin,

Thanks for the clarifications. Most importantly thanks for your energy
in maintaining and continuously improving the lists!

Best regards,

Marc

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] tseries contains a class for irregularly spaced time series

2003-03-04 Thread Adrian Trapletti
A new version of tseries (0.9-10) has been uploaded to CRAN. The new
version contains the class irts for irregularly spaced time series.
Irregular time series are basically time series where each observation
(uni- or multivariate) has a time-stamp represented by an object of
class POSIXct. It provides some basic functionality such as reading
and writing irregular time series from files, or plotting, printing,
subscripting, and interpolating irregular time series. This is a first
version of the class irts and I very much welcome feedback.

best
Adrian

--
Dr. Adrian Trapletti Phone : +41 (0) 1 994 5631
Trapletti Statistical Computing  Mobile: +41 (0)76 370 5631
Wildsbergstrasse 31  Fax   : +41 (0) 1 994 5631
CH-8610 UsterEmail :  mailto:[EMAIL PROTECTED]
Switzerland  WWW   : http://trapletti.homelinux.com

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] Sample size and stepAIC, step, or AIC

2003-03-04 Thread Martin Henry H. Stevens
Do any R functions incorporate a sample sample size correction (e.g., 
Burnham and Anderson 1998).
Thanks,
Hank Stevens

Martin Henry H. Stevens, Assistant Professor
338 Pearson Hall
Botany Department
Miami University
Oxford, OH 45056
Office: (513) 529-4206
Lab: (513) 529-4262
FAX: (513) 529-4243
http://www.cas.muohio.edu/botany/bot/henry.html
http://www.muohio.edu/ecology
__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] problem at installing Rade

2003-03-04 Thread Agustin Lobo

Hi!

I try to install the R version of ade4 on a linux (suse7.3) box
with R 1.6.2 (from rpm)
and get 
/usr/local/bin/R INSTALL ade4_1.00.tar.gz
tar: /usr/local/lib/libpthread.so.0: version `GLIBC_2.2' not found 
(required by /lib/librt.so.1)
ERROR: cannot extract package from 'ade4_1.00.tar.gz'

But I actually have glibc 2.2.4 
(at least this is what my suse's Yast conf. tool 
tells me, perhaps there is a more direct way of checking 
it?)

Having used the rpm for R, should I use the
rpm for ade4? 

Thanks 

Dr. Agustin Lobo
Instituto de Ciencias de la Tierra (CSIC)
Lluis Sole Sabaris s/n
08028 Barcelona SPAIN
tel 34 93409 5410
fax 34 93411 0012
[EMAIL PROTECTED]

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] Parallel programming using R

2003-03-04 Thread lun li
Thank you very much for  the information of  RPvm and RMPI packages. I got 
them.

Best regards,

Lun






From: M.Kondrin [EMAIL PROTECTED]
To: lun li [EMAIL PROTECTED]
CC: R-Help [EMAIL PROTECTED]
Subject: Re: [R] Parallel programming using R
Date: Tue, 04 Mar 2003 14:25:02 -0800
lun li wrote:
Dear All,

For my regression problem(thousand dependents by thousand predictors by 
hundred months observations), I prefer to  parallel computation in R.  Is 
it possible to parallel programming using R as using C and Fortran?

Cheers,

Lun

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help
Look up  RPvm and RMPI packages

Good luck ...

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] Disregard:problem at installing Rade

2003-03-04 Thread Agustin Lobo

Excuses to the list,

this problem is a peculiarity of my system that
I (thought) had solved, but after installing 
a (non-R) package seems to be here again.

Agus

Dr. Agustin Lobo
Instituto de Ciencias de la Tierra (CSIC)
Lluis Sole Sabaris s/n
08028 Barcelona SPAIN
tel 34 93409 5410
fax 34 93411 0012
[EMAIL PROTECTED]

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] Sample size and stepAIC, step, or AIC

2003-03-04 Thread ripley
They use AIC as defined by its author, Akaike.  With good reason: `sample
size' is not a well-defined concept in general, and the R functions are
general.

Most people call the adjusted (I hesitate to call it `corrected') version 
`AICC'.  I believe the `clarified' 2002 edition of Burnham  Anderson uses 
AIC_c.


(It's easy to set up problems, e.g, a binomial glm, in which the number of 
observations relevant to each parameter differs radically.  They occur in 
real life too.)

On Tue, 4 Mar 2003, Martin Henry H. Stevens wrote:

 Do any R functions incorporate a sample sample size correction (e.g., 
 Burnham and Anderson 1998).

-- 
Brian D. Ripley,  [EMAIL PROTECTED]
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel:  +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UKFax:  +44 1865 272595

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] (no subject)

2003-03-04 Thread ripley
On Tue, 4 Mar 2003, Usman Shehu wrote:

 It is possible I run R on MINGW32 MSYS? If so how do I configure the MSYS to do 
 that. 

MSYS is an environment for building Unix packages under Windows.  R runs
under Windows:  you don't need Unix-emulation tools to run it.

We have a much more reliable set of Unix-emulation tools to build 
it, selected from Cygwin (of which MSYS is a derivative, anyway).

-- 
Brian D. Ripley,  [EMAIL PROTECTED]
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel:  +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UKFax:  +44 1865 272595

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] How to extract R{i} from lme object?

2003-03-04 Thread Peng
Hi, lme() users,

Can some one tell me how to do this.
I model Orthodont with the same G for random
variables, but different R{i}'s for boys and girls, so
that I can get sigma1_square_hat for boys and
sigma2_square_hat for girls.

The model is Y{i}=X{i}beta + Z{i}b + e{i}
b ~ iid N(0,G) and e{i} ~ iid N(0,R{i}) i=1,2
orth.lme - lme(distance ~ Sex * age, data=Orthodont,
random=~age|Subject, weights=varIdent(form=~1|Sex),
method=ML)

I can see the numbers I need from summary(), but how
can I extract them? I tried several functions in nlme,
but I cannot find a correct one.

Peng

Peng Liu
--
Peng Liu  |
Division of Statistics|
Northern Illinois University  |
De Kalb, IL 60115, USA|
E-mail: [EMAIL PROTECTED] |
--

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] How to extract R{i} from lme object?

2003-03-04 Thread Sundar Dorai-Raj


Peng wrote:
Hi, lme() users,

Can some one tell me how to do this.
I model Orthodont with the same G for random
variables, but different R{i}'s for boys and girls, so
that I can get sigma1_square_hat for boys and
sigma2_square_hat for girls.
The model is Y{i}=X{i}beta + Z{i}b + e{i}
b ~ iid N(0,G) and e{i} ~ iid N(0,R{i}) i=1,2
orth.lme - lme(distance ~ Sex * age, data=Orthodont,
random=~age|Subject, weights=varIdent(form=~1|Sex),
method=ML)
I can see the numbers I need from summary(), but how
can I extract them? I tried several functions in nlme,
but I cannot find a correct one.
Peng

Peng,

Is this what you need?

R data(Orthodont)
R orth.lme - lme(distance ~ Sex * age,
+  data=Orthodont,
+  random=~age|Subject,
+  weights=varIdent(form=~1|Sex),
+  method=ML)
R orth.lme$modelStruct$varStruct
Variance function structure of class varIdent representing
 MaleFemale
1.000 0.4112708
R rcov.unscaled - as.matrix(orth.lme$modelStruct$reStruct[[1]])
R rcov.scaled - rcov.unscaled * orth.lme$sigma^2
R Stdev - sqrt(diag(rcov.scaled))
R Stdev
(Intercept) age
  1.7914848   0.1408001
R
If there's nesting, look for more elements in $reStruct.

Regards,
Sundar
__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] How to extract R{i} from lme object?

2003-03-04 Thread Peng

--- Sundar Dorai-Raj [EMAIL PROTECTED] wrote:
 
 
 Peng wrote:
  Hi, lme() users,
  
  Can some one tell me how to do this.
  I model Orthodont with the same G for random
  variables, but different R{i}'s for boys and
 girls, so
  that I can get sigma1_square_hat for boys and
  sigma2_square_hat for girls.
  
  The model is Y{i}=X{i}beta + Z{i}b + e{i}
  b ~ iid N(0,G) and e{i} ~ iid N(0,R{i}) i=1,2
  orth.lme - lme(distance ~ Sex * age,
 data=Orthodont,
  random=~age|Subject,
 weights=varIdent(form=~1|Sex),
  method=ML)
  
  I can see the numbers I need from summary(), but
 how
  can I extract them? I tried several functions in
 nlme,
  but I cannot find a correct one.
  
  Peng
  
 
 Peng,
 
 Is this what you need?
 
 R data(Orthodont)
 R orth.lme - lme(distance ~ Sex * age,
 +  data=Orthodont,
 +  random=~age|Subject,
 +  weights=varIdent(form=~1|Sex),
 +  method=ML)
 R orth.lme$modelStruct$varStruct
 Variance function structure of class varIdent
 representing
   MaleFemale
 1.000 0.4112708

Sundar, 

Thanks a lot!
In the above model, I assume the error term for boy
has iid N(0, sigma1), and that for girl has iid
N(0,sigma2). I hope that I wrote the correct lme(). If
so, what I want is:
(1/unique(attributes(orth.lme$modelStruct$varStruct)$weights)*orth.lme$sigma)^2
[1] 2.6292575 0.4447223
, where 2.629 for boys, and 0.445 for girls.

I am still wondering whether there is a function like
getVarCov(), so that I can get R var-cov matrix
directly. I looked through the function list of nlme,
but I cannot find.

Peng

 R rcov.unscaled -
 as.matrix(orth.lme$modelStruct$reStruct[[1]])
 R rcov.scaled - rcov.unscaled * orth.lme$sigma^2
 R Stdev - sqrt(diag(rcov.scaled))
 R Stdev
 (Intercept) age
1.7914848   0.1408001
 R
 
 If there's nesting, look for more elements in
 $reStruct.
 
 Regards,
 Sundar
 
 __
 [EMAIL PROTECTED] mailing list
 http://www.stat.math.ethz.ch/mailman/listinfo/r-help

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] Parallel programming using R

2003-03-04 Thread Michael Na Li
On Tue, 04 Mar 2003, lun li told this:
  Thank you very much for  the information of  RPvm and RMPI packages. I got
  them.
  

You might want to wait for a few days for snow (Simple Network of
Workstations) package to reach CPAN. Or download it from here,

http://www.stat.uiowa.edu/~luke/R/cluster/cluster.html

Cheers,

Michael

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] logistic regression for repeated measurement

2003-03-04 Thread Hiroto Miyoshi
Dear R-users
(B
(BI need your help.
(BI have a data set which was collected from 
(Ban experiment of one between- and one 
(Bwithin-subject design. And the response 
(Bdata is coded by success(1)/failure(0).
(B
(BThe experiment had two groups of subjects:
(BThe one was experimental, and the 
(Bother, control.  The experimental group
(Bgot a task training, and both groups of subjects
(Bwere tested twice, once before the training
(Band once after the training. at the same time.
(BI like to examine the effect of training by
(Bdetecting an interaction effect of the group and 
(Btests. 
(BNow, it seems glm is not appropriate to this 
(Bsituation since it does not deal with stratified
(Berrors.
(B
(BCould you lead me to appropriate functions?
(BSincerely
(B---
(BHiroto Miyoshi
$B;09%90?M(B
(B[EMAIL PROTECTED]
(B
(B__
(B[EMAIL PROTECTED] mailing list
(Bhttp://www.stat.math.ethz.ch/mailman/listinfo/r-help

Re: [R] logistic regression for repeated measurement

2003-03-04 Thread Frank E Harrell Jr
On Wed, 5 Mar 2003 11:14:37 +0900
Hiroto Miyoshi [EMAIL PROTECTED] wrote:

 Dear R-users
 
 I need your help.
 I have a data set which was collected from 
 an experiment of one between- and one 
 within-subject design. And the response 
 data is coded by success(1)/failure(0).
 
 The experiment had two groups of subjects:
 The one was experimental, and the 
 other, control.  The experimental group
 got a task training, and both groups of subjects
 were tested twice, once before the training
 and once after the training. at the same time.
 I like to examine the effect of training by
 detecting an interaction effect of the group and 
 tests. 
 Now, it seems glm is not appropriate to this 
 situation since it does not deal with stratified
 errors.
 
 Could you lead me to appropriate functions?
 Sincerely
 ---
 Hiroto Miyoshi
 [EMAIL PROTECTED]
 

There are several ways to go.  GEE is one, random effects models another.  One other 
approach is to install the Hmisc and Design packages 
(http://hesweb1.med.virginia.edu/biostat/s) and do (assume id is the unique subject 
identifier):

f - lrm(y ~ x1 + x2*x3 + ..., x=T, y=T) # working independence model
g - robcov(f, id) # cluster sandwich variance adjustment
h - bootcov(f, id, B=100) # cluster bootstrap adjustment
summary(g)  # etc.
-- 
Frank E Harrell Jr  Prof. of Biostatistics  Statistics
Div. of Biostatistics  Epidem. Dept. of Health Evaluation Sciences
U. Virginia School of Medicine  http://hesweb1.med.virginia.edu/biostat

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] How to draw several plots in one figure?

2003-03-04 Thread Feng Zhang
Hey,

I want to draw several plots sequently, but have to make them dispaly in one
figure.
So how to achieve this?

Thanks.

Fred

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] How to draw several plots in one figure?

2003-03-04 Thread Feng Zhang
Thanks, Su.

But I want to plot the several plots in the same
x-y axis setting, not in multiple subplots.
- Original Message -
From: Steve Su [EMAIL PROTECTED]
To: Feng Zhang [EMAIL PROTECTED]
Sent: Wednesday, March 05, 2003 12:29 AM
Subject: Re: [R] How to draw several plots in one figure?


 Dear Fred,

 try par(mfrow=c(2,2)) for example, gives four plots per page. Also try
 trellis plots.

 Steve.


 - Original Message -
 From: Feng Zhang [EMAIL PROTECTED]
 To: R-Help [EMAIL PROTECTED]
 Sent: Wednesday, March 05, 2003 4:05 PM
 Subject: [R] How to draw several plots in one figure?


  Hey,
 
  I want to draw several plots sequently, but have to make them dispaly in
 one
  figure.
  So how to achieve this?
 
  Thanks.
 
  Fred
 
  __
  [EMAIL PROTECTED] mailing list
  http://www.stat.math.ethz.ch/mailman/listinfo/r-help


__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] How to draw several plots in one figure?

2003-03-04 Thread Jeremy Z Butler
One way is to make inital plot:
plot()
Then build other data into the plot using:

lines()
points()
etc.
see help files for those
Alternatively:

set par(new=T)
then plot again and the device wont be cleared
J.



From: Feng Zhang [EMAIL PROTECTED]
To: Steve Su [EMAIL PROTECTED]
CC: R-Help [EMAIL PROTECTED]
Subject: Re: [R] How to draw several plots in one figure?
Date: Wed, 5 Mar 2003 00:35:00 -0600
Thanks, Su.

But I want to plot the several plots in the same
x-y axis setting, not in multiple subplots.
- Original Message -
From: Steve Su [EMAIL PROTECTED]
To: Feng Zhang [EMAIL PROTECTED]
Sent: Wednesday, March 05, 2003 12:29 AM
Subject: Re: [R] How to draw several plots in one figure?
 Dear Fred,

 try par(mfrow=c(2,2)) for example, gives four plots per page. Also try
 trellis plots.

 Steve.


 - Original Message -
 From: Feng Zhang [EMAIL PROTECTED]
 To: R-Help [EMAIL PROTECTED]
 Sent: Wednesday, March 05, 2003 4:05 PM
 Subject: [R] How to draw several plots in one figure?


  Hey,
 
  I want to draw several plots sequently, but have to make them dispaly 
in
 one
  figure.
  So how to achieve this?
 
  Thanks.
 
  Fred
 
  __
  [EMAIL PROTECTED] mailing list
  http://www.stat.math.ethz.ch/mailman/listinfo/r-help


__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help
__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help


BSOD with ESS[R] under win2000

2003-03-04 Thread Darryl Greig
I have had a visit from our windows buddy the BSOD twice in one day when
running
R1.6.2 under ESS5.1.24 on a win2000 platform. I wasn't doing anything
special, just
moving or resizing the emacs window, although I did have 4 or 5 emacs
windows
open at the time. Has anyone else observed this problem?

Thanks, Darryl.

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help