Re: [R] High resolution graphics from R Windows Vista

2008-07-31 Thread Prof Brian Ripley
On the TIFF point -- are you using 2.7.1 patched? It is likely that the crash is fixed there (it was intermittent and some of us never see it). Otherwise you have not described what is 'substandard' nor given any examples, and others are not reporting such differences. (I do have XP and

Re: [R] multiple separators in scan() functions

2008-07-31 Thread Prof Brian Ripley
On Thu, 31 Jul 2008, milton ruser wrote: Dear all, Is there a way of I deal with more than one kind of separators when reading a asc file using scan() function? The separators could be :, \t and ;. Yes, pre-process the file (e.g. use a pipe connection and 'tr', or read the file with

Re: [R] stats question

2008-07-31 Thread Moshe Olshansky
Hello Jason, You are not specific enough. What do you mean by significant difference? Let's assume that indeed the incidence in A is 6% and in B is 10% and we are looking for Na and Nb such that with probability of at least 80% the mean of Nb sample from B will be at least, say, 0.03 (=3%)

[R] add string

2008-07-31 Thread Alfredo Alessandrini
Hi, How can I join two string? frequency = 15 I want join the number frequency with a string. Alfredo __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] add string

2008-07-31 Thread john seers (IFR)
Is this what you want: paste(Mystring, frequency, sep=) --- -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Alfredo Alessandrini Sent: 31 July 2008 10:11 To: r-help@r-project.org Subject: [R] add string Hi, How can I join two string? frequency =

Re: [R] bug in 'margins' behavior in reshape - cast

2008-07-31 Thread Michael Frumin
thanks, that workaround, well, works! what are you working on with the rewrite? just efficiency? or major changes in functionality/interface? thanks, mike On Wed, Jul 30, 2008 at 7:10 PM, hadley wickham [EMAIL PROTECTED] wrote: On Wed, Jul 30, 2008 at 12:41 PM, mfrumin [EMAIL PROTECTED]

[R] OT:Selling Data Mining Cloud-A Concept

2008-07-31 Thread Ajay ohri
The Ohri Framework - Data Mining on Demandhttp://decisionstats.com/2008/the-ohri-framework-data-mining-on-demand/ Part of the reason SAS continues to enjoy a profitable lead is 1) Standardized language elements (Data and Procs) 2)Ease of Learning SAS 3) Output Delivery to multiple sources 4)

[R] combinations with replications

2008-07-31 Thread MarinaTarantini
Dear all, Is there a way to compute and list all combinations with replication of two elements in sets of 8 elemnts? For example, I've two elements, 0 and 1, and I would to get all possible combinations with replication such as, for example, , 0001, 0010, and so on. They are 2^8

[R] random effects mixed model, different regressors

2008-07-31 Thread Fensis
Hi everybody, I have built a model that includes subject ID as a random effect, and has a continous variable (time) and I want to test whether the slope of this line differs between treatments (this is tested with the interaction between treatment and time). My question now is that I also want

[R] Identifying common prefixes from a vector of words, and delete those prefixes

2008-07-31 Thread Daren Tan
For example, c(dog.is.an.animal, cat.is.an.animal, rat.is.an.animal). How can I identify the common prefix is .is.an.animal and delete it to give c(dog, cat, rat) ? Thanks _ [[alternative HTML version deleted]]

Re: [R] combinations with replications

2008-07-31 Thread Jim Lemon
On Thu, 2008-07-31 at 02:29 -0700, MarinaTarantini wrote: Dear all, Is there a way to compute and list all combinations with replication of two elements in sets of 8 elemnts? For example, I've two elements, 0 and 1, and I would to get all possible combinations with replication such as, for

Re: [R] combinations with replications

2008-07-31 Thread ONKELINX, Thierry
Or simply Datagrid - expand.grid(rep(list(0:1), 8)) apply(Datagrid, 1, paste, collapse = ) HTH, Thierry ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest Cel

[R] nls weights warning message

2008-07-31 Thread Lisa Readdy
The following warning message occurs when running the nls on some data: ~ Warning message: In is.na(wts) : is.na() applied to non-(list or vector) of type 'NULL' ~ and I would like to know what is

Re: [R] odds ratios in multiway tables (stratified)

2008-07-31 Thread Wim Bertels
On Wed, 2008-07-30 at 12:13 -0500, Marc Schwartz wrote: on 07/30/2008 08:48 AM Wim Bertels wrote: Hi, does anyone know of a function to calculate odds ratios in multiway tables (stratified) (+ the other usual statistics involved) i mean: say we have a table r*c*d, For every d

Re: [R] bug in 'margins' behavior in reshape - cast

2008-07-31 Thread hadley wickham
On Thu, Jul 31, 2008 at 4:59 AM, Michael Frumin [EMAIL PROTECTED] wrote: thanks, that workaround, well, works! what are you working on with the rewrite? just efficiency? or major changes in functionality/interface? I've come up with a strategy that makes reshape _much_ faster - 10-20x for

Re: [R] Sampling two exponentials

2008-07-31 Thread Zhang Yanwei - Princeton-MRAm
Thank you very much. Would you talk more about that? How can I use the copula package to sample two dependent exponentials? Which function shall I use? Sincerely, Yanwei Zhang Department of Actuarial Research and Modeling Munich Re America Tel: 609-275-2176 Email: [EMAIL PROTECTED]

Re: [R] nls weights warning message

2008-07-31 Thread Prof Brian Ripley
On Thu, 31 Jul 2008, Lisa Readdy wrote: The following warning message occurs when running the nls on some data: ~ Warning message: In is.na(wts) : is.na() applied to non-(list or vector) of type 'NULL'

Re: [R] bug in 'margins' behavior in reshape - cast

2008-07-31 Thread mfrumin
sounds great. there are often times I'd like to use reshape on data frames of hundreds of thousands or millions of rows, but I have found that it is just too slow at this point to be convenient. thanks again for everything, Mike hadley wrote: On Thu, Jul 31, 2008 at 4:59 AM, Michael

Re: [R] stats question

2008-07-31 Thread Marc Schwartz
At the risk of oversimplifying the study design, this sounds like a two sample comparison of proportions, in which case power.prop.test() would be the function of interest. This could also be done via Monte Carlo simulation, which would not be difficult to implement. Note that I am also

Re: [R] Sampling two exponentials

2008-07-31 Thread Ben Bolker
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Zhang Yanwei - Princeton-MRAm wrote: | Thank you very much. | Would you talk more about that? How can I use the copula package to sample two dependent exponentials? Which function shall I use? | ~ It's a big topic, and I'm not very familiar with

Re: [R] Identifying common prefixes from a vector of words, and delete those prefixes

2008-07-31 Thread Richard . Cotton
For example, c(dog.is.an.animal, cat.is.an.animal, rat.is.an. animal). How can I identify the common prefix is .is.an.animal and delete it to give c(dog, cat, rat) ? foo - c(dog.is.an.animal, cat.is.an.animal, rat.is.an.animal) sub(.is.an.animal, , foo) Being pedantic, .is.an.animal is a

Re: [R] nls weights warning message

2008-07-31 Thread Prof Brian Ripley
Somewhat coincidentally this will work if you don't supply 'data=Data'. On Thu, 31 Jul 2008, Prof Brian Ripley wrote: On Thu, 31 Jul 2008, Lisa Readdy wrote: The following warning message occurs when running the nls on some data: ~ Warning

Re: [R] equivalent R functions for Numerical Recipes fitxy and fitexy ?

2008-07-31 Thread Liaw, Andy
Not a direct answer to your questions, but for error-in-variables problems, there are newer technologies than what is in NR. For example: install.packages(simex) library(simex) example(simex) Andy From: Marc Fischer Dear Folks, We need to fit the model y~x assuming there are random

Re: [R] Converting to subscripts and superscripts

2008-07-31 Thread Ben Bolker
naw3 at duke.edu writes: Hi, I am reading in a CSV file of chemical reactions where the subscripts and superscripts are encoded in angle brackets, like below: 2 HSUP+/SUP + 2 OSUB2/SUBSUP-/SUP Is there a way to convert these to actual sub/superscripts and save them in another excel

[R] Identifying Points That Belongs to Components of Mixture Models

2008-07-31 Thread Gundala Viswanath
Hi all, Suppose I have the following points: x [1] 1.376545e+00 1.017594e+00 6.734179e+00 2.537501e+00 9.614635e-02 [6] 7.078693e-03 4.067031e-02 1.527854e+00 6.795155e-03 6.247076e-04 [11] 1.934793e-06 1.442637e-02 8.432116e-02 1.126674e-03 6.476260e+00 [16] 1.472769e-01 2.247548e+00

Re: [R] Coefficients of Logistic Regression from bootstrap - how to get them?

2008-07-31 Thread Michal Figurski
Frank and all, The point you were looking for was in a page that was linked from the referenced page - I apologize for confusion. Please take a look at the two last paragraphs here: http://people.revoledu.com/kardi/tutorial/Bootstrap/examples.htm Though, possibly it's my ignorance, maybe

Re: [R] History pruning

2008-07-31 Thread Ken Williams
On 7/30/08 1:59 PM, Marc Schwartz [EMAIL PROTECTED] wrote: I (and many others) use ESS (Emacs Speaks Statistics), in which case, I have an R source buffer in the upper frame and an R session in the lower frame. I also use ESS to edit my R code (inside Aquamacs Emacs), but I usually use the

Re: [R] Significance of confidence intervals in the Non-Linear Least Squares Program.

2008-07-31 Thread bbouling
Dear Pf. Ripley, I use the nls() function coupled with the confint() to obtain the confidence interval of assessed parameters of fitting model. I don't understand why the confidence intervals estimated with such a method aren't consistent? Regards B. Boulinguiez Ph.D. in Chemistry Ecole

Re: [R] stats question

2008-07-31 Thread Iasonas Lamprianou
Amazing responses to ,y question, and very fast, thanks all of you!  Dr. Iasonas Lamprianou Department of Education The University of Manchester Oxford Road, Manchester M13 9PL, UK Tel. 0044 161 275 3485 [EMAIL PROTECTED] - Original Message From: Marc Schwartz [EMAIL PROTECTED] To:

Re: [R] tensor product of equi-spaced B-splines in the unit square

2008-07-31 Thread Patrizio Frederic
dear all, I apologize for a second post on the same subject. I still have the problem. I'm going to describe the problem in a different setting: I have 3 matrix A,B,W such that dim(A) = c(n,M) dim(B) = c(n,M) dim(W) = c(M,M) what I'm searching for is an efficient computation of vector R,

Re: [R] Coefficients of Logistic Regression from bootstrap - how to get them?

2008-07-31 Thread Gustaf Rydevik
On Thu, Jul 31, 2008 at 4:30 PM, Michal Figurski [EMAIL PROTECTED] wrote: Frank and all, The point you were looking for was in a page that was linked from the referenced page - I apologize for confusion. Please take a look at the two last paragraphs here:

Re: [R] odds ratios in multiway tables (stratified)

2008-07-31 Thread Marc Schwartz
on 07/31/2008 07:40 AM Wim Bertels wrote: On Wed, 2008-07-30 at 12:13 -0500, Marc Schwartz wrote: on 07/30/2008 08:48 AM Wim Bertels wrote: Hi, does anyone know of a function to calculate odds ratios in multiway tables (stratified) (+ the other usual statistics involved) i mean: say we have

Re: [R] History pruning

2008-07-31 Thread Marc Schwartz
on 07/31/2008 08:35 AM Ken Williams wrote: On 7/30/08 1:59 PM, Marc Schwartz [EMAIL PROTECTED] wrote: I (and many others) use ESS (Emacs Speaks Statistics), in which case, I have an R source buffer in the upper frame and an R session in the lower frame. I also use ESS to edit my R code

[R] Code to calculate internal rate of return

2008-07-31 Thread Thomas E
Hi all. I am an R newbie and trying to grasp how the simple optimization routines in R work. Specifically, I would like some guidance on how to set up a code to calculate the internal rate of return on an investment project (http://en.wikipedia.org/wiki/Internal_rate_of_return). My main problem

[R] Sperman Correlation with rcorr (Hmisc)

2008-07-31 Thread Birgitle
Hello R-User! I have a data.frame with 82 variables (columns) and 290 rows. The variables are set to classes factor, ordered factor and numeric. I used the following code Matrix.My.data-as.matrix(Df.My.Data[2:82]) Matrix.My.data.rcorr-rcorr(Matrix.My.data, type=spearman) and got the

Re: [R] question about krige code in R

2008-07-31 Thread Roger Bivand
Alessandro alessandro.montaghi at unifi.it writes: Hi All, I have this question. I wish to create a kriging map with R but I haven't a auxiliary map. I have only one txt file with X, Y and Z records. When I use this code: ... subground.ok - krige (Z~sqrt(X+Y), subground,

Re: [R] Identifying common prefixes from a vector of words, and delete those prefixes

2008-07-31 Thread John Kane
There MUST be a better way but this will work. x - c(dog.is.an.animal, cat.is.an.animal, rat.is.an.animal) bb - strsplit(x, \\.) myfun - function(m) m[1] animals - unlist(lapply(bb, myfun)) animals --- On Thu, 7/31/08, Daren Tan [EMAIL PROTECTED] wrote: From: Daren Tan [EMAIL PROTECTED]

[R] Random number generation

2008-07-31 Thread Max
Hi Everyone, I did a quick search of the list and it looks like this may not have been asked before... I'm trying to generate a matrix of random numbers between 0 and 1, with 6 columns, 1 rows. About all I know is that runif(1) gives me the random number I'm looking for. Any help would

Re: [R] Random number generation

2008-07-31 Thread Marc Schwartz
on 07/31/2008 12:24 PM Max wrote: Hi Everyone, I did a quick search of the list and it looks like this may not have been asked before... I'm trying to generate a matrix of random numbers between 0 and 1, with 6 columns, 1 rows. About all I know is that runif(1) gives me the random number

Re: [R] Coefficients of Logistic Regression from bootstrap - how to get them?

2008-07-31 Thread Michal Figurski
Gustaf, Summarizing things I don't understand: - Honestly, I was thinking I can use bootstrap to obtain better estimate of a mean - provided that I want it. So, I can't? - If I can't obtain reliable estimates of CI and variance from a small dataset, but I can do it with bootstrap - isn't it a

Re: [R] Random number generation

2008-07-31 Thread Max
Marc, this is very handy. My next question is, do you know a quick and easy way to transfer all of the output to a txt file? (or .xls)? Thanks, -Max Marc Schwartz explained on 07/31/2008 : on 07/31/2008 12:24 PM Max wrote: Hi Everyone, I did a quick search of the list and it looks like

Re: [R] Sampling two exponentials

2008-07-31 Thread Zhang Yanwei - Princeton-MRAm
I am trying to do this using the copula library and find a possible way out. library(copula) x=mvdc(claytonCopula(.75),c(exp,exp),list(list(rate=1),list(rate=2))) x.sample=rmvdc(x,100) The above code gives a sample with two marginal exponential ditributions. But what does the first argument

Re: [R] Random number generation

2008-07-31 Thread Jorge Ivan Velez
Hi Max, See ?write.table. Perhaps: MAT - matrix(runif(1 * 6), 1, 6) # TXT format write.table(MAT, C:/yourmatrix.txt,col.names=FALSE,row.names=FALSE,quote=FALSE) # XLS format write.table(MAT, C:/yourmatrix.xls,col.names=FALSE,row.names=FALSE,quote=FALSE,sep=\t) HTH, Jorge On Thu,

Re: [R] Random number generation

2008-07-31 Thread Patrizio Frederic
?write.table could help PF 2008/7/31 Max [EMAIL PROTECTED]: Marc, this is very handy. My next question is, do you know a quick and easy way to transfer all of the output to a txt file? (or .xls)? Thanks, -Max Marc Schwartz explained on 07/31/2008 : on 07/31/2008 12:24 PM Max wrote:

[R] anisotrophy in fit.variogram problem

2008-07-31 Thread Alessandro
Hi All, I wish to create a directional variograms at (0, 45, 90, 135) degrees from north (y-axis). v = variogram(Z~1, subground, cutoff=1800, width=80, alpha=c(0, 45, 90, 135)) 90 and 45 deg. look good but I don't understand in the fit.variogram how to create a ansi=c(X) with 90

[R] How to output R image to a file?

2008-07-31 Thread Aiste Aistike
Hello, I would like to ask if anyone could help me. I want to save images I create (e.g. histograms, boxplots, plots, etc.) to a file or files. Does anyone know how to do this? Thank you. Aiste [[alternative HTML version deleted]] __

Re: [R] Create an R package

2008-07-31 Thread FScottDahlgren
Thanks so much for the link. Scott Gabor Grothendieck wrote: Vista works fine for developing packages. You might look at http://batchfiles.googlecode.com which has some batch files for use with Vista and R. On Thu, Jul 24, 2008 at 9:13 PM, FScottDahlgren [EMAIL PROTECTED] wrote:

Re: [R] Random number generation

2008-07-31 Thread David Barron
m - runif(6) mm - matrix(m, ncol=6) HTH Dave On Thu, Jul 31, 2008 at 6:24 PM, Max [EMAIL PROTECTED] wrote: Hi Everyone, I did a quick search of the list and it looks like this may not have been asked before... I'm trying to generate a matrix of random numbers between 0 and 1, with 6

Re: [R] History pruning

2008-07-31 Thread Ken Williams
On 7/31/08 11:01 AM, hadley wickham [EMAIL PROTECTED] wrote: I think that would be a very hard task - Well, at least medium-hard. But I think significant automatic steps could be made, and then a human can take over for the last few steps. That's why I was enquiring about tools rather than

Re: [R] Re gression between adjacent columns - error with NAs

2008-07-31 Thread rcoder
Thank you both very much for your assistance. Both the suggestions worked out in the end, and I managed to achieve what I wanted. There is something else I want to try, which is a slight deviation on the theme. In the code I posted, I export Intercept and Slope regression coefficients to an o/p

Re: [R] How to output R image to a file?

2008-07-31 Thread Erik Iverson
?Devices for a start. Aiste Aistike wrote: Hello, I would like to ask if anyone could help me. I want to save images I create (e.g. histograms, boxplots, plots, etc.) to a file or files. Does anyone know how to do this? Thank you. Aiste [[alternative HTML version deleted]]

Re: [R] Random number generation

2008-07-31 Thread Max
This is perfect, thanks! :) Jorge Ivan Velez pretended : Hi Max, See ?write.table. Perhaps: MAT - matrix(runif(1 * 6), 1, 6) # TXT format write.table(MAT, C:/yourmatrix.txt,col.names=FALSE,row.names=FALSE,quote=FALSE) # XLS format write.table(MAT,

Re: [R] How to output R image to a file?

2008-07-31 Thread Aiste Aistike
Hi, Thank everyone for you help. I got my histogram written to a file. Thanks once again. Aiste 2008/7/31 Jian Zhang [EMAIL PROTECTED] Hi aiste, two easy example produce mypic.ps or mypic.pdf for you x-rnorm(50) postscript(mypic.ps) hist(x) dev.off() pdf( mypic.pdf) hist(x)

Re: [R] History pruning

2008-07-31 Thread Duncan Murdoch
On 7/31/2008 2:08 PM, Ken Williams wrote: On 7/31/08 11:01 AM, hadley wickham [EMAIL PROTECTED] wrote: I think that would be a very hard task - Well, at least medium-hard. But I think significant automatic steps could be made, and then a human can take over for the last few steps. That's

[R] 'system' command through Rserve

2008-07-31 Thread Patil, Prasad
Hello, I've installed an Rserve instance on a remote server. I have no problem interfacing with it and running most commands. I have loaded some R scripts on the remote server, and one of them contains a system command to copy a file (created by the script) onto another server. When I run the

[R] bwplot with Date object

2008-07-31 Thread GOUACHE David
Hello R-helpers, I would like to produce a boxplot for dates, using lattice. Here is a dummy example : dates-as.Date(32768:32895,origin=1900-01-01) plouf-data.frame(days=dates,group=factor(rep(1:2,times=128/2))) bwplot(group~days,data=plouf) # doesn't work, whereas :

[R] vertical bpplot

2008-07-31 Thread GOUACHE David
Hello all, I have a question on bpplot (or more specifically panel.bpplot used with bwplot) from package Hmisc. I would like to produce a box-percentile plot but vertically. And I hqve not been able to use the horizontal argument... using one of the examples from ?panel.bpplot : bwplot(g

Re: [R] Query about finding correlations

2008-07-31 Thread Stefan Grosse
I am sorry if I misinterpreted this but this questions looks/looked very like a very basic one about statistical inference... (have a look at the Manuals/contributed documentation e.g. Faraway) Lalitha Viswanath wrote: Hi This is not a homework assignment :) Me and my manager are trying to

[R] Randomly drawing observations from factors.

2008-07-31 Thread Economics Guy
I have a large data set where one of the columns needs be a unique identifier (ID) for each row. However for a few of the rows they have the same ID. What I need to do is randomly draw one of the rows and keep it in the data frame and drop all the others which have the same ID. For example: v1 -

Re: [R] upgrade to 2.5

2008-07-31 Thread Sundar Dorai-Raj
Iasonas Lamprianou said the following on 5/2/2007 8:25 AM: Hi I am using R version 2.4.1. How can I upgrade to version 2.5 without having to install all the packages again? Thanks Jason You may find the following link relevant.

[R] rollapply() to portions of a matrix

2008-07-31 Thread rcoder
Hi everyone, I have a rollapply statement that applies a function, in steps, over a data matrix as follows: #Code start testm-rollapply(mat, 100, by=100, min, na.rm=F) #Code end This moves down matrix 'mat' and calculates the minimum value over a 100 row range, every 100 rows (i.e. no

Re: [R] History pruning

2008-07-31 Thread Ken Williams
On 7/31/08 2:12 PM, Duncan Murdoch [EMAIL PROTECTED] wrote: expression() returns a list of language objects, and we only asked for one. We can look inside it: Hey, cool. Now let me see if I can do anything useful with that. Thanks. -Ken -- Ken Williams Research Scientist The

Re: [R] reference in article

2008-07-31 Thread Peter Dalgaard
Tomas Mikoviny wrote: Hi all R positive, does anyone know how to refer R in article? thanks tomas Like this (well, almost! This is from R 2.4.1 - 2007 on the most recentversion) citation() To cite R in publications use: R Development Core Team (2006). R: A language and environment for

Re: [R] Query about finding correlations

2008-07-31 Thread Stefan Grosse
How about making your homeworks yourselfes? lalitha viswanath wrote: Hi I have a dataframe which has 3 columns of numeric data A,B,C each of which has been obtained independent of the other. We are trying to find out, which of A or B cause C i.e. We are hypothesising that C is the effect

[R] instal tar.gz package on windows

2008-07-31 Thread georgiana onicescu
Dear HelpeRs, I have created a tar.gz package (the package was created on unix) and I would like to install it in R on a Windows operating system. As far as I know, R in Windows accepts only packages in zip file for installation and I could not find a way to install a tar.gz file. I would like

[R] strip names lattice graphics

2008-07-31 Thread GOUACHE David
R-helpers, I'm having difficulty with customizing strip names for a lattice graphic. Here is an example using the iris data set : xyplot(Sepal.Length+Sepal.Width~Petal.Length,groups=Species,data=iris) ## I'd like to change the 2 strip names to Length and Width for example, this is what

Re: [R] Randomly drawing observations from factors.

2008-07-31 Thread Marc Schwartz
on 07/31/2008 03:20 PM Economics Guy wrote: I have a large data set where one of the columns needs be a unique identifier (ID) for each row. However for a few of the rows they have the same ID. What I need to do is randomly draw one of the rows and keep it in the data frame and drop all the

Re: [R] upgrade to 2.5

2008-07-31 Thread Gabor Grothendieck
On 5/2/07, Sundar Dorai-Raj [EMAIL PROTECTED] wrote: Iasonas Lamprianou said the following on 5/2/2007 8:25 AM: Hi I am using R version 2.4.1. How can I upgrade to version 2.5 without having to install all the packages again? Thanks Jason You may find the following link relevant.

Re: [R] Coarsening the Resolution of a Dataset

2008-07-31 Thread Steve Murray
Please find below my command inputs, subsequent outputs and errors that I've been receiving. crops - read.table(crop2000AD.asc, colClasses = numeric, na=-) str(crops[1:10]) 'data.frame': 2160 obs. of 10 variables: $ V1 : num NA NA NA NA NA NA NA NA NA NA ... $ V2 : num NA NA NA NA

Re: [R] History pruning

2008-07-31 Thread hadley wickham
It's not really equivalent, natural language has ambiguities and subtleties that computer languages, especially functional languages, intentionally don't have. By their nature, computer languages can be turned into parse trees unambiguously and then those trees can be manipulated. But in

Re: [R] Is R's fast fourier transform function different from fft2 in Matlab?

2008-07-31 Thread Li Li
Thanks for both replies. Then I found the ifft2 from Matlab gives different result from fft( , inverse=T) from R. An example: in R: temp - matrix(c(1,4,2, 20), nrow=2) fft(temp) [,1] [,2] [1,] 27+0i -17+0i [2,] -21+0i 15+0i fft(temp,inverse=T) [,1] [,2] [1,] 27+0i -17+0i

Re: [R] instal tar.gz package on windows

2008-07-31 Thread Marc Schwartz
on 07/31/2008 03:07 PM georgiana onicescu wrote: Dear HelpeRs, I have created a tar.gz package (the package was created on unix) and I would like to install it in R on a Windows operating system. As far as I know, R in Windows accepts only packages in zip file for installation and I could not

Re: [R] Problem with the installation of install R on Sun Solaris

2008-07-31 Thread Prof Brian Ripley
You are missing 'ranlib': take a closer look at the output from configure (and the posting guide). My guess is that either /usr/ccs/bin is not in your path or you don't have the tools installed to build applications. This is not an R issue, and you need to seek local Solaris help. On Thu, 3

Re: [R] strip names lattice graphics

2008-07-31 Thread Deepayan Sarkar
On Thu, Jul 31, 2008 at 1:30 PM, GOUACHE David [EMAIL PROTECTED] wrote: R-helpers, I'm having difficulty with customizing strip names for a lattice graphic. Here is an example using the iris data set : xyplot(Sepal.Length+Sepal.Width~Petal.Length,groups=Species,data=iris) ## I'd like to

Re: [R] 'system' command through Rserve

2008-07-31 Thread Jeffrey Horner
This probably has to do with your ssh configuration and nothing to do with R. How are you starting Rserve? Is it run with the same user privileges as when you run R manually? Best, Jeff Patil, Prasad wrote on 07/31/2008 02:51 PM: Hello, I've installed an Rserve instance on a remote

[R] allocMatrix limits

2008-07-31 Thread Vadim Kutsyy
I am getting an error allocMatrix: too many elements specified when I am trying to create large matrix or vector (about 1 billion elements). How can I find out limits on allocMatrix? Can I increase them? Thanks, Vadim. PS: I am running R on SUSE10 on

[R] anisotropy in vgm model. HELP!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

2008-07-31 Thread Alessandro
Hi All, I have this problem. I don’t understand the right code in R when I have an anisotropy in the semivariogram model plot(variogram(Z~1, subground, cutoff=1800, width=80, alpha=c(45, 135, 90, 135))) I have a good model in 90° and eventually in 90° and 45° v = variogram(Z~1,

[R] sort rows of matrix by rows of another matrix

2008-07-31 Thread Timothy W. Hilton
Hello all, I am trying to sort rows of one matrix by rows of another. Given a1 and a2: -- a1 [,1] [,2] [,3] [1,]768 [2,]424 [3,]472 [4,]038 a1 - structure(c(7, 4, 4, 0, 6, 2, 7, 3, 8, 4, 2, 8), .Dim = c(4L, 3L)) a2 [,1] [,2] [,3]

[R] Help with hazard plots

2008-07-31 Thread Alan Cox
Hello.  I am hoping someone will be willing to help me understand something about hazard plots created with muhaz(...).  I have some background in statistics (minor in grad school), but I haven't been able to figure one thing about hazard plots.  I am using hazard plots to track customer

Re: [R] rollapply() to portions of a matrix

2008-07-31 Thread Gabor Grothendieck
rollapply over an index. Modify this appropriately: z - zoo(101:109, 11:19) f - function(ix) if (ix[1] 5) min(z[ix]) else max(z[ix]) rollapply(zoo(seq_along(z)), 3, by = 3, f) 2 5 8 101 104 109 On Thu, Jul 31, 2008 at 3:36 PM, rcoder [EMAIL PROTECTED] wrote: Hi everyone, I have a

[R] dput vs unclass to see what a factor really is composed of

2008-07-31 Thread Jacob Wegelin
I used read.dta() to read in a Stata 9 dataset to R. The Sex01 variable takes on two values in Stata: 0 and 1, and it is labeled M and F respectively, analogous to an R factor. Thus, read.dta reads it in as a factor. Now, I wanted to see what this variable *really* is, in R. For instance,

Re: [R] instal tar.gz package on windows

2008-07-31 Thread Stefan Evert
Dear HelpeRs, I have created a tar.gz package (the package was created on unix) and I would like to install it in R on a Windows operating system. As far as I know, R in Windows accepts only packages in zip file for installation and I could not find a way to install a tar.gz file. I would

Re: [R] Help with hazard plots

2008-07-31 Thread Alan Cox
I clicked Send before making sure I thanked anyone who took the time to help me out.  Sorry about that.  To all who read or respond: thanks!. - Original Message - From: Alan Cox [EMAIL PROTECTED] To: r-help@r-project.org Sent: Thursday, July 31, 2008 5:29:35 PM GMT -05:00 US/Canada

Re: [R] allocMatrix limits

2008-07-31 Thread Prof Brian Ripley
On Thu, 31 Jul 2008, Vadim Kutsyy wrote: I am getting an error allocMatrix: too many elements specified when I am trying to create large matrix or vector (about 1 billion elements). How can I find out limits on allocMatrix? Can I increase them? ?Memory-limits, and you cannot increase them

Re: [R] dput vs unclass to see what a factor really is composed of

2008-07-31 Thread Prof Brian Ripley
See ?NumericConstants. On Thu, 31 Jul 2008, Jacob Wegelin wrote: I used read.dta() to read in a Stata 9 dataset to R. The Sex01 variable takes on two values in Stata: 0 and 1, and it is labeled M and F respectively, analogous to an R factor. Thus, read.dta reads it in as a factor. Now, I

Re: [R] instal tar.gz package on windows

2008-07-31 Thread Duncan Murdoch
On 31/07/2008 4:07 PM, georgiana onicescu wrote: Dear HelpeRs, I have created a tar.gz package (the package was created on unix) and I would like to install it in R on a Windows operating system. As far as I know, R in Windows accepts only packages in zip file for installation and I could not

Re: [R] sort rows of matrix by rows of another matrix

2008-07-31 Thread Johnson, Eric A. (Seattle)
If you're not adverse to cbind-ing a1 and a2, you can use this: a1a2 - cbind(a1, a2) a3 - t(apply(a1a2, 1, function(x) x[order(x[1:ncol(a1)])+ncol(a1)])) Eric -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Timothy W. Hilton Sent: Thursday, July 31,

Re: [R] dput vs unclass to see what a factor really is composed of

2008-07-31 Thread Duncan Murdoch
On 31/07/2008 5:21 PM, Jacob Wegelin wrote: I used read.dta() to read in a Stata 9 dataset to R. The Sex01 variable takes on two values in Stata: 0 and 1, and it is labeled M and F respectively, analogous to an R factor. Thus, read.dta reads it in as a factor. Now, I wanted to see what this

Re: [R] bwplot with Date object

2008-07-31 Thread Deepayan Sarkar
On Thu, Jul 31, 2008 at 12:54 PM, GOUACHE David [EMAIL PROTECTED] wrote: Hello R-helpers, I would like to produce a boxplot for dates, using lattice. Here is a dummy example : dates-as.Date(32768:32895,origin=1900-01-01) plouf-data.frame(days=dates,group=factor(rep(1:2,times=128/2)))

Re: [R] sort rows of matrix by rows of another matrix

2008-07-31 Thread markleeds
below is another way ( maybe the same ? )but with an extra line to make roworder. i'm also not clear on why I have to take the transpose of the result that comes back form the apply call. in ?apply, it says that the function tries to convert the result back to a matrix. that's fine but why

[R] multiple comparison

2008-07-31 Thread Agus Susanto
Dear all, I was trying to understand how multcomp package works by running the examples given in the documentation. However I still don't understand when it comes to multiple comparison set by user (please refer to Ksub in the code). Therefore I run 2 other cases along with the original example

Re: [R] allocMatrix limits

2008-07-31 Thread Vadim Kutsyy
I am getting an error allocMatrix: too many elements specified when I am trying to create large matrix or vector (about 1 billion elements). How can I find out limits on allocMatrix? Can I increase them? ?Memory-limits, and you cannot increase them unless you have a system which has larger

[R] clustering and data-mining...

2008-07-31 Thread losemind
Hi all, I am doing some experiment studies... It seems to me that with different combination of 5 parameters, the end results ultimately converged to two scalars. That's to say, some combinations of the 5 parameters lead to one end result and some other combinations of the 5 parameters lead to

[R] P-I curve

2008-07-31 Thread Bruno P Masi
I am trying to write a function for P-I curve meansurament of Platt (1980). I used the exemple of nonlinear least squares in R-intro.pdf but I can't finish the analysis. Could you help me? Is there another way to do this analysis? Thanks. I- c(0,0,100,100,200,200,500,500) P-

Re: [R] sort rows of matrix by rows of another matrix

2008-07-31 Thread Steven McKinney
I came up with the same solution as Mark Leeds. a1roworders - t(apply(a1, 1, order)) a2ord - t(sapply(seq(nrow(a1)), function(x) a2[x, a1roworders[x,]] )) a2ord [,1] [,2] [,3] [1,] 102 101 103 [2,] 102 101 103 [3,] 103 101 102 [4,] 101 102 103 Mark's question about the

[R] How to abort JGR

2008-07-31 Thread g-h-d
Hi, say I have a loop that cannot complete. how do I abort within JGR? I use unbuntu linux. Thanks. -- View this message in context: http://www.nabble.com/How-to-abort-JGR-tp18766719p18766719.html Sent from the R help mailing list archive at Nabble.com.

[R] Storing Matrices into Hash

2008-07-31 Thread Gundala Viswanath
Hi, Suppose I have these two matrices (could be more). What I need to do is to store these matrices into a hash. So that I can call back any of the matrix back later. Is there a way to do it? mat_1 [,1][,2] [1,] 9.327924e-01 0.067207616 [2,] 9.869321e-01 0.013067929

Re: [R] Storing Matrices into Hash

2008-07-31 Thread Erik Iverson
I think a named list is probably the easiest way to start off, something like: all_mat - list(mat1 = mat1, mat2 = mat2) all_mat$mat2 Gundala Viswanath wrote: Hi, Suppose I have these two matrices (could be more). What I need to do is to store these matrices into a hash. So that I can

Re: [R] Storing Matrices into Hash

2008-07-31 Thread Gundala Viswanath
Thanks so much Eric, But how do you include that in a loop. I tried this, doesn't seem to work. Please advice: __BEGIN__ all_mat - NULL for (matno in 1:10) { mat - process_to_create_matrix(da[matno]) all_mat - list(all_post, matno = mat) } print(all_mat) # it gives funny structure.

Re: [R] Coarsening the Resolution of a Dataset

2008-07-31 Thread jim holtman
In the first case, 'crops' is an object that is 142MB in size, and that alone will take up almost 50% of your memory. You are also probably generating another object of like size in doing the operation. Can you cut down the size of the objects that you need. What operating system do you have and

  1   2   >