[R] Manova and specifying the model

2004-05-24 Thread Stefanie von Felten
Hi,

I would like to conduct a MANOVA. I know that there 's the manova() funciton and the 
summary.manova() function to get the appropriate summary of test statistics.

I just don't manage to specify my model in the manova() call. How to specify a model 
with multiple responses and one explanatory factor?

If I type:
pcor.manova-manova(isol+hcom+habarea+inclin+windprot+shrubcov+herbh+baregr+flowcov~pcor,
 data=pcor.df)

I always get error messages like:
Error in manova(isol + hcom + habarea + inclin + windprot + shrubcov +  : 
need multiple response

Can someone help me?

Thanks
Steffi




Stefanie von Felten
Käppelistrasse 24
4600 Olten
Tel: 062/296 13 14
e-mail: [EMAIL PROTECTED]
[[alternative HTML version deleted]]

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] discriminant analysis

2004-05-24 Thread Stefanie von Felten
Hi,

I have done different discriminant function analysis of multivariat data. With the 
CV=True option I was not able to perform the predict() call. What do I have to do? Or 
is there no possibility at all? You also need the predicted values to produce a plot 
of the analysis, as far as I know.

Here my code:

pcor.lda2-lda(pcor~habarea+hcom+isol+flowcov+herbh+inclin+windprot+shrubcov+baregr, 
data=pcor.df, CV=T)
table2-table(pcor.df$pcor, pcor.lda2$class)
table2

#doesn't work, becoause CV=True?
pcor.ld2-predict(pcor.lda2, dimen=1)$x
plot(pcor.ld2)
plot(pcor.lda2, type=density, dimen=1) #kernel density estimates 

I am happy if I get an answer from somebody!


Stefanie von Felten
Institut für Umweltwissenschaften
Universität Zürich
Winterthurerstrasse 190
CH-8057 Zürich

e-mail: [EMAIL PROTECTED]
telefon: +41 (0)1 635 61 23
felefax: +41 (0)1 635 57 11
[[alternative HTML version deleted]]

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] Manova and specifying the model

2004-05-24 Thread Prof Brian Ripley
See the examples in ?summary.manova, please.

Y - cbind(tear, gloss, opacity)

shows how to make a multiple response.

On Mon, 24 May 2004, Stefanie von Felten wrote:

 I would like to conduct a MANOVA. I know that there 's the manova()
 funciton and the summary.manova() function to get the appropriate
 summary of test statistics.
 
 I just don't manage to specify my model in the manova() call. How to
 specify a model with multiple responses and one explanatory factor?
 
 If I type:
 pcor.manova-manova(isol+hcom+habarea+inclin+windprot+shrubcov+herbh+baregr+flowcov~pcor,
  data=pcor.df)
 
 I always get error messages like:
 Error in manova(isol + hcom + habarea + inclin + windprot + shrubcov +  : 
 need multiple response


-- 
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
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] discriminant analysis

2004-05-24 Thread Prof Brian Ripley
On Mon, 24 May 2004, Stefanie von Felten wrote:

 I have done different discriminant function analysis of multivariat

Using lda in contributed package MASS, uncredited.

 data. With the CV=True option I was not able to perform the predict()
 call. What do I have to do? Or is there no possibility at all? You also

It makes no sense.  You ask for LOO cross-validation, and that is n 
separate fits, not a single fit from which you can predict.

 need the predicted values to produce a plot of the analysis, as far as I
 know.
 
 Here my code:
 
 pcor.lda2-lda(pcor~habarea+hcom+isol+flowcov+herbh+inclin+windprot+shrubcov+baregr, 
 data=pcor.df, CV=T)
 table2-table(pcor.df$pcor, pcor.lda2$class)
 table2
 
 #doesn't work, becoause CV=True?
 pcor.ld2-predict(pcor.lda2, dimen=1)$x
 plot(pcor.ld2)
 plot(pcor.lda2, type=density, dimen=1) #kernel density estimates 
 
 I am happy if I get an answer from somebody!

-- 
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
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] Cannot call R's ISNAN() from a C code in 1.7 versions.

2004-05-24 Thread Prof Brian Ripley
First, there are no versions 1.7, 1.8 and 1.9.

Was your version of R compiled against MSVC++ 6.0?  The binary on CRAN was
not, and binaries for different versions of R were compiled with different
versions of MinGW.  The entry point isnan is part of the statically linked
runtime on modern MinGW.

MSVC++ 6.0 does supply _isnan (as it really should as it is part of the
C99 ISO standard), and you need to link against it appropriately.  Hint:
it may have an extra underline, since it seems it is known to C as _isnan.
You may need to add

#undef ISNAN
#define ISNAN(x) _isnan(x)


We don't support adding extensions to R using a different compiler to the 
one used to build R.  Changes already made for future releases of R will 
make this less likely to work in R 2.0.x.


On Mon, 24 May 2004, Paul Y. Peng wrote:

 Dear R users,
 
 Have you experienced any difficulty in calling R's ISNAN() from a C
 code? I have C codes including ISNAN() calls and they worked well until
 I upgraded my R from 1.7 to later versions. When I tried to compile the
 codes in the version 1.8 and 1.9, I got error messages like this:
 
test.obj : error LNK2001: unresolved external symbol _isnan
.\testR.dll : fatal error LNK1120: 1 unresolved externals
NMAKE : fatal error U1077: 'link.exe' : return code '0x460'
Stop.
 
 I checked Writing R Extensions and did not find any changes in the API
 entry point for ISNAN in the later versions. Could any one enlighten me
 on why _isnan cannot be resolved? The compiler I used is MSVC++ 6.0 and
 the platform is WinXP. Many thanks.

-- 
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
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] bug in cor (..., use= ...)?

2004-05-24 Thread Lorenz . Gygax
Dear R users,

I have not found anything on this in the archives. Does anyone know whehther
the parameter use= is not functioning in cor or enlighten me what it is
supposed to do?

My R version is R version 1.8.1, 2003-11-21 on Windows 2000. I am hoping
to be able to update to 1.9.1 as soon as it has appeared (we are not allowed
here to install software on our own and thus I am trying to be able to have
the .1 versions installed ...).

Test code:

x - 1:10
y - 2:11

x [1] - NA
y [10] - 12

cor (x, y, use= 'all.obs', method= 'kendall')
cor (x, y, use= 'complete.obs', method= 'kendall')
cor (x, y, use= 'pairwise.complete.obs', method= 'kendall')

As I understand, the first one of this should result in an error which it
does not. All the results are the same and seemingly treat the NA as if it
was 0.

Any ideas are appreciated.

Thanks and regards, Lorenz
- 
Lorenz Gygax, Dr. sc. nat.
Tel: +41 (0)52 368 33 84 / [EMAIL PROTECTED]  

Tag der offenen Tür, 11./12. Juni 2004: http://www.fat.ch/2004

Center for proper housing of ruminants and pigs
Swiss Veterinary Office
agroscope FAT Tänikon, CH-8356 Ettenhausen / Switzerland
Fax : +41 (0)52 365 11 90 / Tel: +41 (0)52 368 31 31

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] bug in cor (..., use= ...)?

2004-05-24 Thread Prof Brian Ripley
From the NEWS file for 1.9.0

o   The cor() function did not remove missing values in the
non-Pearson case.

Your example works correctly there.  (I am fairly sure this has been 
discussed on the mailing lists.)

On Mon, 24 May 2004 [EMAIL PROTECTED] wrote:

 Dear R users,
 
 I have not found anything on this in the archives. Does anyone know whehther
 the parameter use= is not functioning in cor or enlighten me what it is
 supposed to do?
 
 My R version is R version 1.8.1, 2003-11-21 on Windows 2000. I am hoping
 to be able to update to 1.9.1 as soon as it has appeared (we are not allowed
 here to install software on our own and thus I am trying to be able to have
 the .1 versions installed ...).
 
 Test code:
 
 x - 1:10
 y - 2:11
 
 x [1] - NA
 y [10] - 12
 
 cor (x, y, use= 'all.obs', method= 'kendall')
 cor (x, y, use= 'complete.obs', method= 'kendall')
 cor (x, y, use= 'pairwise.complete.obs', method= 'kendall')
 
 As I understand, the first one of this should result in an error which it
 does not. All the results are the same and seemingly treat the NA as if it
 was 0.

-- 
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
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] Tramo-seats

2004-05-24 Thread v . demartino2
Working  - among other things- in the field of (short  long term) electricity
forecast,  we are now using too many  too expensive pieces of licensed
software: SAS, SPSS, EViews. This sedimentation is due to the fact that
my predecessors in the past used different consultant companies to manage
each procedure.

Having attended the useR2004! Conference with the aim of assessing if R
ALONE could glue all those fragmented and isolated procedures, I'm almost
convinced now that YES it could do the job!

Now - to start with - a first problem to solve:

we have to comply with the Tramo-seats closed-source procedure 
(http://www.bde.es/informes/be/docs/dt0014e.pdf)
to deal with seasonality of electricity monthly time-series, in line with
the methodology officially adopted by our  National Bureau of Statistics.

Searching in R-help mailing list I didn't find anything about a tramo-seats
R version. 

Does anyone know of unofficial R translation of tramo-seats?
OR What do you suggest?

Vittorio

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] problems with starting R

2004-05-24 Thread Sabine Bader
When I try to start R from my desktop, an information-window pops up:
Fatal error: Invalid HOMEDRIVE.
I already deleted and reinstalled the program, but the problem lasts.
How can I solve this problem?
I would be very glad, if you can give me a hint, what´s going wrong.

Best regards
Sabine Bader

-- 
Sabine Bader
Hamburger Str. 83
44135 Dortmund

0231/7950683

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


RE: [R] problems with starting R

2004-05-24 Thread Rau, Roland
Hi Sabine,

 -Original Message-
 From: Sabine Bader [SMTP:[EMAIL PROTECTED]
 Sent: Monday, May 24, 2004 2:22 PM
 To:   [EMAIL PROTECTED]
 Subject:  [R] problems with starting R
 
 When I try to start R from my desktop, an information-window pops up:
 Fatal error: Invalid HOMEDRIVE.
 
This problem appeared dozens of times during the last few weeks.
Best thing is probably to read one of the previous threads related
to this topic.
If you go, for example, to: 
http://maths.newcastle.edu.au/~rking/R/
there is a searchable R-help archive.
Enter the keywords:
INVALID HOMEDRIVE
and you will get many many results (I just tried it and it gave 128
hits).

Hope this helps,
Roland



+
This mail has been sent through the MPI for Demographic Rese...{{dropped}}

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] problems with starting R

2004-05-24 Thread Wolski
Hi!
This problem was discussed on the list.
For searching the mail archives you can use eg.

http://maths.newcastle.edu.au/~rking/R/

The link to this page can be found at:
http://cran.r-project.org/search.html

Entering 
Fatal error: Invalid HOMEDRIVE.
Will provide you many answers to your question.

Eg.:
http://tolstoy.newcastle.edu.au/R/help/04/04/1171.html

Sincerely
Eryk




*** REPLY SEPARATOR  ***

On 5/24/2004 at 2:22 PM Sabine Bader wrote:

When I try to start R from my desktop, an information-window pops up:
Fatal error: Invalid HOMEDRIVE.
I already deleted and reinstalled the program, but the problem lasts.
How can I solve this problem?
I would be very glad, if you can give me a hint, what´s going wrong.

Best regards
Sabine Bader

-- 
Sabine Bader
Hamburger Str. 83
44135 Dortmund

0231/7950683

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html



Dipl. bio-chem. Eryk Witold Wolski@MPI-Moleculare Genetic   
Ihnestrasse 63-73 14195 Berlin   'v'
tel: 0049-30-83875219   /   \
mail: [EMAIL PROTECTED]---W-Whttp://www.molgen.mpg.de/~wolski

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] non-hierarchical non-exclusive clustering of large data sets

2004-05-24 Thread Murk Wuite
Hi,

I'm trying to use R to cluster words with related meanings. Does anyone
know of a non-hierarchical clustering method in R that produces
non-exclusive clusters? With non-exclusive, I mean that words should be
allowed to be part of multiple clusters. So my data matrix would look
something like:

T1  T2  T3
CLOWN_N 0   1   0
BANK_N  3   0   2
RIVER_N 0   0   2
FLOW_V  0   0   3
MONEY_N 2   0   0
PAY_V   2   0   0

The first line indicates the noun clown occurred only once in my text
collection, namely in text 2. Ideally, the clustering method would
produce the clusters [bank_n,river_n,flow_v], [bank_n,money_n,pay_v] and
[clown_n].
The data matrix I would use would be much bigger than the one above, its
dimensions would be in the order of (10,10). Does anyone know if
this would cause practical problems, perhaps very slow clustering?

Best wishes,

Murk Wuite, MA student
Department of Language and Speech
Katholieke Universiteit Nijmegen, The Netherlands

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] replacing backslashes with slashes using gsub

2004-05-24 Thread witek
Hi!

I am trying to replace backslashes with slashes using gsub (R1.9.0 on XP)

gsub(,/,D:\Prog\R\rw1090\library\cluster\libs)
[1] D:ProgR\rw1090libraryclusterlibs

?

Sincerely Eryk

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] Month names

2004-05-24 Thread BXC (Bendix Carstensen)
This is how I get the month names from within R:

 mon - rep(strptime(01/01/1952, format = %d/%m/%Y), 12)
 mon$mon - mon$mon + 0:11
 mnam - months(mon, abbreviate = F)
 mnam
 [1] januarfebruar   marts april maj   juni
juli  augustseptember
[10] oktober   november  december 

Surely someone on the list can beat that in elegance?

Bendix Carstensen
--
Bendix Carstensen
Senior Statistician
Steno Diabetes Center
Niels Steensens Vej 2
DK-2820 Gentofte
Denmark
tel: +45 44 43 87 38
mob: +45 30 75 87 38
fax: +45 44 43 07 06
[EMAIL PROTECTED]
www.biostat.ku.dk/~bxc

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] replacing backslashes with slashes using gsub

2004-05-24 Thread Jeff Gentry
 gsub(,/,D:\Prog\R\rw1090\library\cluster\libs)
 [1] D:ProgR\rw1090libraryclusterlibs

Probably not the best way, but what about escaping all the backslashes in
the original string?
 gsub(,/,D:\\Prog\\R\\rw1090\\library\\cluster\\libs)

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] (no subject)

2004-05-24 Thread Jeff Gentry
 Hello! Please!How do I download R from the internet?

http://cran.r-project.org/doc/FAQ/R-FAQ.html#How%20can%20R%20be%20obtained%3f

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] (no subject)

2004-05-24 Thread Uwe Ligges
Whittaker, Carly wrote:
Hello! Please!How do I download R from the internet?
Several ways, one ist to right click and say something like save to 
... in your browser when visiting  CRAN at http://cran.r-project.org/ 
and having browsed to the sources or the binary version that is the 
right one for your OS.

Uwe Ligges

Thankyou for your time
carly
[[alternative HTML version deleted]]
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] Stopping the process after a certain time

2004-05-24 Thread Andrew Robinson
Greetings R-community,

I'm running simulations within R that I wrote in C.  The simulations
require fitting that occasionally fails to finish.  I was wondering if
there is any kind of tool for process control in R, such that after e.g.
15 minutes I could kill the process, record the state for post-hoc
analysis, and move to the next simulation?

I'm running FreeBSD and could almost surely do something in Perl but I'd
rather stay inside R if possible.

Thanks,

Andrew

Andrew Robinson  Ph: 208 885 7115
Department of Forest Resources   Fa: 208 885 6226
University of Idaho  E : [EMAIL PROTECTED]
PO Box 441133W : http://www.uidaho.edu/~andrewr
Moscow ID 83843  Or: http://www.biometrics.uidaho.edu
No statement above necessarily represents my employer's opinion.

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] Month names

2004-05-24 Thread Roger D. Peng
How about `month.name'?
-roger
BXC (Bendix Carstensen) wrote:
This is how I get the month names from within R:

mon - rep(strptime(01/01/1952, format = %d/%m/%Y), 12)
mon$mon - mon$mon + 0:11
mnam - months(mon, abbreviate = F)
mnam
 [1] januarfebruar   marts april maj   juni
juli  augustseptember
[10] oktober   november  december 

Surely someone on the list can beat that in elegance?
Bendix Carstensen
--
Bendix Carstensen
Senior Statistician
Steno Diabetes Center
Niels Steensens Vej 2
DK-2820 Gentofte
Denmark
tel: +45 44 43 87 38
mob: +45 30 75 87 38
fax: +45 44 43 07 06
[EMAIL PROTECTED]
www.biostat.ku.dk/~bxc
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] Metafiiles into Word R 1.9.0

2004-05-24 Thread Jesus Frias
Dear R-helpers,

I recently upgraded to R 1.9.0 in my computer at work and at home:

1.-The computer at home has Windows XP and Office XP and it seems to work
perfectly and I copy-paste graphics perfectly.

2.-The computer at work has Microsoft Windows 2000 (5.00.2195 Service Pack
2) and Word 2000 (9.0.4402 SR-1) I cannot copy-paste windows metafiles into
Office applications in my computer at work. The resulting object is empty.

Is there any way of overcoming this other than going back in R versions?
(an indication to patch from Microsoft Office that can take new metafiles
:-).

best regards,

Jesus

P.S.: this mail is probably a similar topic as Patrick Giraudoux email at
the beginning of the month, I apologise for the duplication.


platform i386-pc-mingw32
arch i386
os   mingw32
system   i386, mingw32
status
major1
minor9.0
year 2004
month04
day  12
language R




--
Jesús María Frías Celayeta
School of Food Sci. and Env. Health.
Faculty of Tourism and Food
Dublin Institute of Technology
Cathal Brugha St., Dublin 1. Ireland
t +353 1 4024459 f +353 1 4024495
w www.dit.ie/DIT/tourismfood/science/staff/frias.html
--


-- 
This message has been scanned for content and 
viruses by the DIT Information Services MailScanner 
Service, and is believed to be clean.
http://www.dit.ie

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] Month names

2004-05-24 Thread Martin Maechler
 Roger == Roger D Peng [EMAIL PROTECTED]
 on Mon, 24 May 2004 10:38:09 -0400 writes:

Roger How about `month.name'?

English only.

Note that he got the names in Danish (I think)
Martin

Roger BXC (Bendix Carstensen) wrote:
 This is how I get the month names from within R:
 
 
 mon - rep(strptime(01/01/1952, format = %d/%m/%Y), 12)
 mon$mon - mon$mon + 0:11
 mnam - months(mon, abbreviate = F)
 mnam
 
 [1] januarfebruar   marts april maj   juni
 juli  augustseptember
 [10] oktober   november  december 
 
 Surely someone on the list can beat that in elegance?
 
 Bendix Carstensen
 --
 Bendix Carstensen
 Senior Statistician
 Steno Diabetes Center
 Niels Steensens Vej 2
 DK-2820 Gentofte
 Denmark
 tel: +45 44 43 87 38
 mob: +45 30 75 87 38
 fax: +45 44 43 07 06
 [EMAIL PROTECTED]
 www.biostat.ku.dk/~bxc

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] RMySQL problem

2004-05-24 Thread Ernesto Jardim
Hi,

I'm using R 1.9.0 with RMySQL 0.5-4 and MySQL 3.23.55 on a suse 8.2 box.

I have a simulation study and (as usual for newbies in simulation, I
guess) I have a lot of data that I want to store in MySQL. I want to
write an R script that reads data from RData files and writes it to a
MySQL database. 

I read some R documents (R Data Import/Export and DSC papers) but I'm
finding differences between the documents and the packages (RMySQL and
DBI). I don't find the methods to write data like dbWriteTable...

On the other hand I've tryied to build a sql statement to insert data
but I'm stucked because of ... who knows, my ignorance probably. I want
to take advantage of MySQl INSERT INTO statement that deals with several
rows at once to insert a complete data.frame into a table. I've tryied
to use paste to build the sql string but It works by column and I
need it by row ...

The sql systax should be something like:

INSERT INTO TABLEA(COL1, COL2, COL3) VALUES
(VAL11, VAL12, VAL13),
(VAL21, VAL22, VAL23),
...
(VALN1, VALN2, VALN3);

and I have a data.frame with 3 columns corresponding to that table
columns.

How can I do this ?

Thanks

EJ

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] RPM post-install scripts to update R

2004-05-24 Thread Peter Dalgaard
chris albert [EMAIL PROTECTED] writes:

 Hi,
 
 Paul Johnson asked:
 
 While I'm on the RPM subject, can I ask an R RPM packaging question? I
 want the R RPM to install so that post install then R starts and runs
  update.packages()
 
 as well as
 
 install.packages(c(Design,Hmisc,lmtest,car,rgl,effects))
 
 well, you get the idea. I want to add in more packages, of course. Can
 I ask what might be the best way to package this?
 
 You can add scripts to the %post macro section.

...etc...

Erm, but would you really want the result of installing an RPM to
depend on the install date and the previous configuration of the
system?

I can appreciate Paul's desire to automate, but when designing an RPM,
you also have to consider that you'd be doing things that the user
may not even know what mean. That sort of thing easily becomes a
support nightmare. 

-- 
   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
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] Re: Windows versus Unix packages in CRAN ...

2004-05-24 Thread Thomas Lumley
On Mon, 24 May 2004, Peter Dalgaard wrote:

 Whereas we cannot reasonably require people to test code on platforms
 that they haven't got, we can and should require them to adhere to
 reasonable standards and test procedures (which, mind you, other
 people have invested a serious amount of time in working out).


Well, yes, for CRAN.  That's presumably one reason why the packages in
question aren't on CRAN.

-thomas

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] Stopping the process after a certain time

2004-05-24 Thread Thomas Lumley
On Mon, 24 May 2004, Andrew Robinson wrote:

 Greetings R-community,

 I'm running simulations within R that I wrote in C.  The simulations
 require fitting that occasionally fails to finish.  I was wondering if
 there is any kind of tool for process control in R, such that after e.g.
 15 minutes I could kill the process, record the state for post-hoc
 analysis, and move to the next simulation?

 I'm running FreeBSD and could almost surely do something in Perl but I'd
 rather stay inside R if possible.


As one component of this, if you send an R process SIGUSR1 it will quit
and save .RData (approximately the equivalent of CTRL-C, q(yes)).


-thomas

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] Month names

2004-05-24 Thread Peter Dalgaard
Martin Maechler [EMAIL PROTECTED] writes:

  Roger == Roger D Peng [EMAIL PROTECTED]
  on Mon, 24 May 2004 10:38:09 -0400 writes:
 
 Roger How about `month.name'?
 
 English only.
 
 Note that he got the names in Danish (I think)
 Martin

Yep. Notice, however, that even with Brian's solution, you do need to
set the locale first:

 format(ISOdate(2004,1:12,1),%B)
 [1] January   February  March April May   June
 [7] July  AugustSeptember October   November  December
 Sys.setlocale(LC_TIME,da_DK) ;format(ISOdate(2004,1:12,1),%B)
[1] da_DK
 [1] januarfebruar   marts april maj   juni
 [7] juli  augustseptember oktober   november  december


-- 
   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
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] as.matrix.data.frame() in R 1.9.0 converts to character when it should (?) convert to numeric

2004-05-24 Thread Prof Brian Ripley
I don't think a POSIXt element *is* numeric (that's a basic atomic
vector), so the new behaviour seems right to me.  The Warning is wrong,
though, and will be fixed.

On Mon, 24 May 2004, Don MacQueen wrote:

 Conversion of a data frame to a matrix using as.matrix() when a 
 column of the data frame is POSIXt and all other columns are numeric 
 has changed in R 1.9.0 from R 1.8.1. The new behavior issues a 
 warning message and converts to a character matrix. In R 1.8.1, such 
 an object was converted to a numeric matrix.
 
 Here is an example.
 
  R 1.9.0 
   foo - data.frame( x=1:3,dt=ISOdatetime(2003,1,1:3,0,0,0))
 
   as.matrix(foo)
x   dt 
 1 1 2003-01-01
 2 2 2003-01-02
 3 3 2003-01-03
 Warning message:
 longer object length
  is not a multiple of shorter object length in: cl == 
 c(Date, POSIXct, POSIXlt)
 
   version
   _  
 platform sparc-sun-solaris2.8
 arch sparc  
 os   solaris2.8 
 system   sparc, solaris2.8  
 status   Patched
 major1  
 minor9.0
 year 2004   
 month04 
 day  30 
 language R  
 
 
 
 ### R 1.8.1 
   foo - data.frame( x=1:3,dt=ISOdatetime(2003,1,1:3,0,0,0))
   foo
x dt
 1 1 2003-01-01
 2 2 2003-01-02
 3 3 2003-01-03
 
   as.matrix(foo)
x dt
 1 1 1041408000
 2 2 1041494400
 3 3 1041580800
 
   version
   _  
 platform sparc-sun-solaris2.8
 arch sparc  
 os   solaris2.8 
 system   sparc, solaris2.8  
 status   Patched
 major1  
 minor8.1
 year 2003   
 month12 
 day  03 
 language R  
 
 
 
 In both versions:
   class(foo$dt)
 [1] POSIXt  POSIXct
 
 
 In R 1.8.1, as.matrix.data.frame() has these lines:
  if (length(levels(xj))  0 || !(is.numeric(xj) || is.complex(xj)) ||
  (!is.null(cl - attr(xj, class))  any(cl == c(POSIXct,
  POSIXlt
 
 
 In R 1.9.0 there is instead
  if (length(levels(xj))  0 || !(is.numeric(xj) || is.complex(xj)) ||
  (!is.null(cl - attr(xj, class))  any(cl == c(Date,
  POSIXct, POSIXlt
 
 And that, I think, explains the warning message.
 
 
  From ?as.matrix() in R 1.9.0:
 
   'as.matrix' is a generic function. The method for data frames will
   convert any non-numeric/complex column into a character vector
   using 'format' and so return a character matrix, except that
   all-logical data frames will be coerced to a logical matrix.
 
 The POSIXt element is numeric, and so should be converted to numeric
   is.numeric(foo$dt)
 [1] TRUE
 
 
 
 I think this might qualify for bug status, either in and of itself or 
 relative to documentation. But I'm not, as the posting guide says, 
 completely and utterly sure. So I'm posting to r-help first...I 
 will send a bug report if an R-core member asks me to.
 
 Thanks
 -Don
 

-- 
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
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] R 1.8.1 - 1.9.0 incompatability: Underscore in syntactically valid names

2004-05-24 Thread Michael A. Miller
 Peter == Peter Dalgaard [EMAIL PROTECTED] writes:

 Gah! I could swear we discussed that particular issue
 leading up to 1.9.x and had plans for a compatibility
 option.

 You might file a bug report at least for the docs, since
 the example is clearly wrong...

Done.  

I tried to write some version dependencies into a sample code,
but I'm stumped by the fact that _ is not allowed before 1.9.0.
For example, suppose I have a data file, example.dat, like this: 

a b x   some_factor
1 1 0.4 orange
2 1 0.3 blue
1 1 0.2 dog
2 1 0.1 orange
1 2 0.4 blue
2 2 0.3 dog
1 2 0.2 orange
2 2 0.1 blue

To read and use this in a version independent way, I've tried this:

df - read.table('example.dat',header=T)
if ( version['minor'] == 9.0 ) {
  plot(x ~ some_factor, data=df)
} else {
  plot(x ~ some.factor, data=df)
}


This fails in R 1.8.1, because some_factor throws a syntax error:

   df - read.table('example.dat',header=T)
   if ( version['minor'] == 9.0 ) {
  + if ( version['minor'] == 9.0 ) {
  +   plot(x ~ some_factor, data=df)
  Error: syntax error
   

Ick.  Is there a known idiom for handling this sort of version
dependency in R?  I'm going to avoid R 1.9.x for now, and I
encourage any authors of contributing packages to do their best
to maintain backwards compatibility for those of us who cannot
make the switch quickly.

Mike

-- 
Michael A. Miller   [EMAIL PROTECTED]
  Imaging Sciences, Department of Radiology, IU School of Medicine

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] Stopping the process after a certain time

2004-05-24 Thread Andrew Robinson
Thomas,

that's very interesting - thanks!  That will be helpful.

Andrew

On Mon, 24 May 2004, Thomas Lumley wrote:

 On Mon, 24 May 2004, Andrew Robinson wrote:

  Greetings R-community,
 
  I'm running simulations within R that I wrote in C.  The simulations
  require fitting that occasionally fails to finish.  I was wondering if
  there is any kind of tool for process control in R, such that after e.g.
  15 minutes I could kill the process, record the state for post-hoc
  analysis, and move to the next simulation?
 
  I'm running FreeBSD and could almost surely do something in Perl but I'd
  rather stay inside R if possible.
 

 As one component of this, if you send an R process SIGUSR1 it will quit
 and save .RData (approximately the equivalent of CTRL-C, q(yes)).


   -thomas

 __
 [EMAIL PROTECTED] mailing list
 https://www.stat.math.ethz.ch/mailman/listinfo/r-help
 PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Andrew Robinson  Ph: 208 885 7115
Department of Forest Resources   Fa: 208 885 6226
University of Idaho  E : [EMAIL PROTECTED]
PO Box 441133W : http://www.uidaho.edu/~andrewr
Moscow ID 83843  Or: http://www.biometrics.uidaho.edu
No statement above necessarily represents my employer's opinion.

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] Month names

2004-05-24 Thread Prof Brian Ripley
As far as I know the original solution uses the current locale, and mine 
is just a much-simplified version of the same underlying call.

On 24 May 2004, Peter Dalgaard wrote:

 Martin Maechler [EMAIL PROTECTED] writes:
 
   Roger == Roger D Peng [EMAIL PROTECTED]
   on Mon, 24 May 2004 10:38:09 -0400 writes:
  
  Roger How about `month.name'?
  
  English only.
  
  Note that he got the names in Danish (I think)
  Martin
 
 Yep. Notice, however, that even with Brian's solution, you do need to
 set the locale first:

Just as you did with the original code, which gives English names on my 
machine.

  format(ISOdate(2004,1:12,1),%B)
  [1] January   February  March April May   June
  [7] July  AugustSeptember October   November  December
  Sys.setlocale(LC_TIME,da_DK) ;format(ISOdate(2004,1:12,1),%B)
 [1] da_DK
  [1] januarfebruar   marts april maj   juni
  [7] juli  augustseptember oktober   november  december

-- 
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
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] Mac OS X jpg

2004-05-24 Thread Prof Brian Ripley
I think the message is very clear.  Try ?BATCH

 Using 'R CMD BATCH' sets the GUI to 'none', so none of 'x11',
 'jpeg' and 'png' are available.

and ?jpeg

 R can be compiled without support for either or both of these
 devices: this will be reported if you attempt to use them on a
 system where they are not supported.  They will not be available
 if R has been started with '--gui=none' (and will give a different
 error message), and they may not be usable unless the X11 display
 is available to the owner of the R process.

The bitmap() device will work under R CMD BATCH.

On Mon, 24 May 2004, Rafael Najmanovich wrote:

   I am using R 1.9.0 (2004-04-12) for Mac OS X.   I am having trouble 
 printing to a jpg file. This only happens when I use R in batch (with 
 the BATCH option) mode. If I launch the R GUI or the command line 
 version and run the script from either, I have no problems. Only when 
 running in BATCH mode I get a file with the same name as my script 
 where I find the following message at the end:
 
  Error in jpeg(filename = file, width = 1500, height = 2000, pointsize 
  = 12,  :
  R_X11 module cannot be loaded
  In addition: Warning message:
  X11 module is not available under this GUI
  Execution halted
 
   When I run the same script with the only difference that it prints to 
 a postscript file instead, all works fine.
 
   Any  suggestions as to what could the problem be?

Failure to consult the documentation, I believe.

-- 
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
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] Cannot call R's ISNAN() from a C code in 1.7 versions.

2004-05-24 Thread Prof Brian Ripley
On Mon, 24 May 2004, Paul Y. Peng wrote:

 Prof Brian Ripley wrote:
  First, there are no versions 1.7, 1.8 and 1.9.
 
 Sorry for my misuse of the version numbers.
 
  Was your version of R compiled against MSVC++ 6.0?  The binary on CRAN was
  not, and binaries for different versions of R were compiled with different
  versions of MinGW.  The entry point isnan is part of the statically linked
  runtime on modern MinGW.
 
 I used the binary on CRAN.

Then please use the compilers described in readme.packages which match it.

  MSVC++ 6.0 does supply _isnan (as it really should as it is part of the
  C99 ISO standard), and you need to link against it appropriately.  Hint:
  it may have an extra underline, since it seems it is known to C as _isnan.
  You may need to add
  
  #undef ISNAN
  #define ISNAN(x) _isnan(x)
 
 Many thanks for this suggestion. It works, as always. The extra
 underline is required.
 
  We don't support adding extensions to R using a different compiler to the 
  one used to build R.  Changes already made for future releases of R will 
  make this less likely to work in R 2.0.x.
 
 I wish that the R API entry points documented in Writing R Extensions
 be supported in the future versions of R, because it will make programs
 built with R more portable than directly using compiler-specific
 functions, such as testing and generating the IEEE 754 special values.
 The existing entry points have saved me trouble to chase these values
 whenever I moved to an environment with a different compiler.

Who said they would not supported?  What I said is that using a different
compilers to compile different parts of R is not supported (and never has
been).  And that is even less likely to work in future releases.  All we
guarantee is that ISNAN and R_FINITE produce calls to functions that work
under the system used to configure R.  As on Windows R_FINITE now produces
a call to a MinGW macro, that is not going to work under VC++ 6.

We are not talking about using `compiler-specific functions' here, like
_isnan.  What R core is doing is using *standard* C99 functions where
available, and a problem with VC++ 6.0 is that it is far from compliant
with IEC60559 aka IEEE754.  Given that the recommended compiler is freely 
available, it is unreasonable to expect any support for less capable 
compilers.

-- 
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
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] adress the current index

2004-05-24 Thread Mag. Ferri Leberl
How can I adress the current index of a vector?

I want to work with time series and therefore give the n-th element of a 
vector some value dependent on the value of the n-1th element.

Thank you in advance.

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] what does numeric mean? {was ... as.matrix.data.frame() ...}

2004-05-24 Thread Martin Maechler
 BDR == Prof Brian Ripley [EMAIL PROTECTED]
 on Mon, 24 May 2004 17:02:21 +0100 (BST) writes:

BDR I don't think a POSIXt element *is* numeric (that's a
BDR basic atomic vector), 

well, yes, in some strict sense, but at least

  is.atomic( ISOdatetime(2003,1,1:3,0,0,0))
 [1] TRUE

BDR so the new behaviour seems right to me.  The Warning is
BDR wrong, though, and will be fixed.

Further note that   data.matrix(.) still does convert to
numeric, and  AFAIK all good books on S / R  tell you to rather
use data.matrix(df) in situations you want a numeric matrix from a
data frame df.

BDR On Mon, 24 May 2004, Don MacQueen wrote:

 Conversion of a data frame to a matrix using as.matrix() when a 
 column of the data frame is POSIXt and all other columns are numeric 
 has changed in R 1.9.0 from R 1.8.1. The new behavior issues a 
 warning message and converts to a character matrix. In R 1.8.1, such 
 an object was converted to a numeric matrix.
 
 Here is an example.
 
  R 1.9.0 
  foo - data.frame( x=1:3,dt=ISOdatetime(2003,1,1:3,0,0,0))
 
  as.matrix(foo)
 x   dt 
 1 1 2003-01-01
 2 2 2003-01-02
 3 3 2003-01-03
 Warning message:
 longer object length is not a multiple of shorter object length in: 
 cl == c(Date, POSIXct, POSIXlt)


  


 From ?as.matrix() in R 1.9.0:
 
 'as.matrix' is a generic function. The method for data frames will
 convert any non-numeric/complex column into a character vector
 using 'format' and so return a character matrix, except that
 all-logical data frames will be coerced to a logical matrix.
 
 The POSIXt element is numeric, and so should be converted to numeric
   is.numeric(foo$dt)
 [1] TRUE

that is a point, particularly together with the above
is.atomic(foo$dt)  |- TRUE

 
 I think this might qualify for bug status, either in and of itself or 
 relative to documentation. But I'm not, as the posting guide says, 
 completely and utterly sure. So I'm posting to r-help first...I 
 will send a bug report if an R-core member asks me to.
 
 Thanks
 -Don

Note that with a factor (instead of POSIX*t),
things are similar

  str(ffoo - data.frame(x=1:3, f=gl(3,1)))
 `data.frame':  3 obs. of  2 variables:
  $ x: int  1 2 3
  $ f: Factor w/ 3 levels 1,2,3: 1 2 3
  as.matrix(ffoo)
   x   f  
 1 1 1
 2 2 2
 3 3 3
  data.matrix(ffoo)
   x f
 1 1 1
 2 2 2
 3 3 3

but in that case,

   is.numeric(ffoo$f)
  [1] FALSE

hence, all according to the docs

--

Maybe we should either be more specific in  help(as.matrix)
{and more other places ?} about what numeric means there,
or consider -- and this may well be unfeasible because it
potentially breaks current code -- to redefine
'is.numeric(.)' to be more restrictive:

e.g.,

   if (is.object(x))  ===   ! is.numeric(x) 

?

Martin

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] barplot

2004-05-24 Thread Carlos Guevel
I´ve tried version 1.9.0 barplot  with these (and others) example from  the
help page:

tN - table(Ni - rpois(100, lambda=5))
r - barplot(tN, col='gray') 
I get :


 ...OLE_Obj... 
Same example with version 1.8.1 gives the following result:


 ...OLE_Obj... 

What is wrong with v.1.9.0?

Thanks,

Carlos Guevel

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] identify() in script

2004-05-24 Thread Diogo FC Patrao
hallo,

I want a script of mine to run identify(), and possibly switch to
interactive mode at some point. Is it possible?

Thanks in advance

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] Applications of the Distance from a Point to a Curve

2004-05-24 Thread Dimas Martnez Morera
I would like to know if anyone knows about (real) situations where it is 
necessary to compute the euclidean distance from a point to a curve (or even 
the point where it is reached), you may not assume that the point is very 
close to the curve (unlike the usual fitting problems).
I have special interest in applications to phisics.

Thanks in advance,

Dimas

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] barplot

2004-05-24 Thread Sundar Dorai-Raj

Carlos Guevel wrote:
I´ve tried version 1.9.0 barplot  with these (and others) example from  the
help page:
tN - table(Ni - rpois(100, lambda=5))
r - barplot(tN, col='gray') 
I get :

 ...OLE_Obj... 
Same example with version 1.8.1 gives the following result:

 ...OLE_Obj... 

What is wrong with v.1.9.0?
Thanks,
Carlos Guevel
This has been reported a while back:
http://r-bugs.biostat.ku.dk/cgi-bin/R/Graphics?id=6777;expression=barplot;user=guest
which also suggests a workaround.
tN - table(Ni - rpois(100, lambda=5))
r - barplot(as.vector(tN), col='gray')
Alternatively, barplot2 in the gregmisc package should also work.
library(gregmisc)
tN - table(Ni - rpois(100, lambda=5))
r - barplot2(tN, col='gray')
--sundar
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] Tramo-seats

2004-05-24 Thread Paul Bivand
On Monday 24 May 2004 11:00, [EMAIL PROTECTED] wrote:
 Working  - among other things- in the field of (short  long term)
 electricity forecast,  we are now using too many  too expensive pieces of
 licensed software: SAS, SPSS, EViews. This sedimentation is due to the
 fact that my predecessors in the past used different consultant companies
 to manage each procedure.

 Having attended the useR2004! Conference with the aim of assessing if R
 ALONE could glue all those fragmented and isolated procedures, I'm almost
 convinced now that YES it could do the job!

 Now - to start with - a first problem to solve:

 we have to comply with the Tramo-seats closed-source procedure
 (http://www.bde.es/informes/be/docs/dt0014e.pdf) to deal with seasonality
 of electricity monthly time-series, in line with the methodology officially
 adopted by our  National Bureau of Statistics.

 Searching in R-help mailing list I didn't find anything about a tramo-seats
 R version.

 Does anyone know of unofficial R translation of tramo-seats?
 OR What do you suggest?

 Vittorio

 __
 [EMAIL PROTECTED] mailing list
 https://www.stat.math.ethz.ch/mailman/listinfo/r-help
 PLEASE do read the posting guide!
 http://www.R-project.org/posting-guide.html

Not an R response I'm afraid, but have you looked at Demetra 
http://forum.europa.eu.int/irc/dsis/eurosam/info/data/demetra.htm which 
appears to be an officially sponsored windows-only implementation of 
Tramo-seats and X-12 ARIMA. 

The licence is not clear from the website, but it looks free as in beer but 
probably not open source.

Paul Bivand

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] installing R on Fedora Core 2 test 2

2004-05-24 Thread Ulises Mora Alvarez
Hello:

Please give more info, how did you proceed?, when you install/upgrade 
Fedora what kind on instalation did you made? 

A few days ago I upgrade my computer from FC-1 to FC-2 (workstation). 
After that, I install the binaries (from CRAN) and everything went just 
fine.

On Mon, 24 May 2004, Benjamin Yakir wrote:

 I am new to Linux and just installed Fedora Core 2. Tried to install R from Fedora 
 core 1 binaries and from the source files according to the manual but failed. Would 
 appreciate any help.
 
 Benny Yakir
 Department of Statistics
 The Hebrew University
   [[alternative HTML version deleted]]
 
 __
 [EMAIL PROTECTED] mailing list
 https://www.stat.math.ethz.ch/mailman/listinfo/r-help
 PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
 

-- 
Ulises M. Alvarez
LAB. DE ONDAS DE CHOQUE
FISICA APLICADA Y TECNOLOGIA AVANZADA
UNAM
[EMAIL PROTECTED]

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] Seasonal ARIMA question - stat package (formerly ts)

2004-05-24 Thread Keith Campbell
To whom it may concern:
 
I am trying to better understand the functionality of 'R' when making
arima predictions to avoid any Black Box disadvantages.
 
I'm fitting a seasonal arima model using the following command (having
already loaded 'stat' package).
arimaSeason -
arima(Data,order=c(1,0,1),seasonal=list(order=c(1,0,1),period=12))
 
I can then generate subsequent predictions using the 'predict' function.
However, I can't seem to duplicate these predictions in a separate
program using the model coefficients.  From duplicating simpler models,
I understand the input variables are adjusted by the intercept term.
(e.g. for an arima(1,0,0) the prediction equation is   y(t) = beta1 * (
y(t-1) - beta0 )  +  beta0 where beta0 is the intercept)
 
Currently, I've expected the prediction to follow the equation below:
y(t)  =  beta0  +  [beta1*( y(t-1) - beta0 )]  +  [beta2 * epsilon(t-1)]
+  [beta3 * (y(t-12) - beta0)]  +  [beta4 * epsilon(t-12)]
 
This has proved unsuccessful.  What equation underlies this arima
prediction?  Is there something different that happens
 
Many thanks for your help,
Keith Campbell
 
 
Keith Campbell 
Researcher
ROTELLA CAPITAL MANAGEMENT |180 N. Stetson Suite 5100 I Chicago, IL
60601 |  t: 312.706.0442 |  f: 312.861.0468  | 

_

This message is for the designated recipient only and may contain
privileged, proprietary, or otherwise private information. If you are
not the intended recipient or otherwise believe that you have received
this message in error, please notify the sender immediately and delete
the original. Any other use of this message by you is prohibited.

 

[[alternative HTML version deleted]]

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] Metafiiles into Word R 1.9.0

2004-05-24 Thread Jason Turner
   2.-The computer at work has Microsoft Windows 2000 (5.00.2195 Service
 Pack
 2) and Word 2000 (9.0.4402 SR-1) I cannot copy-paste windows metafiles
 into
 Office applications in my computer at work. The resulting object is empty.


Is it genuinely empty, or can you see the graphs in print layout view, or
print preview?  I've seen that particular Word bug before.  It happens
with Visio drawings too (Visio is a M$ product).

Cheers

Jason

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] installing R on Fedora Core 2 test 2

2004-05-24 Thread Benjamin Yakir
Dear Prof Ripley

Thank you very much. Half the night is over but running the patch
worked.

Benny

On Tue, 2004-05-25 at 00:02, Prof Brian Ripley wrote:
 What was the failure?
 
 For the sources, there is a bug in the X11 headers it ships.  Please use 
 the R-patched tarball from ftp://ftp.stat.math.ethz.ch/Software/R/, which 
 has a workaround.
 
 On Mon, 24 May 2004, Benjamin Yakir wrote:
 
  I am new to Linux and just installed Fedora Core 2. Tried to install R
  from Fedora core 1 binaries and from the source files according to the
  manual but failed. Would appreciate any help.

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] Re: identify() in script

2004-05-24 Thread Diogo FC Patrao
On Mon, 2004-05-24 at 19:43, Diogo FC Patrao wrote:
 On Mon, 2004-05-24 at 17:17, Diogo FC Patrao wrote:
  hallo,
  
  I want a script of mine to run identify(), and possibly switch to
  interactive mode at some point. Is it possible?
  
  Thanks in advance

I think I didn`t explained my point properly. Consider I have a file
named, say, test.R, which contains

x=rnorm(100)
y=rnorm(100)
plot( x,y )
identify(x,y,1:100)

and want to invoke it as a batch, running (in linux)

$ R  test.R

just by doing this, no plot appear. That's my doubt: is it possible to
do this? or, there is a command which halts the execution at some point
and allow me to enter commands, as I've invoked R in interactive mode?

thanks for your patience,

dfcp

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] yeardays

2004-05-24 Thread Steve Su
Dear All,

Is there a R (Windows) equivalent of yeardays in Splus 6.0 for Windows? I am using XP.

**

 Steve Su ([EMAIL PROTECTED])
 PhD student.  

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


 Phone:  +61 7 3864 2017
 Fax:+61 7 3864 1812  
 Mobile: 0421  840  586  
 .   
   _--_|\  
  /  QUT   
  \_.--._/   
v   
  
** 
  
[[alternative HTML version deleted]]

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] Re: Windows versus Unix packages in CRAN ...

2004-05-24 Thread Duncan Murdoch
On 24 May 2004 17:18:52 +0200, Peter Dalgaard
[EMAIL PROTECTED] wrote:

Duncan Murdoch [EMAIL PROTECTED] writes:

 On Fri, 21 May 2004 15:41:18 +0200, Tamas Papp [EMAIL PROTECTED]
 wrote:
 
 You are probably right in saying that they _could_ have done better,
 but I would not use should in this context.  AFAIK the package is
 free (as in beer) software, which means that you are not paying for
 it.  The maintainers probably do not need a Linux version (yet), so it
 was not easy to use it under Linux.  Feel free to contribute.  
 
 Yes, indeed!  

Well, yes and no. Yes, people should feel free to help out with the
maintenance, but no, it is not reasonable to leave cross-platform
issues unaddressed.

There's a list of packages on CRAN in
/bin/windows/contrib/1.9/@ReadMe that fail to build on Windows for one
reason or another.   Should we critcize the authors of those packages
for not addressing cross-platform issues?  I don't think so.

Why should it be any different when the situation is reversed?

Duncan Murdoch

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] loess

2004-05-24 Thread Rupen Shrestha
Hi,
When I was running the function loess(y~x, span=0.0020), I got a warning 
message k-d tree limited by memory. ncmax= 4231

Does that mean the function has not been computed correctly ? If it has not, 
is there any way to adjust it so that it will do correctly ?

Thanks.
Rupen.




***
If you fail to plan, you are planning to fail.
_
Download music tracks from 95c here:
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] Scheire-Ray-Hare extension to Kruskal-Wallis

2004-05-24 Thread Sarah Skikne
Hi,

I was wondering if it was possible to do the Scheirer-Ray-Hare extension of the 
Kruskal-Wallis test using R, or if it was possible to use some other non-parametric 
test for the following situation:

I have continuous data that is grouped by 2 factor variables.  I would like a 
non-parametric test because the group variances are not equal.  I would like to look 
for the effects of either factor variable and their interaction.

Thanks for any help,
Sarah 
[[alternative HTML version deleted]]

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] e1071, R1.9.0, Solaris 2.9, should I be worried?

2004-05-24 Thread Richard A. O'Keefe
In R 1.9.0 running under Solaris 2.9 on a SunBlade 100,
with Sun WorkShop 6 update 2 C++ 5.3 2001/05/15 as the
C++ compiler, I just did
 install.packages(e1071)
The output includes these lines, which I have wrapped to fit nicely in mail:
** libs
cc -I/users/local/lib/R/include -I/usr/local/include  -KPIC -xlibmil \
 -dalign -xO4 -c cmeans.c -o cmeans.o
cc -I/users/local/lib/R/include -I/usr/local/include  -KPIC -xlibmil \
 -dalign  -xO4 -c cshell.c -o cshell.o
cc -I/users/local/lib/R/include -I/usr/local/include  -KPIC -xlibmil \
 -dalign -xO4 -c floyd.c -o floyd.o
cc -I/users/local/lib/R/include -I/usr/local/include  -KPIC -xlibmil \
 -dalign -xO4 -c Rsvm.c -o Rsvm.o
CC -I/users/local/lib/R/include -I/usr/local/include  -KPIC -xlibmil \
 -dalign -xO4 -c svm.cpp -o svm.o
svm.cpp, line 444: Warning: l hides Solver::l.
svm.cpp, line 444: Warning: Q hides Solver::Q.
svm.cpp, line 444: Warning: Cp hides Solver::Cp.
svm.cpp, line 444: Warning: Cn hides Solver::Cn.
svm.cpp, line 444: Warning: eps hides Solver::eps.
svm.cpp, line 507: Warning:
String literal converted to char* in formal argument fmt 
in call to info(char*, ...).
svm.cpp, line 517: Warning:
String literal converted to char* in formal argument fmt 
in call to info(char*, ...).
svm.cpp, line 690: Warning:
String literal converted to char* in formal argument fmt 
in call to info(char*, ...).
svm.cpp, line 881: Warning: l hides Solver::l.
svm.cpp, line 881: Warning: Q hides Solver::Q.
svm.cpp, line 881: Warning: b hides Solver::b.
svm.cpp, line 881: Warning: y hides Solver::y.
svm.cpp, line 881: Warning: alpha hides Solver::alpha.
svm.cpp, line 881: Warning: Cp hides Solver::Cp.
svm.cpp, line 881: Warning: Cn hides Solver::Cn.
svm.cpp, line 881: Warning: eps hides Solver::eps.
svm.cpp, line 881: Warning: si hides Solver_NU::si.
svm.cpp, line 1278: Warning:
String literal converted to char* in formal argument fmt 
in call to info(char*, ...).
svm.cpp, line 1328: Warning:
String literal converted to char* in formal argument fmt 
in call to info(char*, ...).
svm.cpp, line 1404: Warning:
String literal converted to char* in formal argument fmt 
in call to info(char*, ...).
svm.cpp, line 1439: Warning:
String literal converted to char* in formal argument fmt 
in call to info(char*, ...).
svm.cpp, line 1483: Warning:
String literal converted to char* in formal argument fmt 
in call to info(char*, ...).
svm.cpp, line 1507: Warning:
String literal converted to char* in formal argument fmt 
in call to info(char*, ...).
svm.cpp, line 1719: Warning:
String literal converted to char* in formal argument fmt 
in call to info(char*, ...).
24 Warning(s) detected.

How worried should I be?
I guess the Warning: x hides Solver::x warnings related to a deliberate
style choice, but what about the String literal converted to char* ones?

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html