Re: [R] How to read plain text documents into a vector?

2009-10-14 Thread Dieter Menne
Richard Liu wrote: There are actually two vignettes. Both have examples of a vector of characters being made into a tm corpus, but neither shows how to read documents on the file system into the vectors. I tried the other two suggestions, but paste seemed not to glue the separate lines

[R] Error from termplot() with make.panel.svysmooth() for complex survey data

2009-10-14 Thread Christopher T. Moore
Greetings, I am using library(survey) to analyze some complex sample data. After fitting a model I tried to use termplot() with make.panel.svysmooth(), but I received an error (see below). Could someone help me interpret the error message so I can make the necessary corrections? The

Re: [R] attach

2009-10-14 Thread Dieter Menne
Christophe Dutang1 wrote: I would like to know what happens on the memory side when I use attach(inputdata) Is there a second allocation of memory for inputdata? Not, it just guides the syntax. Christophe Dutang1 wrote: Is it better not to use attach function? A qualified

Re: [R] reference on permutation test

2009-10-14 Thread Meyners, Michael, LAUSANNE, AppliedMathematics
It's always worthwhile to look at the articles by Pitman (and maybe the textbook by Fisher, if you have access to it); Welch is a nice paper, too, but might be pretty technical to learn about the area. I don't know any of the textbooks except Edgington (which is in its 4th edition now with

Re: [R] heatmap.2

2009-10-14 Thread Dieter Menne
John Celniker wrote: ... when I try to change the breaks to: br2 [1] 0.000 0.5337751 1.0675502 1.6013253 2.1351003 3.000 3.500 4.000 4.500 [10] 4.8039758 I get the correct heatmap representation but the color key does not update correctly to reflect changes in

Re: [R] update.formula drop interaction terms

2009-10-14 Thread Dieter Menne
Eleni Rapsomaniki-2 wrote: How do I drop multiplication terms from a formula using update? e.g. forml=as.formula(Surv(time, status) ~ x1+x2+A*x3+A*x4+B*x5+strata(sex)) The easiest way is to write the formula again without the A's. Dieter -- View this message in context:

[R] post-hoc test with kruskal.test()

2009-10-14 Thread Robert Kalicki
Dear R users, I would like to know if there is a way in R to execute a post-hoc test (factor levels comparison, like Tukey for ANOVA) of a non-parametric analysis of variance with kruskal.test() function. I am comparing three different groups. The preliminary analysis using the

Re: [R] histogram

2009-10-14 Thread Jim Lemon
On 10/13/2009 10:06 PM, Dmitry Gospodaryov wrote: Dear R developers,How I can build a histogram from matrix: 0 0.5 1 0.25 34 43 65 1 23 35 54 4 22 29 42 10 21 22 29 20 15 17 20 (first string is represented names of columns, first column is represented names of rows) where names of columns

Re: [R] svy / weighted regression

2009-10-14 Thread Laust
Dear Thomas David That makes sense! If I wanted to use survey on the summarized data, I suppose that I could 'de-summarize' or 're-individualize' the data to give the design object the correct information on the number of observations. Or I could revert to using the actual individual-level data.

Re: [R] S4 tutorial

2009-10-14 Thread Robin Hankin
Peng the Brobdingnag package includes a vignette that gives a step-by-step guide to creating a simple package that uses S4. best wishes Robin Peng Yu wrote: I'm looking for some tutorial on S4. I only find the following one, which is not in English. Can somebody let me know if there is any

Re: [R] How to read plain text documents into a vector?

2009-10-14 Thread Dieter Menne
kenhorvath wrote: Dieter Menne wrote: library(tm) filenames = list.files(path=.,pattern=\\.txt) docs = for (filename in filenames){ docs = c(docs,paste(readLines(file(filename)),collapse=\n)) } docs ## continue as in example vs = VectorSource(docs) If in any way

Re: [R] post-hoc test with kruskal.test()

2009-10-14 Thread Meyners, Michael, LAUSANNE, AppliedMathematics
Robert, you can do the corresponding paired comparisons using wilcox.test. As far as I know, there is no such general correction as Tukey's HSD for the Kruskal-Wallis-Test. However, if you have indeed only 3 groups (resulting in 3 paired comparisons), the intersection-union principle and the

Re: [R] update.formula drop interaction terms

2009-10-14 Thread Uwe Ligges
Eleni Rapsomaniki wrote: Dear R users, How do I drop multiplication terms from a formula using update? e.g. forml=as.formula(Surv(time, status) ~ x1+x2+A*x3+A*x4+B*x5+strata(sex)) #I would like to drop all instances of variable A (the main effect and its interactions). The following:

Re: [R] Request update on A (Not So) Short Introduction to S4

2009-10-14 Thread Uwe Ligges
Peng Yu wrote: There are several '?'s on the last page of the following document. Apparently, they are not correct. Could somebody correct it? cran.r-project.org/doc/contrib/Genolini-S4tutorialV0-5en.pdf Please ask the author. Best, Uwe Ligges

[R] Linear Regression Question

2009-10-14 Thread Alexandre Cohen
Dear Sir or Madam, I am a student at MSc Probability and Finance at Paris 6 University/ Ecole Polytechnique. I am using R and I can't find an answer to the following question. I will be very thankful if you can answer it. I have two vectors rendements_CAC40 and rendements_AlcatelLucent. I use

Re: [R] R commander.

2009-10-14 Thread Ilyas .
thank u so much for ur help,,i have tried it before,,and succeded to do so,,the problem is that, in the next stage of the package i am facing a problem which i dont know how to fix,,,i am useing cnvpack ( http://www.meb.ki.se/~yudpaw/) and in the following command i get an error,,which is..

[R] ggplot2 scale_shape question

2009-10-14 Thread Jonathan Bleyhl
Is there a way to have some points solid and some points hollow? I have two classes of points and there are so many points, that it's hard to see just the difference in shapes. I'd like to have one of the classes be hollow in addition to being a different shape. Any help would be grand. Thanks,

Re: [R] How to read plain text documents into a vector?

2009-10-14 Thread kenhorvath
Richard Liu wrote: I tried the other two suggestions, but paste seemed not to glue the separate lines together into one character string. Perhaps I missed something (collapse?). Perhaps I'll have another look. Yes, that is what 'collapse' should do! If you read text using

Re: [R] How to read plain text documents into a vector?

2009-10-14 Thread kenhorvath
Dieter Menne wrote: While I agree that the appending could be more efficiently be done by a list as an intermediate, the docs = c(doc, ljljljl) construct is not recursive, even if not efficient. Yes, of course, that was hastily written, sorry ... but from my experience list

[R] default borders in boxplot and barplot

2009-10-14 Thread Jennifer Young
This is my first post so hopefully I haven't mucked up the rules. I'm trying to change the default borders in either boxplot or barplot so that, at the request of a journal, all of my figures have the same type of border. I've successfully used par(bty=o) using plot(1:10, bty=o), but it seems

[R] currency conversion function?

2009-10-14 Thread Liviu Andronic
Dear all Is there any R function that would perform currency conversion using up-to-date exchange rates? I would be looking for a function that allows to download recent exchange rates (say, from Yahoo!) and then use these in converting currencies (say, USD to EUR). I am not sure whether

Re: [R] default borders in boxplot and barplot

2009-10-14 Thread Barry Rowlingson
On Wed, Oct 14, 2009 at 2:21 AM, Jennifer Young jennifer.yo...@math.mcmaster.ca wrote: This is my first post so hopefully I haven't mucked up the rules. I'm trying to change the default borders in either boxplot or barplot so that, at the request of a journal, all of my figures have the same

Re: [R] Linear Regression Question

2009-10-14 Thread Jakson A. Aquino
On Tue, Oct 13, 2009 at 11:17:11PM +0200, Alexandre Cohen wrote: I have two vectors rendements_CAC40 and rendements_AlcatelLucent. I use the lm function as follows, and then the sumarry function: regression=lm(rendements_CAC40 ~ rendements_AlcatelLucent); sum=summarry(regression); [...] I

Re: [R] Linear Regression Question

2009-10-14 Thread Ted Harding
On 13-Oct-09 21:17:11, Alexandre Cohen wrote: Dear Sir or Madam, I am a student at MSc Probability and Finance at Paris 6 University/ Ecole Polytechnique. I am using R and I can't find an answer to the following question. I will be very thankful if you can answer it. I have two vectors

Re: [R] General means of matching a color specification to an official R color name

2009-10-14 Thread Jim Lemon
Hi Bryan, You can get a near match with the color.id function in the plotrix package. Jim __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] General means of matching a color specification to an official R color name

2009-10-14 Thread Romain Francois
On 10/14/2009 12:05 AM, Barry Rowlingson wrote: On Tue, Oct 13, 2009 at 10:58 PM, Bryan Hansonhan...@depauw.edu wrote: Works perfectly! Thanks Barry. I had actually seen some suggestions on using a distance, but by then I was thinking about hcl spaces and distance isn't so as simple there.

Re: [R] General means of matching a color specification to an official R color name

2009-10-14 Thread Romain Francois
On 10/14/2009 12:26 PM, Romain Francois wrote: On 10/14/2009 12:05 AM, Barry Rowlingson wrote: On Tue, Oct 13, 2009 at 10:58 PM, Bryan Hansonhan...@depauw.edu wrote: Works perfectly! Thanks Barry. I had actually seen some suggestions on using a distance, but by then I was thinking about hcl

Re: [R] rect function

2009-10-14 Thread Rene
Dear all, I have a question about how to load data or (entering data )to each cell of a rectangle created by rect . e.g. I have a matrix rbind(1:2,1:2) I have created a 2x2 rectangle by using: a-0:1/10 b-0:1/10 kk-expand.grid(a,b) plot.new() rect(kk[, 1], kk[, 2], kk[, 1] +

Re: [R] rect function

2009-10-14 Thread Jim Lemon
On 10/14/2009 10:20 PM, Rene wrote: Dear all, I have a question about how to load data or (entering data )to each cell of a rectangle created by rect . e.g. I have a matrix rbind(1:2,1:2) I have created a 2x2 rectangle by using: a-0:1/10 b-0:1/10 kk-expand.grid(a,b) plot.new()

[R] RCMD Rdconv debugging output

2009-10-14 Thread Erich Neuwirth
I am trying (on Windows XP, with R 2.10.0beta) to use RCMD Rdconv -t html myfile.Rd myfile.html to convert some Rd files to html. I get a message that there are warning. How can I tell Rdconv to show me these warnings? -- Erich Neuwirth, University of Vienna Faculty of Computer Science Computer

Re: [R] Selecting initial numerals

2009-10-14 Thread joris meys
On Tue, Oct 13, 2009 at 6:48 PM, PDXRugger j_r...@hotmail.com wrote: II just want to create a new object with the first two numerals of the data. Not sure why this isnt working, consider the following: EmpEst$naics=c(238321, 624410, 484121 ,238911, 81, 531110, 621399, 541613, 524210

Re: [R] RCMD Rdconv debugging output

2009-10-14 Thread Duncan Murdoch
On 10/14/2009 7:45 AM, Erich Neuwirth wrote: I am trying (on Windows XP, with R 2.10.0beta) to use RCMD Rdconv -t html myfile.Rd myfile.html to convert some Rd files to html. I get a message that there are warning. How can I tell Rdconv to show me these warnings? You can do the same

Re: [R] currency conversion function?

2009-10-14 Thread Henrique Dallazuanna
You can try something about like this: foo - function(from, to, date){ url - http://www.oanda.com/convert/classic?script=..%2Fconvert%2Fclassiclanguage=envalue=1; params - sprintf(%sdate=%sexch=%sexch2=margin_fixed=0expr=%sexpr2=SUBMIT=Convert+Nowlang=endate_fmt=us, url,

Re: [R] cdf

2009-10-14 Thread David Winsemius
On Oct 13, 2009, at 6:53 PM, Duncan Murdoch wrote: On 13/10/2009 6:43 PM, David Winsemius wrote: On Oct 13, 2009, at 5:12 PM, maram salem wrote: Dear all, I have the cdf of the following power fuction distribution: F(y)=(y/350)^a ,0y350, where a is some parameter with range

[R] plot discriminant analysis

2009-10-14 Thread Alejo C.S.
I'm confused on how is the right way to plot a discriminant analysis made by lda function (MASS package). (I had attached my data fro reproduction). When I plot a lda object : X - read.table(data, header=T) lda_analysis - lda(formula(X), data=X) plot(lda_analysis) #the above plot is completely

[R] Strange characters that block import

2009-10-14 Thread arnaud Mosnier
Dear useRs, I try to import a text file that contain some strange characters coming from the misinterpretation of foreign language characters by another software (see below). Here is an example of text with a line containing characters that bug the import

Re: [R] post-hoc test with kruskal.test()

2009-10-14 Thread David Winsemius
There is a post hoc test along the lines of the Kruskal-Wallis test. It is implemented on the help page of oneway_test from package coin. The authors of the package, Hothorn, Hornik, van de Wiel, and Zeileis, cite Hollander and Wolfe (1999) for details and say it is called the

[R] pairs

2009-10-14 Thread Seyit Ali Kayis
Dear all, I have two sets of data (say set1 and set2) as follow: set1 x1 x2 x3 0.30 0.43 3.88 0.38 0.59 3.53 0.30 0.42 2.12 0.33 0.53 2.12 0.30 0.47 3.76 set2 y1 y2 y3 0.32 0.47 5.18 0.23 0.26 1.06 0.42 0.65 3.88 0.28 0.38 3.76 0.35 0.47 1.41 The pairs

Re: [R] post-hoc test with kruskal.test()

2009-10-14 Thread Meyners, Michael, LAUSANNE, AppliedMathematics
Robert, What do you mean by not symmetric? If you mean unbalanced in terms of sample size, that's not a problem if you choose the right specifications for wilcox.test. The Kruskal-Wallis-Test is a generalization of the unpaired Wilcoxon test for more than two groups. Not sure whether

Re: [R] RCMD Rdconv debugging output

2009-10-14 Thread Erich Neuwirth
Thanks Duncan, this solved my problem. Here is another thingy I noticed \title{RExcel - Using \R from within Excel} produces RExcel - Using list() from within Excel So the \R macro cannot be used in titles. Is this intentional? Duncan Murdoch wrote: On 10/14/2009 7:45 AM, Erich Neuwirth

Re: [R] SPSS long variable names

2009-10-14 Thread Robert Baer
The problem is the limit of 8 characters long on variable names. And again, my answer is that one approach would be to map SHORT names to long variable LABELS. This was a common use of labels before variable names supported 64 bytes. After reading into R with read.spss() you and easily

Re: [R] currency conversion function?

2009-10-14 Thread David Winsemius
There was such a function in the form of getFX in package=quantmod , but testing it makes me think there might have been a change in how the website on which it depended might be expecting to get requests. library(quantmod) getFX(EUR/USD,from=as.Date(2008-01-01)) Error: oanda.com limits

Re: [R] ggplot2 scale_shape question

2009-10-14 Thread hadley wickham
Is there a way to have some points solid and some points hollow? I have two classes of points and there are so many points, that it's hard to see just the difference in shapes. I'd like to have one of the classes be hollow in addition to being a different shape. Any help would be grand.

Re: [R] plot discriminant analysis

2009-10-14 Thread Alain Guillet
Hi, I did it with Iris - data.frame(rbind(iris3[,,1], iris3[,,2], iris3[,,3]), Sp = rep(c(s,c,v), rep(50,3))) train - sample(1:150, 75) table(Iris$Sp[train]) z - lda(Sp ~ ., Iris, prior = c(1,1,1)/3, subset = train) Then I did plot(z,xlim=c(-10,10),ylim=c(-10,10)) before drawing

Re: [R] using mapply to avoid loops

2009-10-14 Thread David Winsemius
On Oct 14, 2009, at 12:58 AM, Stephen Samaha wrote: Hello, I would like to use mapply to avoid using a loop but for some reason, I can't seem to get it to work. I've included copies of my code below. The first set of code uses a loop (and it works fine), and the second set of code

Re: [R] currency conversion function?

2009-10-14 Thread Andreas Hüsler
Hi Liviu, try yahooSeries from fImport example: library(fImport) yahooSeries(EURUSD=X) Best, andreas -- www.er.ethz.ch/people/ahuesler __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] time grid for survfit Survival function outputs

2009-10-14 Thread Terry Therneau
... is it possible we could make survival function outputs on the pre-specified time grid with fixed increment and fixed length. Look at the help file for summary.survfit. Interpolating the raw data is somewhat harder than you might think for the number at risk component. fit - survfit(

[R] Handle lot of variables - Regression

2009-10-14 Thread anna0102
Hey, I've got a data set (e.g. named Data) which contains a lot of variables, for example: s1, s2, ..., s50 My first question is: It is possible to do this: Data$s1 But is it also possible to do something like this: Data$s1:s50 (I've tried a lot of versions of those without a result) My second

[R] problem when resizing graphics devices

2009-10-14 Thread Lucie Buchi
Dear R users, When I try to resize a graphics device in R, I often get this warning message (mostly when I already have several other graphics devices open) : Not enough memory to modify the size. Alpha blending is desactivated (translation from the french message: Mémoire insuffisante pour

[R] R version of MATLAB symbolic toolbox (variable substitution)

2009-10-14 Thread Jennifer Young
I'm translating some MATLAB code into R and have not found a simple equivalent of the function R = subs(S,old,new). I have, for example, a matrix such as this mx- function(){ matrix( c(0, f1, f2, s1, 0, 0, 0, s2, 0), 3,3, byrow=T) } and a matrix of data

Re: [R] problem when resizing graphics devices

2009-10-14 Thread Erik Iverson
Can you give us a reproducible example of R commands that cause these messages on your system? Like by creating simulated data.frames, opening several devices, etc? -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Lucie Buchi

Re: [R] Handle lot of variables - Regression

2009-10-14 Thread Dieter Menne
anna0102 wrote: I've got a data set (e.g. named Data) which contains a lot of variables, for example: s1, s2, ..., s50 My first question is: It is possible to do this: Data$s1 But is it also possible to do something like this: Data$s1:s50 (I've tried a lot of versions of those without

Re: [R] Strange characters that block import

2009-10-14 Thread Duncan Murdoch
On 10/14/2009 8:25 AM, arnaud Mosnier wrote: Dear useRs, I try to import a text file that contain some strange characters coming from the misinterpretation of foreign language characters by another software (see below). Here is an example of text with a

Re: [R] RCMD Rdconv debugging output

2009-10-14 Thread Duncan Murdoch
On 10/14/2009 8:42 AM, Erich Neuwirth wrote: Thanks Duncan, this solved my problem. Here is another thingy I noticed \title{RExcel - Using \R from within Excel} produces RExcel - Using list() from within Excel So the \R macro cannot be used in titles. Is this intentional? Yes, and it is

[R] Taking specific/timed differences in a zoo timeseries

2009-10-14 Thread Sergey Goriatchev
Hello everyone. I have a specific problem that I have difficulties to solve. Assume I have a zoo object: set.seed(12345) data - round(runif(27)*10+runif(27)*5, 0) dates - as.Date(c(09/03/09, 09/04/09, 09/07/09, 09/09/09, 09/10/09, 09/11/09, 09/14/09, 09/16/09, 09/17/09, 09/18/09, 09/21/09,

Re: [R] Linear Regression Question

2009-10-14 Thread Peter Ehlers
Alexandre, Let me add two small points to Ted's exposition: 1. you can use the extractor function coefficients(), or just coef() on the summary: coef(summary(regression)) which will also give you the matrix of estimates, etc. 2. You will find that using the function str() often is

[R] metaMDS NMDS: use of alternative distances?

2009-10-14 Thread Kim Vanselow
Dear r-helpers! How can I integrate other distances (in the form of a dist object) into function metaMDS? The problem: metaMDS needs the original data.frame for the calculation and only the default distances of function vegdist are allowed. Any suggestions are greatly appreciated! Thank you, Kim

Re: [R] Strange characters that block import

2009-10-14 Thread Prof Brian Ripley
On Wed, 14 Oct 2009, Duncan Murdoch wrote: On 10/14/2009 8:25 AM, arnaud Mosnier wrote: Dear useRs, I try to import a text file that contain some strange characters coming from the misinterpretation of foreign language characters by another software (see below).

Re: [R] Taking specific/timed differences in a zoo timeseries

2009-10-14 Thread Gabor Grothendieck
Try this: library(zoo) # temp - ... from post asking question # create a day sequence, dt, with no missing days # and create a 0 width series with those times. # merge that with original series giving original # series plus a bunch of times having NA values. # Use na.locf to fill in those values

Re: [R] metaMDS NMDS: use of alternative distances?

2009-10-14 Thread Gavin Simpson
On Wed, 2009-10-14 at 16:57 +0200, Kim Vanselow wrote: Dear r-helpers! How can I integrate other distances (in the form of a dist object) into function metaMDS? The problem: metaMDS needs the original data.frame for the calculation and only the default distances of function vegdist are

Re: [R] Selecting initial numerals

2009-10-14 Thread joris meys
Josh, One way would be to convert the numeric vector to a character and use the function substr(). Following code returns a numeric vector with the 2 first digits of every element. naics=c(238321, 624410, 484121 ,238911, 81, 531110, 621399,541613,524210 ,236115 ,811121 ,236115 ,236115

[R] axis label

2009-10-14 Thread maram salem
Hi all, I want the y-axis label to be ( in symbols) g(sigma given alpha) where given is the conditional sign. I've tried ylab=expression(g(sigma|alpha))) but it gave me g(|(sigma,alpha)) where the sigma and alpha are in greek but the conditional sign is misplaced (before the bracket) Any help

Re: [R] axis label

2009-10-14 Thread Jorge Ivan Velez
Hi Maram, How about this? plot(1, ylab = expression(sigma*|*alpha)) HTH, Jorge On Wed, Oct 14, 2009 at 11:12 AM, maram salem wrote: Hi all, I want the y-axis label to be ( in symbols) g(sigma given alpha) where given is the conditional sign. I've tried ylab=expression(g(sigma|alpha)))

Re: [R] axis label

2009-10-14 Thread Gabor Grothendieck
Try this: plot(0, main = ~ g(sigma * | * alpha)) On Wed, Oct 14, 2009 at 11:12 AM, maram salem marammagdysa...@yahoo.com wrote: Hi all, I want the y-axis label to be ( in symbols) g(sigma given alpha) where given is the conditional sign. I've tried ylab=expression(g(sigma|alpha))) but it

Re: [R] SPSS long variable names

2009-10-14 Thread joris meys
On Wed, Oct 14, 2009 at 2:45 PM, Robert Baer rb...@atsu.edu wrote: The problem is the limit of 8 characters long on variable names. And again, my answer is that one approach would be to map SHORT names to long variable LABELS.  This was a common use of labels before variable names supported

Re: [R] Creating a list of empty lists

2009-10-14 Thread Magnus Torfason
On 10/13/2009 10:06 AM, Henrique Dallazuanna wrote: Try this: replicate(3, list()) Thanks! I now have three ways to achieve my goal: 1: rep(list(list()), 3) 2: replicate(3, list()) 3: Due to the way R recycles arguments, I found that it is enough to have construct a list(list()), and

Re: [R] post-hoc test with kruskal.test()

2009-10-14 Thread Thomas Lumley
On Wed, 14 Oct 2009, Meyners, Michael, LAUSANNE, AppliedMathematics wrote: Robert, What do you mean by not symmetric? If you mean unbalanced in terms of sample size, that's not a problem if you choose the right specifications for wilcox.test. The Kruskal-Wallis-Test is a generalization of

Re: [R] SPSS long variable names

2009-10-14 Thread Orvalho Augusto
Hi The .dat file is a tab delimited file with the long variables names on it. The .sps file has the instructions to read the .dat and place all the variable and value labels. The ideia of reading the dat directely is good but I need the labels placed. Yes I could read the dat file and parese

Re: [R] Time Dependent Cox Model

2009-10-14 Thread quaildoc
Does anyone have suggestions? Thanks! quaildoc wrote: I am having trouble formatting some survival data to use in a time dependent cox model. My time dep. variable is habitat and I have it recorded for every day (with some NAs). I think it is working properly except for calculating the

Re: [R] Taking specific/timed differences in a zoo timeseries

2009-10-14 Thread Sergey Goriatchev
Dear Gabor, Thank you very much for your help! I'm now using your suggestion with my data. May I ask a stupid question? The output's index now has format 2009-10-14. How can I transform it back into original 10/14/09 and use this in a zoo object? Regards, Sergey On Wed, Oct 14, 2009 at 17:03,

[R] Problem with NLSstClosestX; and suggested fix

2009-10-14 Thread Keith Jewell
Problem is demonstrated with this code, intended to find the approximate 'x' at which the 'y' is midway between the left and right asymptotes. This particular data set returns NA, which is a bit silly! -- sXY - structure(list(x = c(0, 24, 27, 48, 51, 72, 75, 96, 99), y = c(4.98227,

[R] Getting indeices of intersecting elements.

2009-10-14 Thread Praveen Surendran
Hi, Is there a command to get the indices of intersecting elements of two vectors as intersect() will give the elements and not its indices. Thanks in advance. Praveen Surendran School of Medicine and Medical Sciences University College Dublin Belfield, Dublin 4 Ireland.

Re: [R] plot discriminant analysis

2009-10-14 Thread Alejo C.S.
Hi Alain, thanks for the fast response. I've the same results with iris data, but when I use my data (mentioned in the first message), I have different results. Regards, Alejo 2009/10/14 Alain Guillet alain.guil...@uclouvain.be Hi, I did it with Iris - data.frame(rbind(iris3[,,1],

Re: [R] Time Dependent Cox Model

2009-10-14 Thread joris meys
Well, it might be wise to elaborate a bit more about the variables and what exactly you want e.g. death-time to be. I'd interprete it as time of death, but the fact that it is 0/1, means it is a logical (?) binary variable of some sort. Please ask your question in such a way that somebody who

[R] Survival and nonparametric

2009-10-14 Thread Ashta
Hi all, Has any body the exprience to iclude a nonparametric component into the survival analysis using R package? *Can someone recommend *me * some ** references? * Thanks a lot Ashta [[alternative HTML version deleted]] __

Re: [R] plot discriminant analysis

2009-10-14 Thread David Winsemius
On Oct 14, 2009, at 12:24 PM, Alejo C.S. wrote: Hi Alain, thanks for the fast response. I've the same results with iris data, but when I use my data (mentioned in the first message), You are apparently under the false impression that the data made it through the listserv. Read the

Re: [R] Getting indeices of intersecting elements.

2009-10-14 Thread David Winsemius
On Oct 14, 2009, at 12:15 PM, Praveen Surendran wrote: Hi, Is there a command to get the indices of intersecting elements of two vectors as intersect() will give the elements and not its indices. ?which samp1 - sample(seq(3,198, by=3), 20); samp2 - sample(seq(3,198, by=3), 20) int

Re: [R] currency conversion function?

2009-10-14 Thread Liviu Andronic
Hello On 10/14/09, Henrique Dallazuanna www...@gmail.com wrote: foo('BRL', 'USD', '2009-10-14') Nice function, thank you. Two issues, though: - it seems to provide reverse output. Example: ## how many dollars do you get from one euro? foo('EUR', 'USD', '2009-10-14') [1] 0.67544 ## however,

Re: [R] Survival and nonparametric

2009-10-14 Thread David Winsemius
http://finzi.psych.upenn.edu/views/Survival.html Kalbfleisch Prentice The Statistical Analysis of Failure Time Data Therneau Grambsch Modeling Survival Data Harrell Regression Modeling Strategies On Oct 14, 2009, at 12:35 PM, Ashta wrote: Hi all, Has any body the exprience to iclude a

Re: [R] SPSS long variable names

2009-10-14 Thread joris meys
You could read in the sps file in R with readLines() for example, and then use the tools for regular expression and substring to find the label statements. Then you can just use R to add the labels to it, without having to pass through PSPP. So you could actually just create an R script that

Re: [R] Survival and nonparametric

2009-10-14 Thread David Winsemius
On Oct 14, 2009, at 12:47 PM, David Winsemius wrote: http://finzi.psych.upenn.edu/views/Survival.html # lots of 404 link errors The finzi.psych server does something weird to the task view pages, so you would get greater linkability with the actual CRAN version:

Re: [R] Tunnelling X for R graphics

2009-10-14 Thread Santosh
In an ssh client, under connections section, there is an option for tunneling. Please ensure that the tunneling options are turned on.. and if applicable, incoming/outgoing tunnels, listen/destination ports, etc. are set. Thanks, Santosh On Tue, Feb 3, 2009 at 2:03 PM, Adam D. I. Kramer

[R] Error: testing 'stats' failed - R 2.9.2 on Linux

2009-10-14 Thread John Murdie
I've just built R 2.9.2 from source on Slackware Linux 13.0 - 32-bit (will try 64-bit also next) - and seen: Collecting examples for package 'stats' Running examples in package 'stats' Error: testing 'stats' failed Execution halted make[3]: *** [test-Examples-Base] Error 1 Looking at

Re: [R] pairs

2009-10-14 Thread Greg Snow
Does the pairs2 function in the TeachingDemos package do what you want? -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- project.org] On

Re: [R] SPSS long variable names

2009-10-14 Thread Orvalho Augusto
Ok. I will try that then. Caveman On Wed, Oct 14, 2009 at 6:53 PM, joris meys jorism...@gmail.com wrote: You could read in the sps file in R with readLines() for example, and then use the tools for regular expression and substring to find the label statements. Then you can just use R to add

[R] Large data sets with high dimensional fixed effects

2009-10-14 Thread dcervone
Hi, I have a data set that consists of about 2 million observations and several high dimensional fixed effects (2 factors at around 1000 levels each, and others with a few hundred levels). I'm looking to run linear and logit regressions. I've tried packages such as filehash and biglm to store

Re: [R] Getting indeices of intersecting elements.

2009-10-14 Thread William Dunlap
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Praveen Surendran Sent: Wednesday, October 14, 2009 9:15 AM To: r-help@r-project.org Subject: [R] Getting indeices of intersecting elements. Hi, Is there a command to get

[R] RPostgreSQL: unable to load shared library

2009-10-14 Thread Fanfaar
Hello list, I'm using R 2.9.2 on a WinXP system, and I installed the RPostgreSQL library using the package installer. When trying to load it, I get the following error: library('RPostgreSQL') Error in inDL(x, as.logical(local), as.logical(now), ...) : unable to load shared library

[R] puzzle using gsub (and encodings maybe)

2009-10-14 Thread Adrian Dragulescu
Hello, Below is some output that shows my issue. I have a variable x that I read from a file (more on this below) x [1] NEW YORK NEW ENGLAND gsub( -, -, x)# this does not work! [1] NEW YORK NEW ENGLAND Encoding(x) # is x in a special encoding? no [1]

Re: [R] puzzle using gsub (and encodings maybe)

2009-10-14 Thread Duncan Murdoch
On 10/14/2009 1:30 PM, Adrian Dragulescu wrote: Hello, Below is some output that shows my issue. I have a variable x that I read from a file (more on this below) x [1] NEW YORK NEW ENGLAND gsub( -, -, x)# this does not work! [1] NEW YORK NEW ENGLAND It looks as though it

Re: [R] puzzle using gsub (and encodings maybe)

2009-10-14 Thread Adrian Dragulescu
charToRaw(x) [1] 4e 45 57 20 59 4f 52 4b 20 ad 4e 45 57 20 45 4e 47 4c 41 4e 44 charToRaw(y) [1] 4e 45 57 20 59 4f 52 4b 20 2d 4e 45 57 20 45 4e 47 4c 41 4e 44 So they are different. Adrian I use R 2.8.1 on WinXP On Wed, 14 Oct 2009, Duncan Murdoch wrote: On 10/14/2009 1:30 PM,

Re: [R] puzzle using gsub (and encodings maybe)

2009-10-14 Thread Duncan Murdoch
On 10/14/2009 1:41 PM, Adrian Dragulescu wrote: charToRaw(x) [1] 4e 45 57 20 59 4f 52 4b 20 ad 4e 45 57 20 45 4e 47 4c 41 4e 44 charToRaw(y) [1] 4e 45 57 20 59 4f 52 4b 20 2d 4e 45 57 20 45 4e 47 4c 41 4e 44 So they are different. Adrian I use R 2.8.1 on WinXP But that's

[R] Scatter plot using icons (from a gif) instaed of points - is it possible ?

2009-10-14 Thread Tal Galili
Hello dear R-help group. I wish to plot a scatter plot using icons (or images) instead of points. Is it possible? and how so? Thanks, Tal -- My contact information: Tal Galili E-mail: tal.gal...@gmail.com Phone number: 972-52-7275845 FaceBook:

Re: [R] Handle lot of variables - Regression

2009-10-14 Thread Frank E Harrell Jr
anna0102 wrote: Hey, I've got a data set (e.g. named Data) which contains a lot of variables, for example: s1, s2, ..., s50 My first question is: It is possible to do this: Data$s1 But is it also possible to do something like this: Data$s1:s50 (I've tried a lot of versions of those without a

Re: [R] puzzle using gsub (and encodings maybe)

2009-10-14 Thread Prof Brian Ripley
On Wed, 14 Oct 2009, Adrian Dragulescu wrote: charToRaw(x) [1] 4e 45 57 20 59 4f 52 4b 20 ad 4e 45 57 20 45 4e 47 4c 41 4e 44 charToRaw(y) [1] 4e 45 57 20 59 4f 52 4b 20 2d 4e 45 57 20 45 4e 47 4c 41 4e 44 So they are different. We really do need the 'at a minimum' information we asked

Re: [R] puzzle using gsub (and encodings maybe)

2009-10-14 Thread Adrian Dragulescu
I get the same results (not working) using R 2.9.2 and R.10.0 beta. Thank you for looking at this. On Wed, 14 Oct 2009, Duncan Murdoch wrote: On 10/14/2009 1:41 PM, Adrian Dragulescu wrote: charToRaw(x) [1] 4e 45 57 20 59 4f 52 4b 20 ad 4e 45 57 20 45 4e 47 4c 41 4e 44 charToRaw(y)

Re: [R] puzzle using gsub (and encodings maybe)

2009-10-14 Thread Duncan Murdoch
On 10/14/2009 2:16 PM, Adrian Dragulescu wrote: I get the same results (not working) using R 2.9.2 and R.10.0 beta. But it is working: the dash is an ad in x, not a 2d. You need to ask to substitute for the ad character, e.g. by spacelongdash - rawToChar(as.raw(c(0x20, 0xad)))

[R] Cacheing computationally expensive getter methods for S4 objects

2009-10-14 Thread Steve Lianoglou
Hi, I was wondering if there was a way to store the results of a computationally expensive getter call on an S4 object, so that it is only calculated once for each object. Trivial example: let's say I want to cache the expensive area calculation of a square object. setClass(Square,

Re: [R] puzzle using gsub (and encodings maybe)

2009-10-14 Thread Adrian Dragulescu
Thank you. If I use gsub( \xad, -, x) [1] NEW YORK-NEW ENGLAND I get what I want. Adrian sessionInfo() R version 2.9.2 (2009-08-24) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United

Re: [R] RPostgreSQL: unable to load shared library

2009-10-14 Thread Prof Brian Ripley
On Wed, 14 Oct 2009, Fanfaar wrote: Hello list, I'm using R 2.9.2 on a WinXP system, and I installed the RPostgreSQL library using the package installer. When trying to load it, I get the following error: library('RPostgreSQL') Error in inDL(x, as.logical(local), as.logical(now), ...) :

  1   2   >