Re: [R] Fitting a line to a qqplot's points?

2007-03-23 Thread Uwe Ligges
Paul Lynch wrote: I've made some normal plots of my data using qqplot, and now I would like to fit a line to the points on the plot and check the correlation coefficient to have a more objective measure of how straight the line is. Is there a simple way of doing that? (I'm still pretty new

Re: [R] Fitting a line to a qqplot's points?

2007-03-23 Thread Uwe Ligges
Uwe Ligges wrote: Paul Lynch wrote: I've made some normal plots of my data using qqplot, and now I would like to fit a line to the points on the plot and check the correlation coefficient to have a more objective measure of how straight the line is. Is there a simple way of doing that?

Re: [R] Truncated x-axis values

2007-03-23 Thread Stephen Tucker
You can try playing around with oma, omi, mai, mar, etc. in par(): myMai - par(mai) myMai[1] - max(nchar(y))*par(cin)[1] par(mfrow=c(2, 1),mai=myMai,oma=rep(0,4),las = 2) # your plot --- Urban, Alexander [EMAIL PROTECTED] wrote: John Thanks for your reply and sorry for not beein specific

[R] Im, Paseran and Shin test in R

2007-03-23 Thread justin bem
Dear all, Does R have a function to implement Im, Paseran and Shin test in R ? Justin BEM Elève Ingénieur Statisticien Economiste BP 294 Yaoundé. Tél (00237)9597295. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

[R] concatenate 2 data.frames

2007-03-23 Thread João Fadista
Dear all, I would like to know how can I concatenate 2 data.frames into a single one. Both data frames have the same number of columns and the same class type in each correspondent column. So what I want is to have a new data.frame where I have first the values from one data.frame and then

[R] Updating a worksheet in Excel file using RODBC

2007-03-23 Thread Moshe Olshansky
Hello! I have no problem reading Excel files (each worksheet in the file is a table which can be read - at least in my case). What I would like to do is to read such a table, change it (just the contents, not the format) and write it back, and this I can not do. I am getting the following

Re: [R] Updating a worksheet in Excel file using RODBC

2007-03-23 Thread Prof Brian Ripley
The problem is that way the ODBC driver exposes table names is not valid SQL, and nor is the way quoting has to be used. You can get around this via direct SQL sent by sqlQuery. In addition, by default the Excel ODBC driver gives you read-only access to worksheets. Searching the list

[R] landscape pdf

2007-03-23 Thread Sebastian Weber
Hello together! How can I plot a landscape letter-format plot? With postscript, I just use the horizontal option and I get what I want, but it seems that the pdf lacks this option. Well, I could do a ps2pdf conversion of the generated ps-file. But is there a way to directly produce landscape

Re: [R] landscape pdf

2007-03-23 Thread Uwe Ligges
Sebastian Weber wrote: Hello together! How can I plot a landscape letter-format plot? With postscript, I just use the horizontal option and I get what I want, but it seems that the pdf lacks this option. Well, I could do a ps2pdf conversion of the generated ps-file. But is there a way to

Re: [R] landscape pdf

2007-03-23 Thread Prof Brian Ripley
PDF does not have orientation: just set the width height in pdf() to get a 'landscape' plot. The 'horizontal' option in postscript() was originally (in S/S-PLUS) designed for sending plots direct to a printer: pdf() does not have that option either. On Fri, 23 Mar 2007, Sebastian Weber

Re: [R] non-linear curve fitting

2007-03-23 Thread Hufkens Koen
I'll give you the equation of the reference I based my thinking upon. this the link: http://users.pandora.be/requested/images/equation.png It's a sigmoid error function, but I thought I simplified things by picking a less complex form, although I think It won't matter that much. Anyway, as

Re: [R] concatenate 2 data.frames

2007-03-23 Thread Gavin Simpson
On Fri, 2007-03-23 at 08:51 +0100, João Fadista wrote: Dear all, I would like to know how can I concatenate 2 data.frames into a single one. Both data frames have the same number of columns and the same class type in each correspondent column. So what I want is to have a new data.frame

[R] Translating actxserver server commands to rcom

2007-03-23 Thread W Eryk Wolski
Dear R and R com user I have the following matlab which acesses a com server, code which I am trying to translate to R using the Rcom package: clear all; clear analysis; analysis = actxserver('EDAL.MSAnalysis');

Re: [R] Updating a worksheet in Excel file using RODBC

2007-03-23 Thread Prof Brian Ripley
On Fri, 23 Mar 2007, Prof Brian Ripley wrote: The problem is that way the ODBC driver exposes table names is not valid SQL, and nor is the way quoting has to be used. You can get around this via direct SQL sent by sqlQuery. In addition, by default the Excel ODBC driver gives you read-only

Re: [R] Fitting a line to a qqplot's points?

2007-03-23 Thread ian white
See also ?shapiro.test which is effectively based on this correlation coefficient. On Fri, 2007-03-23 at 00:13 -0400, Paul Lynch wrote: I've made some normal plots of my data using qqplot, and now I would like to fit a line to the points on the plot and check the correlation coefficient to

Re: [R] concatenate 2 data.frames

2007-03-23 Thread John Kane
?rbind rbind(df1, df2) This is a very basic question. You probably would save a lot of time if you read some of the documentation on the CRAN site (not the home R page). Click on Contributed on the left side of the screen and have a look at some of the documentation there. There are some

[R] Effect display of proportional odds model

2007-03-23 Thread Jan Wijffels
Dear useRs, I very much like the effect display of the proportional odds model on page 29 (Figure 8) of the following paper by John Fox: http://socserv.mcmaster.ca/jfox/Papers/logit-effect-displays.pdf It really gives a very concise overview of the model. I would like to use it to illustrate the

[R] distribution graph

2007-03-23 Thread Plessen, Christian von
I am looking for a way to produce a distribution graph as in the example: (http://cecsweb.dartmouth.edu/release1.1/datatools/dgraph.php?year=2003geotype=STD_HRRevent=A01_DISeventtype=UTIL Anybody who can help? Christian von Plessen Department of Pulmonary Medicine Haukeland university

Re: [R] Cohen's Kappa

2007-03-23 Thread Christian Schulz
...many thanks for all the answers and clarity!!! regards, christian Hi, Chance is not .5 in your data, it's a function of the expected values for presence and absence: (((7792*10855)/11974) + ((4182*1119)/11974))/11974 [1] 0.6225686 (.6862368-.6225686)/(1-.6225686)

[R] Creating new directory/folder from R script on run time.

2007-03-23 Thread d. sarthi maheshwari
Hi, Please correct me if this request does not belong to this discussion forum. Kindly suggest me if there is a way to create a new folder from R on runtime? for e.g. From one of my R script I want to store the PDF result file in a perticular folder let say c:\tt for that initially i am

Re: [R] Creating new directory/folder from R script on run time.

2007-03-23 Thread Henrique Dallazuanna
Perhaps you can try: ?dir.create dir.create('C:\tt') -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40 S 49° 16' 22 Ohttp://maps.google.com/maps?f=qhl=enq=Curitiba,+Brazillayer=ie=UTF8z=18ll=-25.448315,-49.276916spn=0.002054,0.005407t=kom=1 On 23/03/07, d. sarthi maheshwari [EMAIL

Re: [R] Truncated x-axis values

2007-03-23 Thread John Kane
Also think aboutpar( cex.axis = something small) Depending on what you're doing you might want to consider usin horizontal=TRUE and reduce the plot size a bit to accomomodate the label lengths. I don't know the packages well enough to be sure but I suspect that lattice or ggplot may handle this

[R] sort argument in mosaicplot

2007-03-23 Thread Jean lobry
Dear R-help, I do not understand how to use the sort argument in mosaicplot(). From the documentation sort is a vector ordering of the variables, containing a permutation of the integers 1:length(dim(x)) (the default). x - matrix(1:4,2,2) mosaicplot(x) # This one is OK mosaicplot(x, sort =

Re: [R] Effect display of proportional odds model

2007-03-23 Thread Gabor Grothendieck
Looks like there is code in the appendix. On 3/23/07, Jan Wijffels [EMAIL PROTECTED] wrote: Dear useRs, I very much like the effect display of the proportional odds model on page 29 (Figure 8) of the following paper by John Fox: http://socserv.mcmaster.ca/jfox/Papers/logit-effect-displays.pdf

[R] Binary information convert into hexadecimal

2007-03-23 Thread Schmitt, Corinna
Hallo, I just started with programming with R. I have the following problem: Given is binary information and should be translated into integer and afterwards into hexadecimal: 0 0 0001 1 1 0010 2 2 0011 3 3 0100 4 4 0101 5 5 0110 6 6 0111 7 7 1000 8 8 1001 9 9 1010 10

[R] dimension reduction and multinomial responses

2007-03-23 Thread Laurent Valdes
Hi, I'm trying to fit a pls model with a multinomial response. Do you have some ideas on the right package ? Laurent -- We are drowning in information, but starved for knowledge. «Germain» @http://www.le-valdo.com [[alternative HTML version deleted]]

[R] Get home directory and simple I/O

2007-03-23 Thread Alberto Monteiro
Is there any generic function that gets the home directory? This should return /home/user in Linux and x:/Documents and Settings/user (or whatever) in Windows XP. Another (unrelated) question: what is the _simplest_ way to read and write R variables to/from files such that they are stored in a

[R] Change Axis Size

2007-03-23 Thread webb
Sorry if this is an obvious question but I have not been able to find the answer. I wish to plot 3 lines on the same plot. However, whichever one I plot first, the axis does not have a big enough range for the other two to be shown in the plot (they get cut off at the top and the bottom). Is there

Re: [R] distribution graph

2007-03-23 Thread Michael Kubovy
?violinplot (You need to install the UsingR package first.) On Mar 23, 2007, at 4:06 AM, Plessen, Christian von wrote: I am looking for a way to produce a distribution graph as in the example: (http://cecsweb.dartmouth.edu/release1.1/datatools/dgraph.php?

Re: [R] Binary information convert into hexadecimal

2007-03-23 Thread Alberto Monteiro
Corinna Schmitt wrote: Given is binary information and should be translated into integer and afterwards into hexadecimal: 0 0 If you have a _string_ of 0s and 1s, and you want to convert it to an integer, I think the best way should be: (1) convert to a vector of 0s and 1s (2)

Re: [R] Effect display of proportional odds model

2007-03-23 Thread Ernst Hansen
Gabor Grothendieck writes: Looks like there is code in the appendix. The appendix even has a URL where the code is available, namely http://socserv.mcmaster.ca/jfox/Papers/polytomous-effect-displays.html On 3/23/07, Jan Wijffels [EMAIL PROTECTED] wrote: Dear useRs, I very much

Re: [R] Change Axis Size

2007-03-23 Thread ONKELINX, Thierry
You need to sed the ylim parameter in the plot function to the range that you need. Cheers, Thierry ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Reseach Institute for Nature and Forest Cel biometrie,

Re: [R] Get home directory and simple I/O

2007-03-23 Thread Andris Jankevics
If you want get a username of user currently running R on Linux,you can use a system command and read enviroment variables: paste(/home/,system (whoami,intern=TRUE),sep=) Andris Jankevics On Piektdiena, 23. Marts 2007 14:30, Alberto Monteiro wrote: Is there any generic function that gets the

Re: [R] Cohen's Kappa

2007-03-23 Thread Jonathan Baron
See also section 6.4 of http://www.psych.upenn.edu/~baron/rpsych/rpsych.html which also points to a few packages that have kappa code in them. -- Jonathan Baron, Professor of Psychology, University of Pennsylvania Home page: http://www.sas.upenn.edu/~baron

Re: [R] can't load just saved R object ReadItem: unknown type 65

2007-03-23 Thread Luke Tierney
According to the logs nothing at all has changed in the serialization code in a month and nothing of consequence for much longer than that. To track this down we will need a complete, reproducible, and preferably minimal example. Best, luke On Thu, 22 Mar 2007, Mark W Kimpel wrote: I have run

Re: [R] distribution graph

2007-03-23 Thread Ted Harding
On 23-Mar-07 11:06:49, Plessen, Christian von wrote: I am looking for a way to produce a distribution graph as in the example: (http://cecsweb.dartmouth.edu/release1.1/datatools/dgraph.php?year=2003; geotype=STD_HRRevent=A01_DISeventtype=UTIL Anybody who can help? Christian von

Re: [R] Change Axis Size

2007-03-23 Thread Gavin Simpson
On Fri, 2007-03-23 at 12:40 +, [EMAIL PROTECTED] wrote: Sorry if this is an obvious question but I have not been able to find the answer. I wish to plot 3 lines on the same plot. However, whichever one I plot first, the axis does not have a big enough range for the other two to be shown

Re: [R] Get home directory and simple I/O

2007-03-23 Thread Henrik Andersson
Sys.getenv(HOME) works in Linux at least - Henrik Andris Jankevics wrote: If you want get a username of user currently running R on Linux,you can use a system command and read enviroment variables: paste(/home/,system (whoami,intern=TRUE),sep=) Andris Jankevics On Piektdiena, 23. Marts

Re: [R] Effect display of proportional odds model

2007-03-23 Thread John Fox
Dear Jan, First, I inadvertently removed material on these displays from my web site when the paper was published in Sociological Methodology 2006. I'll update and repost the material some time in the next couple of days, including a copy of the published paper (with a link on my home page).

Re: [R] Get home directory and simple I/O

2007-03-23 Thread jim holtman
On 3/23/07, Alberto Monteiro [EMAIL PROTECTED] wrote: Is there any generic function that gets the home directory? This should return /home/user in Linux and x:/Documents and Settings/user (or whatever) in Windows XP. Another (unrelated) question: what is the _simplest_ way to read and write

Re: [R] Change Axis Size

2007-03-23 Thread John Kane
I think that technically the answer is no. What you need to do is set the size of the appropriate axis with the first plot command ( xlim = or ylim = Have a look at ?plot.default --- [EMAIL PROTECTED] wrote: Sorry if this is an obvious question but I have not been able to find the answer.

Re: [R] Effect display of proportional odds model

2007-03-23 Thread Jeff Miller
I REALLY found this paper to be helpful. Will you please let the list know once you have made the update? Thank you, Jeff Miller University of Florida -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of John Fox Sent: Friday, March 23, 2007 9:07 AM To:

Re: [R] Get home directory and simple I/O

2007-03-23 Thread Gabor Grothendieck
See: ?R.home ?dput On 3/23/07, Alberto Monteiro [EMAIL PROTECTED] wrote: Is there any generic function that gets the home directory? This should return /home/user in Linux and x:/Documents and Settings/user (or whatever) in Windows XP. Another (unrelated) question: what is the _simplest_

Re: [R] Get home directory and simple I/O

2007-03-23 Thread Liaw, Andy
From: Gabor Grothendieck See: ?R.home That's not what Alberto wanted: It gives the location of the R installation, not where user's home directory is. AFAIK Windows does not set the HOME environment variable by default. ?dput On 3/23/07, Alberto Monteiro [EMAIL PROTECTED] wrote: Is

[R] generating lognormal variables with given correlation

2007-03-23 Thread Mollet, Fabian
Dear R users I use simulated data to evaluate a model by sampling the parameters in my model from lognormal distributions. I would like these (lognormal distributed) parameters to be correlated, that is, I would like to have pairwise samples of 2 parameters with a given correlation

Re: [R] generating lognormal variables with given correlation

2007-03-23 Thread Karl Ove Hufthammer
Mollet, Fabian: I would like these (lognormal distributed) parameters to be correlated, that is, I would like to have pairwise samples of 2 parameters with a given correlation coefficient. I have seen that a covariance matrix can be fixed when generating random variables from a

Re: [R] Get home directory and simple I/O

2007-03-23 Thread Prof Brian Ripley
But the request was for a *generic* solution. On Windows there might not be anything corresponding to a home directory (and the rw-FAQ discusses the concept and how R resolves this). The best answer I know of is path.expand(~). On Fri, 23 Mar 2007, Henrik Andersson wrote: Sys.getenv(HOME)

[R] simulate from a multivariate lognormal distribution with fixed covariance/correlation structure

2007-03-23 Thread Mollet, Fabian
Dear R users I use simulated data to evaluate a model by sampling the parameters in my model from lognormal distributions. I would like these (lognormal distributed) parameters to be correlated, that is, I would like to have pairwise samples of 2 parameters with a given correlation

Re: [R] Get home directory and simple I/O

2007-03-23 Thread Gabor Grothendieck
On 3/23/07, Liaw, Andy [EMAIL PROTECTED] wrote: From: Gabor Grothendieck See: ?R.home That's not what Alberto wanted: It gives the location of the R installation, not where user's home directory is. AFAIK Windows does not set the HOME environment variable by default. ok. Try this,

Re: [R] Get home directory and simple I/O

2007-03-23 Thread Gabor Grothendieck
path.expand(~) also seems to work on Windows XP. On 3/23/07, Prof Brian Ripley [EMAIL PROTECTED] wrote: But the request was for a *generic* solution. On Windows there might not be anything corresponding to a home directory (and the rw-FAQ discusses the concept and how R resolves this). The

Re: [R] Get home directory and simple I/O

2007-03-23 Thread Alberto Monteiro
Prof Brian Ripley wrote: But the request was for a *generic* solution. On Windows there might not be anything corresponding to a home directory (and the rw-FAQ discusses the concept and how R resolves this). The best answer I know of is path.expand(~). Thanks, this is the only solution

[R] Logistic Regression

2007-03-23 Thread Sergio Della Franca
Dear R-Helpers, I want to perform a logistic regression on my dataset (y). I used the following code: logistic-glm(formula=interest_variable~.,family = binomial(link = logit), data = y) This run correctly. Then i want to develop the logistic regression with three different method: -forward

Re: [R] Effect display of proportional odds model

2007-03-23 Thread John Fox
Dear Jan, I've placed a copy of the published version of the paper at http://socserv.socsci.mcmaster.ca/jfox/Misc/polytomous-effect-displays/inde x.html, along with the R code for computing effects and their standard errors and R code for the graphs in the paper. As you'll see, the functions

[R] subtotal for same row data

2007-03-23 Thread Yuan, Qiaoping \(NIH/NIAAA\) [E]
Hi, There, I would like to subtotal the number in a specified column for all rows having the same data for specified columns. The following is the simple example: x=matrix(c(1,2,2,0.3,2,2,2,0.5,1,2,1,0.2),3,4,byrow=T) rownames(x)=c(R1,R2,R3) colnames(x)=c(C1,C2,C3,F) x C1 C2

Re: [R] Effect display of proportional odds model

2007-03-23 Thread Jan Wijffels
Great! That was all I needed! Best regards, Jan Jan Wijffels University Center for Statistics W. de Croylaan 54 3001 Heverlee Belgium tel: +32 (0)16 322784 fax: +32 (0)16 322831 http://www.kuleuven.be/ucs -Original Message- From: John Fox [mailto:[EMAIL PROTECTED] Sent: vrijdag 23

[R] subtotal for same row data

2007-03-23 Thread Yuan, Qiaoping \(NIH/NIAAA\) [E]
Hi, There, I would like to subtotal the number in a specified column for all rows having the same data for specified columns. The following is the simple example: x=matrix(c(1,2,2,0.3,2,2,2,0.5,1,2,1,0.2),3,4,byrow=T) rownames(x)=c(R1,R2,R3) colnames(x)=c(C1,C2,C3,F) x    C1 C2 C3   F R1 

[R] Simple bar plot question

2007-03-23 Thread Janet
Dear all, This is a simple question, but as far as I can tell, it's not in MASS or the R help archives or par man pages. Are there any cex.* parameters which allow me to set the text size of the labels for the axes? I want to plot the following table in which the Definitely not, etc. are

[R] plotting symbol

2007-03-23 Thread Federico Calboli
Hi All, can I have a plot where the symbol for the dots is smaller than pch =20 but bigger than pch = '.'? Best, Fede -- Federico C. F. Calboli Department of Epidemiology and Public Health Imperial College, St. Mary's Campus Norfolk Place, London W2 1PG Tel +44 (0)20 75941602 Fax +44

[R] Simple bar plot question

2007-03-23 Thread Janet
Dear all, This is a simple question, but as far as I can tell, it's not in MASS or the R help archives or par man pages. Are there any cex.* parameters which allow me to set the text size of the labels for the axes and the size of the legend? I want to plot the following table in which the

[R] Subset

2007-03-23 Thread Sergio Della Franca
Dear R-Helpers, I have this dataset: YEARPRODUCTS cluster 1 10 2 2 42 3 3 25 2 4 42 3 5 40 3 6 45 1 7 44 1 8 47 1

Re: [R] plotting symbol

2007-03-23 Thread Chuck Cleland
Federico Calboli wrote: Hi All, can I have a plot where the symbol for the dots is smaller than pch =20 but bigger than pch = '.'? Have you considered using the cex argument to reduce the size of pch=20? X - rnorm(20) par(mfrow=c(2,2)) plot(X, pch=20, cex=1.0) plot(X, pch=20, cex=0.6)

Re: [R] Get home directory and simple I/O

2007-03-23 Thread Peter Dalgaard
Andris Jankevics wrote: If you want get a username of user currently running R on Linux,you can use a system command and read enviroment variables: paste(/home/,system (whoami,intern=TRUE),sep=) This would be more to the point (and safer), I think. Sys.getenv(HOME) HOME /home/bs/pd

[R] simple bar plot question

2007-03-23 Thread Janet
Literally 60 seconds after I sent my question, I found the cex.names parameter to barplot. I haven't found a parameter for the size of the text in the legend. Apologies for clogging inboxes semi-unnecessarily. Janet __ R-help@stat.math.ethz.ch

Re: [R] Updating a worksheet in Excel file using RODBC

2007-03-23 Thread Hans-Peter
Hi, 2007/3/23, Moshe Olshansky [EMAIL PROTECTED]: Hello! I have no problem reading Excel files (each worksheet in the file is a table which can be read - at least in my case). What I would like to do is to read such a table, change it (just the contents, not the format) and write it back, and

Re: [R] R and clinical studies

2007-03-23 Thread Cody_Hamilton
Thanks for the tip. I will look forward to trying this package out soon! Regards, -Cody Hans-Peter [EMAIL PROTECTED]

Re: [R] distribution graph

2007-03-23 Thread Marc Schwartz
On Fri, 2007-03-23 at 14:22 +, [EMAIL PROTECTED] wrote: [Apologies -- there were errors in the code I posted previously. A corrected version is below] On 23-Mar-07 11:06:49, Plessen, Christian von wrote: I am looking for a way to produce a distribution graph as in the example:

Re: [R] generating lognormal variables with given correlation

2007-03-23 Thread Francisco J. Zagmutt
This reference may be relevant for you: Connover, W.J., Iman, R.L. A distribution-free approach to inducing rank correlation among input variables. Technometric, 3, 311-334, 1982. Also, you may want to look at a more modern approach implemented in the copula package: install.packages(copula)

[R] Course on Survival Analysis in Amsterdam

2007-03-23 Thread Dick Verkerk
Survival Analysis By Dr. Ronald Geskus May 24-25, 2007 Amsterdam, The Netherlands http://www.can.nl/events/details.php?id=30 This course is aimed at everyone who wants to analyze data in which the time to the occurrence of some event, and its dependence on covariates, is of interest. Survival

[R] course on Analysis of Repeated Measurements by Ronald Geskus

2007-03-23 Thread Dick Verkerk
Analysis of Repeated Measurements By Dr. Ronald Geskus April 27, 2007 Amsterdam, The Netherlands http://www.can.nl/events/details.php?id=31 This course is aimed at everyone who is working with data that contain repeated measurements on persons or otherwise related data and

Re: [R] plotting symbol

2007-03-23 Thread Peter Dalgaard
Chuck Cleland wrote: Federico Calboli wrote: Hi All, can I have a plot where the symbol for the dots is smaller than pch =20 but bigger than pch = '.'? Have you considered using the cex argument to reduce the size of pch=20? X - rnorm(20) par(mfrow=c(2,2)) plot(X, pch=20,

[R] using \t in mtext

2007-03-23 Thread Thomas Kaliwe
Hi, Using tab spaces in mtext e.g. mtext(\t) little squares are plotted. Is there a way to use \t without getting squares displayed? Thanks Thomas [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list

Re: [R] Subset

2007-03-23 Thread Mark Wardle
Sergio Della Franca wrote: Dear R-Helpers, I have this dataset: YEARPRODUCTS cluster 1 10 2 2 42 3 3 25 2 4 42 3 5 40 3 6 45 1 7 44

Re: [R] Subset

2007-03-23 Thread Uwe Ligges
See ?subset !!! And please do read the psoting guide. Uwe Ligges Sergio Della Franca wrote: Dear R-Helpers, I have this dataset: YEARPRODUCTS cluster 1 10 2 2 42 3 3 25 2 4 42 3 5

[R] creating R packs for all

2007-03-23 Thread Erin Hodgess
Dear R People: I am in the process of creating an R package via Windows. If I would decide to submit in to CRAN, what would I need to do in order to make it run for the Linux or Mac People, please? Thanks in advance! Sincerely, Erin Erin Hodgess Associate Professor Department of Computer and

Re: [R] plotting symbol

2007-03-23 Thread Chuck Cleland
Peter Dalgaard wrote: Chuck Cleland wrote: Federico Calboli wrote: Hi All, can I have a plot where the symbol for the dots is smaller than pch =20 but bigger than pch = '.'? Have you considered using the cex argument to reduce the size of pch=20? X - rnorm(20)

[R] R at the Language Log

2007-03-23 Thread davidr
I was pleasantly surprised in my visit to the Language Log to see R used: These days, the most straightforward way to answer a question like this is to do a Monte Carlo simulation. This is easy to do, e.g. in the (free software) statistics language R http://www.r-project.org/ :

[R] test ANOVA/ANCOVA

2007-03-23 Thread sebastien boutry
Hello everybody, I search a test for compare the k means but I have one quantitative variable and two groups date, traitement. And I suppose my samples are dependant with the date. What the statistical test would I use? Thank you. my data: dateEU DW wk1 EU1 5,324547829 wk1

[R] p-values for GLMMs

2007-03-23 Thread Cristina Gomes
Hi there, I have a question about the GLMM that I'm doing, that a statistician friend suggested I should have for my analysis. I would like to know if there's any way of obtaining a p value and R square for the full model (and not each variable separately) as to asses whether this model is

Re: [R] using \t in mtext

2007-03-23 Thread Marc Schwartz
On Fri, 2007-03-23 at 18:29 +0100, Thomas Kaliwe wrote: Hi, Using tab spaces in mtext e.g. mtext(\t) little squares are plotted. Is there a way to use \t without getting squares displayed? Thanks Thomas See this post:

Re: [R] Subset

2007-03-23 Thread Petr Klasterecky
You want to obtain a subset of your data, so what about to use subset()...?? You didn't even consider doing some basic search for the solution as the posting guide asks you... Petr Sergio Della Franca napsal(a): Dear R-Helpers, I have this dataset: YEARPRODUCTS cluster 1

Re: [R] creating R packs for all

2007-03-23 Thread Prof Brian Ripley
On Fri, 23 Mar 2007, Erin Hodgess wrote: Dear R People: I am in the process of creating an R package via Windows. If I would decide to submit in to CRAN, what would I need to do in order to make it run for the Linux or Mac People, please? If it passes R CMD check on Windows it should work

[R] Completely off topic, but amusing?

2007-03-23 Thread Bert Gunter
Folks: Thought that many on this list might find this amusing, perhaps even a bit relevant. Hope it's OK: WASHINGTON - The government's estimate of the number of Americans without health insurance fell by nearly 2 million Friday, but not because anyone got health coverage.

[R] Rweb - calculating mean for an excel table opened in Rweb

2007-03-23 Thread Natalie O'Toole
Hi, Would anyone happen to know the steps and code for calculating mean on one column in a multicolumn exel spreadsheet. I would like to do this calculation on a file opened from my computer into Rweb using the Rweb open file interface? Any help would be greatly appreciated! Thank-you, Nat

[R] memory, speed, and assigning results into new v. existing variable

2007-03-23 Thread David L. Van Brunt, Ph.D.
I have a very large data frame, and I'm doing a conversion of all columns into factors. Takes a while (thanks to folks here though, for making faster!), but am wondering about optimization from a memory perspective... Internally, am I better off assigning into a new data frame, or doing one of

Re: [R] memory, speed, and assigning results into new v. existing variable

2007-03-23 Thread David L. Van Brunt, Ph.D.
Yes, at the end... but I'm getting memory allocation errors and outright crashes DURING the operation when assigning into the same variable... still waiting for the alternative to survive or not... On 3/23/07, Bert Gunter [EMAIL PROTECTED] wrote: Actually... there can be more than one working

[R] objects of class matrix and mode list?

2007-03-23 Thread Stephen Tucker
Hello everyone, I cannot seem to find information about objects of class matrix and mode list, and how to handle them (apart from flattening the list). I get this type of object from using sapply(). Sorry for the long example, but the code below illustrates how I get this type of object. Is

Re: [R] simple bar plot question

2007-03-23 Thread Stephen Tucker
I think you can set legend=FALSE in barplot() and add your own legend, in which you have a lot more control: barplot(#your arguments#,legend=FALSE) legend(x=topleft,cex=yourCex) etc. --- Janet [EMAIL PROTECTED] wrote: Literally 60 seconds after I sent my question, I found the cex.names

Re: [R] distribution graph

2007-03-23 Thread Ted Harding
On 23-Mar-07 16:55:40, Marc Schwartz wrote: [...] How about something like this: DistPlot - function(x, digits = 1, ...) { x - round(x, digits) Tab - table(x) Vals - sapply(Tab, function(x) seq(x) - mean(seq(x))) X.Vals - unlist(Vals, use.names = FALSE) tmp -

Re: [R] simple bar plot question

2007-03-23 Thread Stephen Tucker
Sorry, that was a bit premature - you probably want other arguments in legend as well; in particular 'fill' is an argument you'd probably be interested in for legends of barplots: legend(x=topleft,cex=yourCex,fill=yourColors,legend.text=yourText, lty=NA,pch=NA,#and so on#) --- Stephen

Re: [R] can't load just saved R object ReadItem: unknown type 65

2007-03-23 Thread Mark W Kimpel
Luke, I'll be gone for about 2 weeks but will work on getting you a reproducible example when I get back. If this topic comes up with anyone else, please copy me on your responses as I may miss it in the 600 emails I'll have to delete on my return :) Mark Luke Tierney wrote: According to the

[R] Space for X and Y axis labels

2007-03-23 Thread Joseph Wakeling
Hello all, I'm having a bit of a problem with x and y axis labels. Two things: first, if I want to create a plot with, plot.new() plot.window(.) axis(1) axis(2) lines(...) points(...) [etc.] ... where do I introduce the xlab=... and ylab=... commands? I

[R] [R-pkgs] RTisean 3.0-7 released

2007-03-23 Thread Antonio, Fabio Di Narzo
Dear R users, I've just uploaded to CRAN a new version of RTisean, the TISEAN-to-R interface. This is now compatible with the recent, new 3.0.1 release of TISEAN [1]. This new TISEAN version is explicitely GPL-ed, and has some more routines handling multivariate time series. Bests, Antonio, Fabio

Re: [R] subtotal for same row data

2007-03-23 Thread jim holtman
try this: x - as.data.frame(x) x C1 C2 C3 F R1 1 2 2 0.3 R2 2 2 2 0.5 R3 1 2 1 0.2 do.call('rbind',by(x, list(x$C1, x$C2), function(z){z$F - sum(z$F); z[1,]})) C1 C2 C3 F R1 1 2 2 0.5 R2 2 2 2 0.5 On 3/23/07, Yuan, Qiaoping (NIH/NIAAA) [E] [EMAIL PROTECTED] wrote:

[R] Conversion from string to date type

2007-03-23 Thread Andreas Tille
Hi, I'm on my very first steps to R and so I hope that I do not ask a really stupid questions but I did not found it via R-Search, in the FAQ or Google (BTW, the name R is not a really good seekiong criterion ;-) ). I have a data file containing a table that containes dates and values like

[R] plotting dnorm() issued from mclust models

2007-03-23 Thread Frederic Jean
Dear all I have a problem in fitting lines() of the normal distributions identified with Mclust on a histogram or a mclust1Dplot. Here is some sample code to explain : set.seed(22) foo - c(rnorm(400, 10, 2), rnorm(500, 17, 4)) mcl - Mclust(foo, G=2) mcl.sd -

Re: [R] Conversion from string to date type

2007-03-23 Thread jim holtman
Here is how you can convert them to a Date object: x - c('01.03.2007','02.03.2007','03.03.2007') y - as.Date(x, format=%d.%m.%Y) y [1] 2007-03-01 2007-03-02 2007-03-03 str(y) Class 'Date' num [1:3] 13573 13574 13575 On 3/23/07, Andreas Tille [EMAIL PROTECTED] wrote: Hi, I'm on my

[R] Optimal matching analysis package

2007-03-23 Thread Julien Barnier
Hi, I recently heard about application of sequencing algorithms to social data, which seems very interesting. There seems to be at least two programs whith whom you can do this kind of analysis : Optimize and TDA. http://home.uchicago.edu/~aabbott/om.html

Re: [R] Conversion from string to date type

2007-03-23 Thread Andreas Tille
On Fri, 23 Mar 2007, jim holtman wrote: Here is how you can convert them to a Date object: x - c('01.03.2007','02.03.2007','03.03.2007') y - as.Date(x, format=%d.%m.%Y) y Well, this is what I tried when reading the docs, but mydata - read.csv(file='mydata.dat', sep = '\t', quote='',

Re: [R] Conversion from string to date type

2007-03-23 Thread Gabor Grothendieck
Read the help desk article in R-News 4/1 and see ?as.Date ?strptime (for setting the as.Date format= argument) Also, you might be interested in the zoo package library(zoo) ?read.zoo vignette(zoo) vignette(zoo-quickref) On 3/23/07, Andreas Tille [EMAIL PROTECTED] wrote: Hi, I'm on my very

Re: [R] Conversion from string to date type

2007-03-23 Thread Dirk Eddelbuettel
Hi Andreas, Welcome to r-help :) On 23 March 2007 at 22:21, Andreas Tille wrote: | I'm on my very first steps to R and so I hope that I do | not ask a really stupid questions but I did not found it | via R-Search, in the FAQ or Google (BTW, the name R is | not a really good seekiong criterion

  1   2   >