[R] Scale set of 0 values returns NAN??

Hi, More questions in my ongoing quest to convert from RapidMiner to R. One thing has become VERY CLEAR: None of the issues I'm asking about here are addressed in RapidMiner. How it handles misisng values, scaling, etc. is hidden within the black box. Using R is forcing me to take a much

Re: [R] Strange column shifting with read.table

Hi, Thanks for the continued support. I've been working on this all night, and have learned some things: 1) Since I'm really committed to using an SVM, I need to skip the examples with missing data. I have a training set of approximately 22,000 examples of which about 500 have missing

Re: [R] i'm so stuck with text file and contour plot

Hi David, I think he may also need to add the header=TRUE argument: No! The argument header= is not required in this case. ## TDat - read.csv(small.txt, sep=\t) str(TDat[,1:3]) 'data.frame': 10 obs. of 3 variables: $ Placename: Factor w/ 10 levels Aankoms,Aapieshoek,..: 1 2 3 4 5 6 7 8

Re: [R] i'm so stuck with text file and contour plot

And I meant to add, but somehow forgot, that the default for read.csv is header=TRUE (which is different from read.table, where it is FALSE). Regards, Mark. Mark Difford wrote: Hi David, I think he may also need to add the header=TRUE argument: No! The argument header= is not required

[R] [R-pkgs] New package for interactive brushing/cropping/deleting points from RGL : sculpt3d v 0.2

sculpt3d is a GTK+ toolbar that allows for more interactive control of a dataset inside the RGL plot window. Controls for simple brushing, highlighting, labeling, and mouseMode changes are provided by point-and-click rather than through the R terminal interface. More discussion and demo movie:

Re: [R] How to stop an R script when running JGR on a Linux/SuSE system

PD == Peter Dalgaard p.dalga...@biostat.ku.dk on Sat, 01 Aug 2009 14:08:30 +0200 writes: PD Bernd Kreuss wrote: sorry for the eventual double posting, but i got a strange error from a versatel(???) server about not enough quota when replying to the message PD Yes, I

Re: [R] .... running JGR on a Linux/SuSE system

BK == Bernd Kreuss prof7...@googlemail.com on Sat, 01 Aug 2009 12:40:45 +0200 writes: BK mau...@alice.it wrote: I wonder whether there is a more gentle way to stop an R script running on top of JGR aother than ... unplugging the power cord. BK there must be a bug in

Re: [R] Scale set of 0 values returns NAN??

Dear Noah, You get Nan values because you devide by zero is all values are equal. So it does not only happen when you have only zero's. Try scale(rep(1, 10)) A workaround is to turn the scaling off when all values are identical (have zero variance). You can rearrange your nested loops with the

Re: [R] What does .[foo] really mean?

It's a place holder. See http://en.wikipedia.org/wiki/Foo HTH KJ Simon Tian simon...@gmail.com wrote in message news:2119d1440908012150j373bfaf5v9d498bded94f...@mail.gmail.com... Hi R users, I really want to know what exactly .[foo] means. Thanks in advance. -Simon [[alternative HTML

Re: [R] .... running JGR on a Linux/SuSE system

Martin Maechler wrote: But there are quite a few others (I do not really know from personal experience), see, e.g., all the projects (listed on the left hand side) of http://www.r-project.org/GUI I personally can recommend rkward which is now available on many platforms and which is

[R] R GUI page on the wiki (was: Re: .... running JGR on a Linux/SuSE system)

Hi, I've started the process of moving the page R gui page [1] to the R wiki [2]. The idea is that instead of being Philippe Grosjean's burden to update the static page, it becomes a collaborative exercize. So that it can be more up to date and can also contain things as use case, feature

Re: [R] combine venn diagram and pie chart

kevinchang wrote: Hi R users, I am wondering if it there is R's function that can help integrate venn diagram and pie chart to compare two related datasets. I know the package limma has something built-in for making venn diagram, but I guess it would be very painful to use line and text to

Re: [R] Problem with RGtk2 Rattle

Hi I am pleased to report that I used a windows solution (deleted everything then re-installed it all again) this time there are no problems - thanks once again Regards Wayne Graham Williams wrote: 2009/8/1 Felix Andrews fe...@nfrac.org Um, it sounds as if you are trying to install

[R] [R-pkgs] new package: perm

Hi all, I have uploaded a new package called perm which does linear permutation tests. It does two-sample, k-sample and trend tests. The package is redundant, since everything it does can be done in the coin package, and it is usually faster in coin. The main motivation is to provide another

[R] [R-pkgs] new package: interval

I have just uploaded a new package called interval. It does weighted logrank tests (e.g., logrank, Wilcoxon-type tests) on interval censored data. It allows Surv objects and can calculated exact p-values. It calls the perm package, but also allows calculation of scores that may be output to

[R] apply cut function on various values (factor constrained)

hi all, presumably i have a small request, but i do indeed struggle. i want to use the cut function as follows. 1) imagine two vectors X1 and X2 2) code values of X1 due to it's quantile values (as this is ordinary use of the cut function) 3) code values of X2 due to it's quantile values, by

Re: [R] xyplot: superpose 2 time series with different time intervals

Thank you Deepayan. Your suggestion worked perfectly. Thanks too to Gabor for your suggestion and for your help trying to get this help request through to r-help before I subscribed. I'm all set now. Gary On Sat, Aug 1, 2009 at 1:17 PM, Deepayan Sarkardeepayan.sar...@gmail.com wrote: On Sat,

[R] How to get w and b in SVR? (package e1071)

Dear R users, I'm running a SVR in package e1071 but I did not able to calculate the parameters w and b of the regression. I don't know how to do that and if it is possible to do it with this package. Someone have some idea. Any help would be much appreciated. Marlene [[alternative

[R] apply function to named numeric object

Maybe a very basic question. How can I get mean of a object like this str(object) num [1:53, 1:2] 734 734 736 734 736 ... - attr(*, dimnames)=List of 2 ..$ : NULL ..$ : chr [1:2] x y I want the mean of both x and y ___ Paulo E. Cardoso [[alternative HTML version

Re: [R] diagonal LDA and QDA

cG == cindy Guo cindy.g...@gmail.com on Sat, 1 Aug 2009 13:03:08 -0700 writes: cG Hi, all, I am wondering if there is any package doing cG lda and qda which allows assuming diagonal covariance cG matrices. I checked the lda function in MASS, and it cG seems it does not

Re: [R] apply function to named numeric object

On Aug 3, 2009, at 9:07 AM, Paulo Eduardo Cardoso wrote: Maybe a very basic question. How can I get mean of a object like this str(object) num [1:53, 1:2] 734 734 736 734 736 ... - attr(*, dimnames)=List of 2 ..$ : NULL ..$ : chr [1:2] x y I want the mean of both x and y apply(object,

Re: [R] diagonal LDA and QDA

On Aug 1, 2009, at 4:03 PM, cindy Guo wrote: Hi, all, I am wondering if there is any package doing lda and qda which allows assuming diagonal covariance matrices. I checked the lda function in MASS, and it seems it does not support this. Did you try one of the r-search sites? e.g.:

[R] apply cut function on various values (factor constrained)

hi all, presumably i have a small request, but i do indeed struggle. i want to use the cut function as follows. 1) imagine two vectors X1 and X2 2) code values of X1 due to it's quantile values (as this is ordinary use of the cut function) 3) code values of X2 due to it's quantile values, by

Re: [R] apply function to named numeric object

Thanks. Works A doubt. I asked that because I have a loop that can generate objects like that or this one str(object) Named num [1:2] 670 314 - attr(*, names)= chr [1:2] x y To avoid using a ifelse() to calculate means of x and y, I'd like a function that handle with this named numeric object

[R] calculate glm slope

Hello list-members! Is there a function that automatically calculates the slope and the confidence intervals for the different explanatory variables from the glm output? Many thanks in advance for every hint. B. = Birgit Lemcke Institut of Systematic Botany

Re: [R] How to get w and b in SVR? (package e1071)

Hi Marlene, On Aug 3, 2009, at 8:59 AM, marlene marchena wrote: Dear R users, I'm running a SVR in package e1071 but I did not able to calculate the parameters w and b of the regression. I don't know how to do that and if it is possible to do it with this package. Someone have some idea.

Re: [R] apply function to named numeric object

On Aug 3, 2009, at 9:52 AM, Paulo Eduardo Cardoso wrote: Thanks. Works A doubt. I asked that because I have a loop that can generate objects like that or this one str(object) Named num [1:2] 670 314 - attr(*, names)= chr [1:2] x y To avoid using a ifelse() to calculate means of x

Re: [R] calculate glm slope

On Aug 3, 2009, at 9:58 AM, Birgit Lemcke wrote: Hello list-members! Is there a function that automatically calculates the slope and the confidence intervals for the different explanatory variables from the glm output? Many thanks in advance for every hint. B.; meet ?? ,the quick way

Re: [R] apply function to named numeric object

Both objects are numeric, but they differ in structure and I don't know how to handle both with a single function to get means. When I try to apply colMeans to a named numeric object like this: Named num [1:2] 670 314 - attr(*, names)= chr [1:2] x y I get the error messages: Error in

Re: [R] Cox ridge regression

Question 1. Consider the following example from help(ridge): fit1 - coxph(Surv(futime, fustat) ~ rx + ridge(age, ecog.ps, theta=1), ovarian) As I understand, this builds a model in which `rx' is the predictor, whereas ridge penalty term contains variables `age' and `ph.ecog'. Could someone

Re: [R] calculate glm slope

On Aug 3, 2009, at 10:13 AM, David Winsemius wrote: On Aug 3, 2009, at 9:58 AM, Birgit Lemcke wrote: Hello list-members! Is there a function that automatically calculates the slope and the confidence intervals for the different explanatory variables from the glm output? Many thanks

Re: [R] apply function to named numeric object

So you need to show us how you made the second object. (The Posting Guide's phrase reproducible code is hereby invoked.) On Aug 3, 2009, at 10:18 AM, Paulo Eduardo Cardoso wrote: Both objects are numeric, but they differ in structure and I don't know how to handle both with a single

Re: [R] apply function to named numeric object

That may just be a vector with two elements, so functions designed for matrices are not going to work. Why are you taking means of such small data objects anyway? obj - c(x=640, y=410) obj x y 640 410 str(obj) Named num [1:2] 640 410 - attr(*, names)= chr [1:2] x y mean(obj)

Re: [R] apply function to named numeric object

The objects are centroid x y 669.5 313.5 str(centroid) Named num [1:2] 670 314 - attr(*, names)= chr [1:2] x y or centroid x y [1,] 659.5 472.5 [2,] 659.5 473.5 str(centroid) num [1:2, 1:2] 660 660 472 474 - attr(*, dimnames)=List of 2 ..$ : NULL ..$ : chr [1:2]

Re: [R] apply function to named numeric object

Dear Paulo, Your object is a named vector. Apply() and colMeans() only work on matrices and dataframes. So you need to convert the vector to a matrix or data.frame. See the example below. Please not that calculating the mean of one element is a waiste of time. # a named numerical vector object -

[R] Table to Data-frame

Hello How to pass from table to data-frame ? Thanks __ 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 and provide commented, minimal,

Re: [R] Table to Data-frame

On Aug 3, 2009, at 10:35 AM, Olivier Peron wrote: Hello How to pass from table to data-frame ? Copied from the help page for table: The as.data.frame method for objects inheriting from class table can be used to convert the array-based representation of a contingency table to a data

Re: [R] apply function to named numeric object

Yes, mean of a one element object is meanless, but I'm trying to avoid a ifelse in the sentence. The result of mean values of x and y is added at each loop to a data.frame (or any other object that could speed up the loop. 2009/8/3 ONKELINX, Thierry thierry.onkel...@inbo.be Dear Paulo, Your

Re: [R] Table to Data-frame

as.data.frame(table())? 2009/8/3 Olivier Peron olivier.pe...@univ-pau.fr Hello How to pass from table to data-frame ? Thanks __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] Table to Data-frame

Hi, On Aug 3, 2009, at 10:35 AM, Olivier Peron wrote: Hello How to pass from table to data-frame ? R a - sample(LETTERS[1:10], 20, replace=T) R b - table(a) R c - as.data.frame(b) R c a Freq 1 A4 2 B2 3 C1 4 D4 5 E1 6 F1 7 H1 8 J6 -steve -- Steve

Re: [R] Table to Data-frame

Hi Olivier, How about this? # some data x - rpois(100,5) # the table table(x) # converting the table as.data.frame.table(table(x)) HTH, Jorge On Mon, Aug 3, 2009 at 10:35 AM, Olivier Peron olivier.pe...@univ-pau.frwrote: Hello How to pass from table to data-frame ? Thanks

Re: [R] min frequencies of categorical predictor variables in GLM

On Aug 3, 2009, at 12:06 AM, Thomas Mang wrote: Hi, Suppose a binomial GLM with both continuous as well as categorical predictors (sometimes referred to as GLM-ANCOVA, if I remember correctly). For the categorical predictors = indicator variables, is then there a suggested minimum

[R] selectively altering variable value

Hello, I have run an eye-tracking experiment for which I now like to analyse the saccades. Participants looked from a fixation cross (ia = 5) to the target area (ia = 4) in following example of a data frame. ia = 9 stands for everything else. A saccade is indicated by saccade = 1. Sometimes the

[R] problem with predict(...,type=terms)

Hello R users, I'm trying to understand the meaning of the results returned by the command 'predict' with type=terms. I thought it was the matrix with the fitted values for each term (value of predictor * related coef), but the calculation is different. I include an example below: y -

[R] plm summary error

Dear plm-Package insiders, [I posted the following observation is April already but unfortunately I am not aware of any answers. With the hope that someone found an answer in the mean time, I ask again:] I realized the following difficulty with the summary.plm function (demonstrated with

[R] Variable names as inputs...

Hi, I have been using R for a bit, but never came across this simple issue... So, I wonder if anyone could give me a hint on it. I have calculated a bunch o models and now need to compare them using AIC(). The issue seem to be that I have placed each output in a list object. So, the question is

[R] penalized logistic regression

Hi, R users, Is there any package for penalized logistic regression with more than two response classes? I read the manual for stepPlr, but it seems it's only for binary case. Thank you, Annie [[alternative HTML version deleted]] __

[R] How to catch an error using try

Sometimes the following function call causes a database exception: gene.seq - getSequence (id=gene.map[,ensembl_transcript_id], type=ensembl_transcript_id, + seqType=3utr, mart=hmart) I understand the above function must be called by try to capture the eventual

[R] survdiff for left-truncated data?

Hi Does anyone know if there is a function like survdiff which can also handle left-truncated and right-censored data? When I use it on left-truncated and right-censored data I get an error message saying Right censored data only. Many thanks Rajen [[alternative HTML version deleted]]

Re: [R] plm summary error

Ulrich, as usual with such specific questions, it is a good idea to have the package authors at least Cc (done now). In April I had already forwarded your message privately to the plm authors with the following comment: comment I've quickly checked what is going on. If you do something like

Re: [R] selectively altering variable value

This should set the first two ia==4 entries to 1 in saccade column: x - read.table('/tempxx.txt.r', header=TRUE) x vp trial ia saccade 1 1 1 5 0 2 1 1 5 0 3 1 1 5 0 4 1 1 5 1 5 1 1 5 1 6 1 1 9 1 7 1 1 9

Re: [R] How to catch an error using try

Hi, On Aug 3, 2009, at 11:53 AM, mau...@alice.it wrote: Sometimes the following function call causes a database exception: gene.seq - getSequence (id=gene.map[,ensembl_transcript_id], type=ensembl_transcript_id, + seqType=3utr, mart=hmart) I understand the

Re: [R] How to catch an error using try

gotError - try(gene.seq - getSequence (id=gene.map[,ensembl_transcript_id], type=ensembl_transcript_id, seqType=3utr, mart=hmart)) if (inherits(gotError, 'try-error')){.error occurred, do recovery} On Mon, Aug 3, 2009 at 11:53 AM, mau...@alice.it wrote: Sometimes the following

[R] For Loop Iteration: Variable Counter

Hey guys, How do I iterate such that I add 100 to the counter every time? Suppose: for (i in c(1:100)) I want i to be 1, 10, 20, 30, ... instead of 1,2,3,4,5 ... How can this be done? Thanks, Vivek __ R-help@r-project.org mailing list

Re: [R] For Loop Iteration: Variable Counter

Hi, On Aug 3, 2009, at 12:09 PM, Vivek Ayer wrote: Hey guys, How do I iterate such that I add 100 to the counter every time? Suppose: for (i in c(1:100)) I want i to be 1, 10, 20, 30, ... instead of 1,2,3,4,5 ... How can this be done? ?seq -steve -- Steve Lianoglou Graduate Student:

Re: [R] For Loop Iteration: Variable Counter

Got it..Thanks Vivek On Mon, Aug 3, 2009 at 9:13 AM, Steve Lianogloumailinglist.honey...@gmail.com wrote: Hi, On Aug 3, 2009, at 12:09 PM, Vivek Ayer wrote: Hey guys, How do I iterate such that I add 100 to the counter every time? Suppose: for (i in c(1:100)) I want i to be 1, 10, 20,

Re: [R] Cox ridge regression

Thank you Terry, that answered all questions. As a suggestion, help page for ridge() might indicate that the ridge term simultaneously introduces predictors and penalizes them. Ljubomir From: Terry Therneau thern...@mayo.edu To: Ljubomir Buturovic ljubo...@sfsu.edu Cc: r-help@r-project.org

[R] principal component analysis for class variables

Dear Forum, I have a class variable 1 (populations A-E), and two other class variables, variable 2 and variable 3. What I want is to see if the combination of var 2 and var 3, will give me a pattern that allows to distinguish populations. I found several packages like ade4, with pcaiv function

[R] wireframe color

Hello, hopefully simple question, but I cannot find the answer. I need to change the color from the standard default. Still want the scaled colors, but need different colors for different graphs. Code is: wireframe(z ~ y*x, mat.df, drape = TRUE, zlab = list(Water mass error

Re: [R] principal component analysis for class variables

The PCA doesn't work with class variables so the error is normal. You should try to work with a discriminant factorial analysis (see discrimin.coa in ade4). Alain andreiabb wrote: Dear Forum, I have a class variable 1 (populations A-E), and two other class variables, variable 2 and

Re: [R] principal component analysis for class variables

Hi, R users, I'm using the lme function in R to estimate a 2 level mixed effects model, in which the size of the groups are different. It turned out that It takes forever for R to converge. I also tried the same thing in SPSS and SPSS can give the results out within 20 minutes. Anyone can give

[R] lme funcion in R

Hi, R users, I'm using the lme function in R to estimate a 2 level mixed effects model, in which the size of the subject groups are different. It turned out that It takes forever for R to converge. I also tried the same thing in SPSS and SPSS can give the results out within 20 minutes. Anyone

Re: [R] lme funcion in R

On 8/3/2009 1:15 PM, Hongwei Dong wrote: Hi, R users, I'm using the lme function in R to estimate a 2 level mixed effects model, in which the size of the subject groups are different. It turned out that It takes forever for R to converge. I also tried the same thing in SPSS and SPSS can

Re: [R] Lattice xyplot: same scales within one factor

On Mon, Jul 13, 2009 at 8:33 AM, OBo.borderl...@gmail.com wrote: I am using R 2.8.1 and lattice to produce xyplots conditioned on two factors. What I would like is to have the scales be free between values of one factor, but some within. Thus, in this example, xyplot(mpg ~ disp | factor(gear)

Re: [R] lme funcion in R

On 2009.08.03 10:15:46, Hongwei Dong wrote: Hi, R users, I'm using the lme function in R to estimate a 2 level mixed effects model, in which the size of the subject groups are different. It turned out that It takes forever for R to converge. I also tried the same thing in SPSS and SPSS can

Re: [R] lme funcion in R

Thanks for the replies above. Here are my script and data structure: library(nlme) tlevel-lme(fixed = LN_unitlandval ~ MH_D+APT_D+ResOth_D+NonRes_D+Vacant_D+access_emp1+pct_vacant+transit_D+park_dum,data=lusdrdata,random = ~MH_D+APT_D+ResOth_D+NonRes_D+Vacant_D | TAZ) str: $ TAZ : int 100 100

Re: [R] principal component analysis for class variables:somenone has experience with made 4 and plotgenes?

Alain Guillet-2 wrote: The PCA doesn't work with class variables so the error is normal. You should try to work with a discriminant factorial analysis (see discrimin.coa in ade4). Alain thanks for the tip but it did not work, well then I manaege to put it to work but the result is

Re: [R] dot plot with several points for 2 categories

On Wed, Jul 29, 2009 at 9:11 AM, jaregisuck...@mpi-cbg.de wrote: Hi Michael, Steve, and 1Rnwb, I'm very impressed by the quick replies on the mailer. Thanks a lot for your suggestions. They worked very well. In general, I have to say that I'm a bit disappointed that in R, like in Excel,

Re: [R] boxplot( ) headers with Greek letters, values, and text

Johnson, Roger W. wrote: Any suggestions on how I can produce the same title using the boxplot command? If you use LaTeX you could try the tikzDevice package ( https://r-forge.r-project.org/projects/tikzdevice/ ). In that case you would use: main = paste('\\mu =',mu,', \\lambda

Re: [R] lattice shingle plot axis annotation

On Thu, Jul 30, 2009 at 12:43 AM, Armin Goralczykagoralc...@gmail.com wrote: Hello (R-)Experts I hope someone can help with this problem concerning axis annotation of a lattice shingle plot. I want a plot with three shingles to display some laboratory value over time. In the first panel over

Re: [R] Defining different objects within a loop

On Thu, Jul 30, 2009 at 3:09 PM, jimdarejamesdar...@gmail.com wrote: Dear R Users, I have a dataset that I have split by group.  For each group I plot a figure using: for (i in splitdf){ plot-xyplot() print(plot) savePlot(filename=paste(i$Group[1]),type=pdf) } This

[R] Open data-file in R - from Rserve

How can I open a data file from the client on the RServe server using Java? -- View this message in context: http://www.nabble.com/Open-data-file-in-R---from-Rserve-tp24795557p24795557.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] wireframe color

On Mon, Aug 3, 2009 at 9:36 AM, Marlin Keith Coxmarlink...@gmail.com wrote: Hello, hopefully simple question, but I cannot find the answer.  I need to change the color from the standard default.  Still want the scaled colors, but need different colors for different graphs. Code is:

Re: [R] Lattice xyplot: same scales within one factor

Thank you! I ended up hunting down a few previous posts you made and getting it to to work, although yes, it was not pretty. -Orion On 8/3/09 10:30 am, Deepayan Sarkar deepayan.sar...@gmail.com wrote: On Mon, Jul 13, 2009 at 8:33 AM, OBo.borderl...@gmail.com wrote: I am using R 2.8.1 and

Re: [R] boxplot( ) headers with Greek letters, values, and text

Try this: boxplot(1:5, main = bquote(mu == .(mu) * , ~ lambda == .(lambda) * , ~ trunction == .(truncation) * , ~ N[T] == .(n))) On Sun, Aug 2, 2009 at 8:37 PM, Johnson, Roger W.roger.john...@sdsmt.edu wrote: Hi - I've been using the option main=bquote(paste(mu==.(mu),,

[R] session logging

Consider all the text that one sees on the console during an R session. Is there a way, within R, to make all this text--both the output and the messages--automatically get copied to a single text file, in addition to seeing it on the console? If I remember to save the console to a file at

[R] MSM package and qmatrix

Hi R-helpers, I am having a problem understanding how to construct the qmatrix, i have read the help menu for msm but I am still a bit lost. Below is an example of one of my transition matrices: statetable.msm(BEH, ID, data = A1) to from 1 2 3 4 5 23 200 201 203 999 1 86 11 2

Re: [R] status of boolean?

Ronggui Huang wrote: package of boolean has been removed from CRAN. Anyone know why? Thanks. Yes, because the maintainer has not fixed the issues the package had with R CMD check and got ORPHANED. You can try to contact the maontainer in order to take over maintainer ship. Unfortunately,

[R] session logging

Consider all the text that one sees on the console during an R session. Is there a way, within R, that all this text--i.e., both the output and the messages--can be automatically copied to a single file, in addition to its being visible on the console? If I remember to save the console to a file

[R] Installing package for i386

I am trying to install a package from a collaborator that was given to me as a zip file. I'm running R 2.9.0 using Mac OSX 10.5.7 [R.app GUI 1.28 (5395) i386-apple-darwin8.11.1] foo.zip contains all the usual package directories (man, R, R-ex, etc) After unzipping foo, I used the following

[R] Help with Ecdf function

Dear R users, I'm using Ecdf (Hmisc library) to plot four cdf in a same graphic. In this graphic I also plot the 0.99 quantile for these cdf. I successfully plot cdfs using different types of line to distinguish them, but I can't determine the type of lines showing 0.99 quantile. Is there a way

Re: [R] session logging

Maybe I'm missing something in the history help file. I read it and subsequently tried savehistory(trythis.txt) but the resulting file does not contain the text that scrolled past the console. It contains only a few uninteresting timestamps. Is there some more useful application of history() or

[R] if confusion

Simple question: Why doesn't the following work? Or what 'R' rule am I missing? tclass - Testing 1 2 3 if(tclass == Testing 1 2 3) { cat(Testing, tclass, \n) } else { cat(tclass, \n) } I get an error 'else' is unexpected. Thank you. Kevin

[R] hist2d

I'd like to use the hist2d function. What library contains it? Thanks! Simon Firestone Economist, Fair Lending Enforcement Division of Consumer and Community Affairs Federal Reserve Board Washington, DC 20551 simon.b.firest...@frb.gov (202) 785-6056 [[alternative HTML version

Re: [R] if confusion

Because the line if... Is syntactically complete. Move the else statement to the line above it like tclass - Testing 1 2 3 if ( tclass == Testing 1 2 3 ) { cat ( Testing , tclass , \n ) } else { cat ( tclass , \n ) } david -- David  

Re: [R] if confusion

This is because you could imagine a program like the following: tclass - Testing 1 2 3 if(tclass == Testing 1 2 3) { cat(Testing, tclass, \n) } cat(more stuff) However, in your program, the next line is else, which is not a command you can just put in the middle of your program. The

Re: [R] if confusion

On 4/08/2009, at 8:40 AM, rkevinbur...@charter.net wrote: Simple question: Why doesn't the following work? Or what 'R' rule am I missing? tclass - Testing 1 2 3 if(tclass == Testing 1 2 3) { cat(Testing, tclass, \n) } else { cat(tclass, \n) } I get an error 'else' is unexpected.

Re: [R] hist2d

Hi, On Aug 3, 2009, at 4:44 PM, simon.b.firest...@frb.gov wrote: I'd like to use the hist2d function. What library contains it? Thanks! Try: R RSiteSearch('hist2d') Several signs point to the gplots package ... try that one. -steve -- Steve Lianoglou Graduate Student: Computational

Re: [R] Installing package for i386

On Aug 3, 2009, at 3:01 PM, Sarah Hawley wrote: I am trying to install a package from a collaborator that was given to me as a zip file. I'm running R 2.9.0 using Mac OSX 10.5.7 [R.app GUI 1.28 (5395) i386-apple-darwin8.11.1] foo.zip contains all the usual package directories (man, R, R-ex,

Re: [R] Installing package for i386

On Aug 3, 2009, at 4:01 PM, Sarah Hawley wrote: I am trying to install a package from a collaborator that was given to me as a zip file. I'm running R 2.9.0 using Mac OSX 10.5.7 [R.app GUI 1.28 (5395) i386-apple-darwin8.11.1] foo.zip contains all the usual package directories (man, R, R-ex,

Re: [R] Open data-file in R - from Rserve

Do I have to use to openFile method? How can I assign this file to a dataframe? joaodaniel wrote: How can I open a data file from the client on the RServe server using Java? -- View this message in context:

Re: [R] session logging

Jacob, As far as I am aware, the OSX GUI has no logging/transcript feature per se, short of saving the console output to an external file. The history file will, as you note, only contain the commands entered at the R prompt, not the additional output. One of the solutions would be to

[R] Help with reshaping data.frame

I'm having trouble reshaping a data.frame from long to wide. (I think that's the right terminology; feel free to educate me.) I've looked at the reshape function and package and plyr package, but I can't quite figure out how to do this after a dozen variations. I have a data.frame with more

[R] What does this error message mean?

Hi, I used R to run a linear regression and keep getting the following error message. I do not understand it very well. Anyone can help out? Thanks. Error in storage.mode(y) - double : invalid to change the storage mode of a factor In addition: Warning message: In model.response(mf, numeric) :

Re: [R] Variable names as inputs...

That's a great guide!!! Thanks I guess, I wasn't 2 clear about my question. I have succeed naming my vars, but I can't seem to be able to create a c() of objects named after the string. So far I have MyList$lm1 MyList$lm2, ... but can't use them as an input to AIC() that is, AIC(Mylist)

[R] Help with data type

Hi there, Using a quantmod function, I calculate the daily change between two points in a time series. However, I don't think I am using the data correctly. Code: getSymbols(^GSPC, src=yahoo) CloseData - Cl(GSPC) Delta - diff(CloseData, lag=1) for (i in 3:length(Delta)) { if

[R] Some SQL Challenges

Hi, Im trying to complete a list of jobs using SQL Querries and some if else commands but im stucked in some steps. Could any of you give me some help? -in COLUMN1 change the format 20JAN2000:00:00:00 to 20JAN2000 and exclude every row that date is different from 20. -extract the first

[R] Merge partially duplicated rows

Dear all, I have a dataset, and I wanted to merge the rows with duplicated IDs by calculating the means or medians from the duplicate rows. I tried using the command duplicated(x), but it only tells where the duplicated rows are. Any suggestions will be appreciated. -- View this message in

[R] RServe - How to use 'createReference' method?

I need to input a txt, or xls, file from a client to R, using RServe. From what I've been reading, the best way to do this, is using the 'createReference' method, from REngine package. But I couldn't find any documents exemplifying it's use. I got to upload a file from java? And then? How do i

Re: [R] how do i retain decimal values

Thanks Peter, the options(digits=) works. The matrix I send as example is the x values in the correlation I will do. The Y values are different and are not tied. Thanks again. manisha -Original Message- From: P Ehlers [mailto:ehl...@ucalgary.ca] Sent: Sunday, August 02, 2009 10:32 AM

  1   2   >