[R] Takagi-Sugeno fuzzy model

Dear R-help. I am learning fuzzy algorithm, at this moment about Takagi-Sugeno fuzzy model. I have been searching on mailing list r-help about this, but I didn't find any list about Takagi-Sugeno fuzzy model. I want to know, are there any R code (or library/package) for Takagi-Sugeno fuzzy model.

Re: [R] fitting a lognormal distribution using cumulative probabilities

Dear Prof Ripley and Dimitris (cc R-list), thank you very much for your very insightful responses. I've been checking how to use survreg{survival} to fit a left-censored lognormal distribution, and I was surprised to find that results are exactly the same as with fitdistr{MASS}. Here is an exampl

Re: [R] fitting a lognormal distribution using cumulative probabilities

Your example code is asserting that the events occurred at the times in 'obs.time', not before those times. Surv(event = 1) means uncensored. If you try event = 0, the fitter diverges towards an exact fit, as I said it should. Yes, you will get a good fit to the ECDF, but you are modelling the

[R] Newbie: Measuring distance between clusters.

Hi, I had 1 genes, and I clustered them using K-means clustering in R. kc<-kmeans(data.sub,7) kc n cluster sum of squares by cluster: [1] 60631.76 135886.19 159049.71 101783.27 90040.72 183335.60 158867.81 Available components: [1] "cluster" "centers" "withinss" "size" I am very new t

Re: [R] fitting a lognormal distribution using cumulative probabilities

Dear Professor Ripley, thanks again for your mail and useful comments. Surv(event = 1) means uncensored. If you try event = 0, the fitter diverges towards an exact fit, as I said it should. Sorry, I misunderstood this. (BTW, using a .com address suggests you are a COMpany and in the absence of

[R] Hiding a function

Hi the list Is it possible to 'hide' a function from the user ? I cut a big fonction in sub function and I would like to hide the sub function, just like if I declare them in the big function : -- a <- function(x){ b <- function(y){y^2} d <- function(y){y^3} b(x)+d(x)+2 } a(2

Re: [R] addtable2plot(plotrix)

Markus Didion wrote: > Dear List > > adding a table to a plot using addtable2plot(plotrix) does not seem to > work when using a logarithmic axis. The table is then reduced to one > line. Is there an argument to indicate that a log-scale is used, or an > alternative to add a bunch of informati

Re: [R] Hiding a function

On 23/02/2008 5:15 AM, Christophe Genolini wrote: > Hi the list > > Is it possible to 'hide' a function from the user ? I cut a big > fonction in sub > function and I would like to hide the sub function, just like if I > declare them > in the big function : > > -- > a <- function(x

Re: [R] Hiding a function

Duncan Murdoch a écrit : > On 23/02/2008 5:15 AM, Christophe Genolini wrote: >> Hi the list >> >> Is it possible to 'hide' a function from the user ? I cut a big >> fonction in sub >> function and I would like to hide the sub function, just like if I >> declare them >> in the big function : >> >

[R] Set without argument

Hi the list, I am defining S4 objet. Is it possbile to define a method that change the slot of an object without using <- ? My object contain a numeric and a matrix. At some point, I would like to impute the missing value in the matrix. So I would like to use something like : -

Re: [R] Hiding a function

On 23/02/2008 5:58 AM, Christophe Genolini wrote: > Duncan Murdoch a écrit : >> On 23/02/2008 5:15 AM, Christophe Genolini wrote: >>> Hi the list >>> >>> Is it possible to 'hide' a function from the user ? I cut a big >>> fonction in sub >>> function and I would like to hide the sub function, jus

Re: [R] using subset() in data frame

On 2/22/2008 8:01 PM, Robert Walters wrote: > R folks, > As an R novice, I struggle with the mystery of subsetting. Textbook and > online examples of this seem quite straightforward yet I cannot get my > mind around it. For practice, I'm using the code in MASS Ch. 6, > "whiteside data" to analyz

Re: [R] Hiding a function

It depends what you mean by 'hiding', you can start the function names with a dot and then they are not listed by ls(), so this is kind of hiding. > .a <- function() TRUE > ls() character(0) > .a function() TRUE Personally i would not do this though. G. On Sat, Feb 23, 2008 at 11:58:57AM +0100

Re: [R] using subset() in data frame

On 2/23/2008 6:09 AM, Chuck Cleland wrote: > On 2/22/2008 8:01 PM, Robert Walters wrote: >> R folks, >> As an R novice, I struggle with the mystery of subsetting. Textbook >> and online examples of this seem quite straightforward yet I cannot >> get my mind around it. For practice, I'm using the

Re: [R] Hiding a function

Gabor Csardi a écrit : > It depends what you mean by 'hiding', Well, the main idea was the to limit the existance of a function. a( ) need b( ) but no other function will need b(). So I would have like to let b( ) exists localy only when a( ) is called, not elsewhere. Christophe > you can start

Re: [R] Hiding a function

Just to clarify, what Duncan was referring to as the alternative was nesting the definition of one function in another, e.g. look at FUNx in by.data.frame -- FUNx is available to by.data.frame but is not visible outside of by.data.frame . On Sat, Feb 23, 2008 at 6:09 AM, Duncan Murdoch <[EMAIL PRO

Re: [R] Set without argument

Christophe Genolini <[EMAIL PROTECTED]> writes: > Hi the list, > > I am defining S4 objet. Is it possbile to define a method that change > the slot of an object without using <- ? > My object contain a numeric and a matrix. At some point, I would like to > impute the missing value in the matrix.

Re: [R] Hiding a function

Ok, I saw FUNx So, I reformulate my question : is there a way, without nesting b( ) in a( ), to make b( ) available only in a( ) ? > Just to clarify, what Duncan was referring to as the > alternative was nesting the definition of one function > in another, e.g. look at FUNx in by.data.frame -- > F

Re: [R] using subset() in data frame

Erin Hodgess wrote: > Hi Robert! > > Could you please check > str(data.b) > > and see what you have for porosity? It needs to be a factor > > Thanks, > Erin > > Erin, Yes, porosity is a factor. See output below. > str(data.b) 'data.frame': 96 obs. of 7 variables: $ system : Factor w/

Re: [R] Hiding a function

On 23/02/2008 8:22 AM, Christophe Genolini wrote: > Ok, I saw FUNx > So, I reformulate my question : is there a way, without nesting b( ) in > a( ), to make b( ) available only in a( ) ? Not that I know of. It will be available to anything using the same environment as a(), which in the case of

[R] Standard method for S4 object

Hi the list, I am defining a new class MyClass. Shortly, I will submit a package with it. Before, I would like to know if there is a kind of "non official list" of what method a new S4 object have. More precisely, personnaly, I use 'print', 'summary' and 'plot' a lot. So for my new class, I def

Re: [R] using subset() in data frame

Chuck Cleland wrote: > On 2/23/2008 6:09 AM, Chuck Cleland wrote: >> On 2/22/2008 8:01 PM, Robert Walters wrote: Chuck, Thanks for the pointers on subset(). When I submit the two variants you suggested, below: > fit1 <- lm(y ~ x, subset(data.b, porosity == "macro")) Error in eval(expr, envir, e

Re: [R] Set without argument

Thanks a lot > The 'usual' way to write the above code is > > >> a <- imputeMyObj(a) >> :-( I was hopping something like the use of the superassignator <<- inside imputeMyObj (I was hoping even if I did not find it myself)... To bad. But thanks for teaching me 'usual' pratice. It is somet

[R] Aranda-Ordaz

Hi all, Does anyone know R code or SAS code for Aranda-Ordaz link family? thanks, xiao yue - [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/

Re: [R] Standard method for S4 object

Christophe Genolini wrote: >Hi the list, > >I am defining a new class MyClass. Shortly, I will submit a package with >it. Before, I would like to know if there is a kind of "non official >list" of what method a new S4 object have. >More precisely, personnaly, I use 'print', 'summary' and 'plot'

Re: [R] Hiding a function

After some private email with some high level programmer (Gabor), I have the solution and I can answer to my own question :-) The problem was 1° to define a function B( ) available only in a function A( ) (for code safeness) and 2° to define B( ) elsewhere than in A( ) (for code readability)

Re: [R] Aranda-Ordaz

I have no idea if it is helpful, but a quick google search turned up: LINKINF # Two S-PLUS functions to compute influence diagnostics ...The assumed logit link is embedded within the Aranda-Ordaz parametric # family of link functions. # Written by John Yick and Andy H. Lee ... phase.hpcc.jp/mirror

Re: [R] counting sequence mismatches

it is pretty enough for me. Thanks - Original Message From: Martin Morgan <[EMAIL PROTECTED]> To: joseph <[EMAIL PROTECTED]> Cc: r-help@r-project.org Sent: Friday, February 22, 2008 6:41:41 PM Subject: Re: [R] counting sequence mismatches One kind of ugly solution > d.f=data.fr

[R] exprs function download

Hi, I am supposed to use exprs as a function. Where can i download exprs function? I tried searching at bioconductor and seach engine but no luck. Is it located in one of the library in R? thanks. C -- View this message in context: http://www.nabble.com/exprs-function-download-tp15654560p15654

Re: [R] exprs function download

Cristophe -- 'exprs' is a generic function defined in the Biobase package, but with methods in additional packages in Bioconductor. Mostly exprs is used to extract a matrix of expression-like values from a more complicated data structure, the precise package you need depends on what data object yo

Re: [R] exprs function download

> I am supposed to use exprs as a function. Where can i download exprs > function? I tried searching at bioconductor and seach engine but no luck. Is > it located in one of the library in R? source("http://bioconductor.org/biocLite.R";) biocLite("Biobase") library(Biobase) ?exprs HTH, Tobias __

Re: [R] overdrawing a plot

> I had a similar problem, and Paul Murrell sent me a workaround: > > http://finzi.psych.upenn.edu/R/Rhelp02a/archive/104038.html > > "You can capture the ggplot drawing as a grid grob (gTree), edit that (no > drawing occurs to this point), and then draw it ..." > > Hadley realized that this i

Re: [R] plotEst

On Tue, Feb 19, 2008 at 3:53 PM, sigalit mangut-leiba <[EMAIL PROTECTED]> wrote: > Sorry to disturb, > I managed to plot 2 together with 'layer' like this: > > > > qplot(se, or, min=lcl1, max=ucl1, data=df1, geom="pointrange")+layer(data = > df2, mapping = *aes*(x = se, y = OR2,min=lcl2,max=ucl2

Re: [R] Viable Approach to Parallel R?

Hi Dan, I've had pretty good luck using Snow with with Rpvm. It's definitely not what you'd call "plug and play," but it does work. I'm using it on a single computer to just take advantage of multiple processors, and it does a pretty good job of keeping them busy. The main gotchas I've found wi

Re: [R] using subset() in data frame

Look carefully at the output and commands below. The first level of porosity is "macro " (notice the space at the end) but you are asking for "macro" (without the space). Computers are very literal, so "macro " is not equal to "macro". Try it again with the space (or renaming the level in th

Re: [R] Hiding a function

> > After some private email with some high level programmer (Gabor), I have shouldn't that rather be: high in knowledge but low in level? Best, Hans-Peter [[alternative HTML version deleted]] ___

Re: [R] Viable Approach to Parallel R?

Similar experience, with snow & MPI (LAM). Actually, plug and play. G. On Sat, Feb 23, 2008 at 10:57:22AM -0700, Eric W Anderson wrote: > Hi Dan, > > I've had pretty good luck using Snow with with Rpvm. It's definitely > not what you'd call "plug and play," but it does work. I'm using it on >

[R] Error in ma.svd(X, 0, 0) : 0 extent dimensions

Hi, I run a maanova analysis and found this message error: Error in ma.svd(X, 0, 0) : 0 extent dimensions I did a google search and found this: \item ma.svd: function to compute the sigular-value decomposition of a rectangular matrix by using LAPACK routines DEGSVD AND ZGESVD. \ite

Re: [R] using subset() in data frame

Greg Snow wrote: > Look carefully at the output and commands below. > > The first level of porosity is "macro " (notice the space at the > end) but you are asking for "macro" (without the space). Computers are > very literal, so "macro " is not equal to "macro". T Greg and Chuck, Thanks for

Re: [R] bootstrap: definition of original statistic

> Gad Abraham wrote: >> Hi, >> In the boot package, the original statistic is simply the statistic function evaluated on the original data (called t0). >> However, in Harrell et al 1996 "Multivariable prognostic models..." Stats Med vol 15, pp. 361--387, it is different (p. 372): >> The statistic f

[R] Bimodal deconvolution

Hi Everyone- After searching through posts and my favorite R-help websites I'm still confused about a problem. I have data which is bimodal in nature, but there is no clearly obvious separation between the two peaks. In programs such as Origin, I can deconvolute the two distributions and have

Re: [R] using subset() in data frame

Chuck Cleland wrote: > On 2/23/2008 9:13 AM, Robert Walters wrote: >> Chuck Cleland wrote: Chuck, For the record, I might add that that the following two variants for subsetting worked equally well: fit1 <- lm(pore.pct ~ Db, subset(data.b, porosity == "macro ")) fit1 <- lm(pore.pct ~ Db, data.

Re: [R] Bimodal deconvolution

Seth Imhoff wrote: > Hi Everyone- > > After searching through posts and my favorite R-help websites I'm still > confused about a problem. I have data which is bimodal in nature, but there > is no clearly obvious separation between the two peaks. In programs such as > Origin, I can deconvolute

[R] clarification about glm

Hello, I have a question about glm: if i have a binary covariate (unit=1,0) the reference group would be 0? (prediction for unit=1) example: dat1<-data.frame(y,unit,x1,x2) log_u <- glm(y~.,family=binomial,data=dat1) summary(log_u) Estimate Std. Error z value Pr(>|z|) (Intercept) -0.

Re: [R] clarification about glm

Hi Sigalit, yes, you can see this from the fact that the table says unit"1" meaning that it compares 1 to 0 and not vice versa. Everytime you regress on dummies the label will have added this to the original variable name. Say you have gender "male" and "female". Then "gendermale" in the label of

Re: [R] How to Include greek symbol in axis label?

Thanks for your help. I tried par(family="serif") and it worked well. I'm running R 2.6.1, on windows XP machine. font=6 appears to be similar to times new roman or serif, they appear to be very similar (to my untrained eye). Serif seems to work fine for my needs! Thanks, Ken Kenneth Tak

Re: [R] clarification about glm

Yes, that's what I thought, thanks for the clarification! Sigalit. On 2/23/08, Daniel Malter <[EMAIL PROTECTED]> wrote: > > Hi Sigalit, > > yes, you can see this from the fact that the table says unit"1" meaning > that > it compares 1 to 0 and not vice versa. Everytime you regress on dummies > the

Re: [R] plotEst

Thank you for your suggestion, now I get 2 lines but confidence limits just for the first one.. do you know what I'm doing wrong?! Thanks again, Sigalit. On 2/23/08, hadley wickham <[EMAIL PROTECTED]> wrote: > > On Tue, Feb 19, 2008 at 3:53 PM, sigalit mangut-leiba <[EMAIL PROTECTED]> > wrote: >

Re: [R] How to Include greek symbol in axis label?

On Sat, 23 Feb 2008, Kenneth Takagi wrote: > Thanks for your help. I tried par(family="serif") and it worked well. > I'm running R 2.6.1, on windows XP machine. font=6 appears to be > similar to times new roman or serif, they appear to be very similar (to > my untrained eye). Serif seems to

[R] Calculate missing value using a correlation metric

Could someone please let me know how can I calculate missing values using a correlation metric? I know that the knn function is used for calculating missing values using Euclidean distance. Thanks -- View this message in context: http://www.nabble.com/Calculate-missing-value-using-a-correlati

Re: [R] Newbie: Measuring distance between clusters.

One way to do it is to find the distances between ther "centers" (= centres in English) of the clusters. dist(kc$centers) It rather depends on how you define distances between clusters, though. There are many possibilities. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROT

Re: [R] exprs function download

RSiteSearch("exprs") will start you off. If not, ask whoever is requiring you to use it. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Keizer_71 Sent: Sunday, 24 February 2008 3:15 AM To: r-help@r-project.org Subject: [R] exprs function download

[R] where can I find source code for particle filters applied to stochastic volatilities?

Hi all, Could anybody point me to some overview/survey papers about using particle filters and sequential monte carlo methods to estimate stochastic volatilities? I couldn't find any such articles giving a big-picture view of the literature. How do these estimation methods compare to EMM and othe

Re: [R] Corrected : Efficient writing of calculation involving each element of 2 data frames

Vikas, Please provide reproducible code when posting to this list (read the posting guide). You need to vectorise your code to make it run faster. Here are some timing results for 3 different vectorised methods. All 3 methods take under 30 seconds to run on a pair of 30,000-vectors. The results