Re: [R] poisson fit for histogram

Hi, see: http://cran.r-project.org/doc/contrib/Ricci-distributions-en.pdf Regards, Vito Thomas Isenbarger isen at plantpath.wisc.edu wrote: I haven't been an R lister for a bit, but I hope to enlist someone's help here. I think this is a simple question, so I hope the answer is not

Re: [R] Clustered standard errors in a panel

Have you considered lmer in library(lme4)? If you are interested in this, you may want to check the article by Doug Bates in the latest R news, www.r-project.org - Documentation: Newsletter. spencer graves Thomas Davidoff wrote: I want to do the following: glm(y ~ x1

Re: [R] family

You want to aim to write a family for R, not find the equivalents of S constructs -- they are different and so exact equivalents do not exist. In particular, an R family has several components which an S family does not. There are lots of example families for you to follow (e.g. see ?family

[R] again, a question between R and C++

Dear R Users, I want to make a call from R into C++. My inputs are List1, List2, List3, IntegerID. The amount of elements of the lists and their type depend on IntegerID. Typical elements of a given list can be vectors, doubles, and even other lists. I want to return also a list (whose nature

Re: [R] cut in R

Steve Su [EMAIL PROTECTED] writes: Dear All, I wonder whether it is still valid to use the following R code for cut. All I have done is changed: if (is.na(breaks) | breaks 2) to: if (is.na(breaks) | breaks 1) so that it covers interval of 1? It seems okay for my

[R] dpill in KernSmooth package

Hi, just a quick question does dpill computes the bandwidth or half-bandwidth? The help says bandwidth, but in the literature there is often confusion between the bandwidth and half-bandwidth. thanks, Giacomo [[alternative HTML version deleted]]

[R] Question about 'text' (add lm summary to a plot)

I would like to annotate my plot with a little box containing the slope, intercept and R^2 of a lm on the data. I would like it to look like... ++ | Slope : 3.45 +- 0.34 | | Intercept : -10.43 +- 1.42 | | R^2 : 0.78 |

[R] bubble.plot() - standardize size of unit circle

Hello, I wrote a wrapper for symbols() that produces a bivariate bubble plot, for use when plot(x,y) hides multiple occurrences of the same x,y combination (e.g. if x,y are integers). Circle area ~ counts per bin, and circle size is controlled by 'scale'. Question: how can I automatically make

Re: [R] Is it possible to create highly customized report in *.xls format by using R/S+?

Thank you all for the replies. It is very eye-opening for me. I probably need something like RDCOMClient. I've tried it last night. Very nice package!!! On 7/20/05, Gabor Grothendieck [EMAIL PROTECTED] wrote: Here is an example where R is the client and Excel is the server so that R is issuing

Re: [R] again, a question between R and C++

Jordi, The place to ask this question is probably the r-devel list; it's a little too heavy for r-help. This is fairly easy to do using the .Call interface. Have a look at lapply2 in the Writing R Extensions manual. http://cran.r-project.org/doc/manuals/R-exts.html#Evaluating-R-expressio

Re: [R] The steps of building library in R 2.1.1

Ivy_Li wrote: Dear All, With the warm support of every R expert, I have built my R library successfully. Especially thanks: Duncan Murdoch Gabor Grothendieck Henrik Bengtsson Uwe Ligges You are welcome. The following is

[R] Problem with read.table()

Dear all, I have encountered a strange problem with read.table(). When I try to read a tab delimited file I get an error message for line 260 not being equal to 14 (see below). Using count.fields() suggests that a number of lines have length not equal to 14, but not 260. Looking at the

[R] R graphics

Hi I am trying to set up 16 graphs on one graphics page in R. I have used the mfrow=c(4,4) command. However I get a lot of white space between each graph. Does anyone know how I can reduce this? Thanks Sam __ R-help@stat.math.ethz.ch mailing list

Re: [R] Is it possible to create highly customized report in *.xls format by using R/S+?

So your conclusion is that the only choice is to make mistakes and get in trouble. (That's what Excel excels at.) Two options I haven't seen mentioned are: 1. Create your deliverables in HTML format, and change the extension from .htm to .xls; Excel will import them automatically. The way the

Re: [R] Question about 'text' (add lm summary to a plot)

Dear Dan I can only help you with your third problem, expression and paste. You can use: plot(1:5,1:5, type = n) text(2,4,expression(paste(Slope : , 3.45%+-%0.34, sep = )), pos = 4) text(2,3.8,expression(paste(Intercept : , -10.43%+-%1.42)), pos = 4) text(2,3.6,expression(paste(R^2,: , 0.78, sep

[R] RandomForest question

Hello, I'm trying to find out the optimal number of splits (mtry parameter) for a randomForest classification. The classification is binary and there are 32 explanatory variables (mostly factors with each up to 4 levels but also some numeric variables) and 575 cases. I've seen that although

[R] heatmap color distribution

Hi all, I've got a set of gene expression data, and I'm plotting several heatmaps for subsets of the whole set. I'd like the heatmaps to have the same color distribution, so that comparisons may be made (roughly) across heatmaps; this would require that the color distribution and

Re: [R] Chemoinformatic people

Just with R, or via another tool integrating R, such as Pipeline Pilot? best, -tony On 7/20/05, Frédéric Ooms [EMAIL PROTECTED] wrote: Dear colleague, Just an e-mail to know if they are people working in the field of chemoinformatic that are using R in their work. If yes I was wondering if

Re: [R] R graphics

Sam Baxter wrote: Hi I am trying to set up 16 graphs on one graphics page in R. I have used the mfrow=c(4,4) command. However I get a lot of white space between each graph. Does anyone know how I can reduce this? Thanks Sam __

Re: [R] The steps of building library in R 2.1.1

I think you have been using R too long. Something like this is very much needed. There are two problems: 1. the process itself is too complex (need to get rid of perl, integrate package development tools with package installation procedure [it should be as easy as downloading a

Re: [R] R graphics

Sam Baxter wrote: Hi I am trying to set up 16 graphs on one graphics page in R. I have used the mfrow=c(4,4) command. However I get a lot of white space between each graph. Does anyone know how I can reduce this? Thanks Sam Two options: 1. play around with the `mar' parameter

Re: [R] Chemoinformatic people

I am looking for both. Fred -Original Message- From: A.J. Rossini [mailto:[EMAIL PROTECTED] Sent: Thursday, July 21, 2005 3:36 PM To: Frédéric Ooms Cc: r-help@stat.math.ethz.ch Subject: Re: [R] Chemoinformatic people Just with R, or via another tool integrating R, such as Pipeline

Re: [R] bubble.plot() - standardize size of unit circle

Thanks- 'sizeplot' didn't come up in any of my searches. Dan --- Jim Lemon [EMAIL PROTECTED] wrote: Dan Bebber wrote: Hello, I wrote a wrapper for symbols() that produces a bivariate bubble plot, for use when plot(x,y) hides multiple occurrences of the same x,y combination (e.g.

Re: [R] Problem with read.table()

On Thu, 21 Jul 2005, Kristian Skrede Gleditsch wrote: Dear all, I have encountered a strange problem with read.table(). Most `strange problems' are user error, so please try not to blame your tools. When I try to read a tab delimited file I get an error message for line 260 not being

Re: [R] Clustered standard errors in a panel

No, he wants to fit a glm() and get the right standard errors. For linear models the best way to do this is to model some random effects, but doing so in glm changes the meanings of the parameters. To estimate the same parameters you want to use the sandwich standard errors variously

[R] cutomized link function in R

Hello! I am trying to run my S+ code in R (version 2.1.0). I've created a customized link function, namely my.binomial where parameter theta has to be given. I'm considering theta to be .05. Unfortunately, R is giving an error (I had used MASS in adjusting the S+ code to R). I would very much

Re: [R] Chemoinformatic people

I am just curious why I always want to have a position like that but never find one. Am I lazy or unlucky for job huntering?(^%$$%*^( weiwei On 7/21/05, Frédéric Ooms [EMAIL PROTECTED] wrote: I am looking for both. Fred -Original Message- From: A.J. Rossini [mailto:[EMAIL

Re: [R] heatmap color distribution

You can use the breaks argument in image to do this. (You don't specify a function you're using, but other heatmap functions probably have a similar parameter.) Look across all your data, figure out the ranges you want to have different colors, and specify the appropriate break points in each

[R] Rprof fails in combination with RMySQL

Dear R community, I tried to optimized my R code by using Rprof. In my R code I'm using MySQL database connections intensively. After a bunch of queries R fails with the following error message: Error in .Call(RS_MySQL_newConnection, drvId, con.params, groups, PACKAGE = .MySQLPkgName) :

Re: [R] Problem with read.table()

I don't really understand it, but the problem seems to come down to the presence of apostrophes (single right quotes ') in the text strings. The first of these occurs in line 149 (not counting the header line). If one tries to scan just that line, one gets a vector of length 10. Fields 10 to

Re: [R] Chemoinformatic people

I don't. THere is an address an email at novartis in the ASA directory ID 068970 NameAnthony J. Rossini Company Novartis Pharma AG Address Biostatistics WSJ-27.1.012 City State Zip CH-4002 Basel Country Switzerland Phone (206) 543-2005 Email

Re: [R] Chemoinformatic people

I know of a good number of companies who use R via pipeline pilot (and have looked into it a bit recently), but not R by itself. One of the big I wish items that I've got is seemless handling of large data. Some of the RDBMS will do it, but not quite seemlessly. SPLUS 7.0 does it for a

Re: [R] RandomForest question

[EMAIL PROTECTED] wrote: Hello, I'm trying to find out the optimal number of splits (mtry parameter) for a randomForest classification. The classification is binary and there are 32 explanatory variables (mostly factors with each up to 4 levels but also some numeric variables) and 575

Re: [R] heatmap color distribution

Thanks for the reply. As I understand it, breaks only controls the binning. The problem I'm having is that each subset heatmap has slightly different min and max log2 intensities. I'd like the colors to be based on the overall (complete set) max and min, not the subsets' max and min --

[R] debian vcd package

[Apologies if you have already read this message sent from another email address] Hi R-Help, I have been using R in Linux (Debian) for the past month. The usual way I install packages is through apt. Recently, a new packages vcd became available on CRAN. I tried installing it today and

Re: [R] heatmap color distribution

Breaks affects the binning into colors. Try this. Assume that temp is one of your data sets. It's values are restricted to 0.25 - 0.75, and we'll assume that the full data set goes from 0 to 1. temp - matrix(runif(60, 0.25, 0.75), nc = 6) breaks - seq(from = 0, to = 1, length = 11)

Re: [R] Rprof fails in combination with RMySQL

I think you're barking up the wrong tree. Optimize the MySQL code separately from optimizing the R code. A very nice reference about the former is http://highperformancemysql.com/. Also, if possible, do everything in MySQL. hth, b. -Original Message- From: Thieme, Lutz [mailto:[EMAIL

Re: [R] The steps of building library in R 2.1.1

An article like that would be really great. On 7/21/05, Uwe Ligges [EMAIL PROTECTED] wrote: Gabor Grothendieck wrote: I think you have been using R too long. Something like this is very much needed. There are two problems: 1. the process itself is too complex (need to get rid of

[R] principal component analysis in affy

Hi, I have been using the prcomp function to perform PCA on my example microarray data, (stored in metric text files) which looks like this: 1a 1b 1c 1d 1e 1f ...4r 4s 4t g11.2705 1.2766

Re: [R] RandomForest question

Hi, I found the following lines from Leo's randomForest, and I am not sure if it can be applied here but just tried to help: mtry0 = the number of variables to split on at each node. Default is the square root of mdim. ATTENTION! DO NOT USE THE DEFAULT VALUES OF MTRY0 IF YOU WANT TO OPTIMIZE THE

Re: [R] The steps of building library in R 2.1.1

On 7/21/2005 9:43 AM, Gabor Grothendieck wrote: I think you have been using R too long. Something like this is very much needed. There are two problems: 1. the process itself is too complex (need to get rid of perl, integrate package development tools with package installation

Re: [R] The steps of building library in R 2.1.1

On 7/21/2005 10:29 AM, Uwe Ligges wrote: Gabor Grothendieck wrote: I think you have been using R too long. Something like this is very much needed. There are two problems: 1. the process itself is too complex (need to get rid of perl, integrate package development tools with package

Re: [R] heatmap color distribution

-Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Behalf Of Jacob Michaelson Sent: 21 July 2005 12:26 To: r-help@stat.math.ethz.ch Subject: [R] heatmap color distribution Hi all, I've got a set of gene expression data, and I'm plotting several heatmaps

[R] output of variance estimate of random effect from a gamma frailty model using Coxph in R

Hi, I have a question about the output for variance of random effect from a gamma frailty model using coxph in R. Is it the vairance of frailties themselves or variance of log frailties? Thanks. Guanghui __ R-help@stat.math.ethz.ch mailing list

[R] Concatenate 2 functions

hi all I need to concatenate 2 functions into one like temp-1:1000 for(i=0;i1000;i++) { func- func function(beta) dweibull(temp[i],beta,eta) } Any idee on this? thks guillaume. // Webmail Oreka : http://www.oreka.com

Re: [R] R:plot and dots

On Thu, 2005-07-21 at 16:18 +0200, Clark Allan wrote: hi all a very simple question. i have plot(x,y) but i would like to add in on the plot the observation number associated with each point. how can this be done? / allan If you mean the unique observation number associated with

[R] reorder bug in heatmap.2?

I want to plot a heatmap without reordering the columns. This works fine in heatmap: heatmap(meanX[selected,], col=cm.colors(256), Colv=NA) But in heatmap.2 I get: heatmap.2(meanX[selected,], col=cm.colors(256), Colv=NA) Error in if (!is.logical(Colv) || Colv) ddc - reorder(ddc, Colv) :

Re: [R] debian vcd package

On Thu, Jul 21, 2005 at 03:55:29PM +0100, Peter Ho wrote: [Apologies if you have already read this message sent from another email address] Hi R-Help, I have been using R in Linux (Debian) for the past month. The usual way I install packages is through apt. Recently, a new packages

Re: [R] The steps of building library in R 2.1.1

Duncan Murdoch wrote: On 7/21/2005 10:29 AM, Uwe Ligges wrote: Gabor Grothendieck wrote: I think you have been using R too long. Something like this is very much needed. There are two problems: 1. the process itself is too complex (need to get rid of perl, integrate package

Re: [R] Concatenate 2 functions

[EMAIL PROTECTED] wrote: hi all I need to concatenate 2 functions into one like temp-1:1000 for(i=0;i1000;i++) { func- func function(beta) dweibull(temp[i],beta,eta) } Please read An Introduction to R. Please read the posting guide. What do you expect to be in func? This is

Re: [R] R:plot and dots

Clark Allan wrote: hi all a very simple question. i have plot(x,y) but i would like to add in on the plot the observation number associated with each point. how can this be done? See ?text Uwe Ligges / allan

[R] About object of class mle returned by user defined functions

Hi, There is something I don't get with object of class mle returned by a function I wrote. More precisely it's about the behaviour of method confint and profile applied to these object. I've written a short function (see below) whose arguments are: 1) A univariate sample (arising from a

[R] R:plot and dots

hi all a very simple question. i have plot(x,y) but i would like to add in on the plot the observation number associated with each point. how can this be done? / allan__ R-help@stat.math.ethz.ch mailing list

Re: [R] The steps of building library in R 2.1.1

On 7/21/05, Duncan Murdoch [EMAIL PROTECTED] wrote: On 7/21/2005 9:43 AM, Gabor Grothendieck wrote: I think you have been using R too long. Something like this is very much needed. There are two problems: 1. the process itself is too complex (need to get rid of perl, integrate

[R] opening RDB files

Hi all, I've recently upgraded to R version 2.1.1 and when trying to inspect the contents of many packages in the library (for instance library\MASS\R) I've realized wordpad, or the notepad, won't open them since they have *.RDB and *.RDX extensions which these editors cannot recognize.

Re: [R] opening RDB files

Emili Tortosa-Ausina wrote: Hi all, I've recently upgraded to R version 2.1.1 and when trying to inspect the contents of many packages in the library (for instance library\MASS\R) I've realized wordpad, or the notepad, won't open them since they have *.RDB and *.RDX extensions which

[R] normal reference intervals

I am interested in calculating Age-Specific normal reverence intervals, using non-parametric methods - or ideally something called the LMS method (which as I understand it uses cubic splines fitted to the data). Any packages in R that you think might help me? Any other advice gratefully received.

Re: [R] opening RDB files

Emili Tortosa-Ausina wrote: Hi all, I've recently upgraded to R version 2.1.1 and when trying to inspect the contents of many packages in the library (for instance library\MASS\R) I've realized wordpad, or the notepad, won't open them since they have *.RDB and *.RDX extensions which

Re: [R] colnames

Hi Adai, Your diagnosis is absolutely right; class(r1) returned data.frame and your suggested solution worked perfectly. Your assumption is also right; both x and y are positive. If I want to compare the performance of the my old function with yours, are there some functions in R I could use

Re: [R] RandomForest question

See the tuneRF() function in the package for an implementation of the strategy recommended by Breiman Cutler. BTW, randomForest is only for the R package. See Breiman's web page for notice on trademarks. Andy From: Weiwei Shi Hi, I found the following lines from Leo's randomForest,

Re: [R] RandomForest question

From: [EMAIL PROTECTED] Hello, I'm trying to find out the optimal number of splits (mtry parameter) for a randomForest classification. The classification is binary and there are 32 explanatory variables (mostly factors with each up to 4 levels but also some numeric variables) and

Re: [R] R graphics

Sam Baxter [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] I am trying to set up 16 graphs on one graphics page in R. I have used the mfrow=c(4,4) command. However I get a lot of white space between each graph. Does anyone know how I can reduce this? The default par()$mar is

Re: [R] colnames

See help(system.time). On Thu, 2005-07-21 at 17:56 +0100, Gilbert Wu wrote: Hi Adai, Your diagnosis is absolutely right; class(r1) returned data.frame and your suggested solution worked perfectly. Your assumption is also right; both x and y are positive. If I want to compare the

Re: [R] Problem with read.table()

Thanks to all who responded to my earlier message. The problem lies in that apostrophes (i.e., ') in some of the text fields are read as quotes. The file can be read without problems setting quotes= in read.table. Incidently, read.delim() also works, even without setting quotes= explicitly.

Re: [R] help with a xyplot legend

On 7/20/05, Ronaldo Reis-Jr. [EMAIL PROTECTED] wrote: Hi, I try to put a legend in a xyplot graphic. xyplot(y~x|g,ylim=c(0,80),xlim=c(10,40),as.table=T,layout=c(2,3), ylab=Número de machos capturados,xlab=expression(paste(Temperatura (,degree,C))), key=list(corner=c(0,0),x=0, y=0,

Re: [R] The steps of building library in R 2.1.1

On 7/21/2005 12:03 PM, Gabor Grothendieck wrote: 2. there is too much material to absorb just to create a package. The manuals are insufficient. The first sentence here is basically a repetition of the process is too complex. I think the second sentence is incorrect. Could you

[R] problem running R with perl Statistics::R

__ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

[R] Using Perl Statistics::R module, not running the

__ R-help@stat.math.ethz.ch 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] Chemoinformatic people

Sure, but Luke, I DO NOT currently use R at work... (now, that's not to say I won't be using it in a few months, but currently...). best, -tony On 7/21/05, Luke Tierney [EMAIL PROTECTED] wrote: I don't. THere is an address an email at novartis in the ASA directory ID 068970

Re: [R] The steps of building library in R 2.1.1

Gabor Grothendieck wrote: I think you have been using R too long. Something like this is very much needed. There are two problems: 1. the process itself is too complex (need to get rid of perl, integrate package development tools with package installation procedure [it should be

Re: [R] Question about 'text' (add lm summary to a plot)

On Thu, 21 Jul 2005, Christoph Buser wrote: Dear Dan I can only help you with your third problem, expression and paste. You can use: plot(1:5,1:5, type = n) text(2,4,expression(paste(Slope : , 3.45%+-%0.34, sep = )), pos = 4) text(2,3.8,expression(paste(Intercept : , -10.43%+-%1.42)), pos = 4)

Re: [R] The steps of building library in R 2.1.1

Just a few thoughts... Good documentation helps everybody - the beginners and the experts (less beginner questions if there is thorough and accessible documentation. I fully appreciate that this is a volunteer effort - I'm just trying to pin down some places where we have documentation issues.

Re: [R] unable to call R t-test from Java

In my last post, I left off version info: Java: jdk1.5.0_03 R: 2.1.1 SJava: 0.68 OS: SunOs 5.8 By trying either the eval or the call method to execute a t-test results in a core dump. e.g. eval method /* produces a core */ System.err.println(eval

Re: [R] output of variance estimate of random effect from a gamma frailty model using Coxph in R

On Thu, 21 Jul 2005 [EMAIL PROTECTED] wrote: Hi, I have a question about the output for variance of random effect from a gamma frailty model using coxph in R. Is it the vairance of frailties themselves or variance of log frailties? Thanks. For a Gamma frailty model it is the variance of

Re: [R] debian vcd package

I just installed it on a Debian 2.6.8.1 and the following R platform i386-pc-linux-gnu arch i386 os linux-gnu system i386, linux-gnu status major2 minor0.1 year 2004 month11 day 15 language R By the way are you using apt-get install vcd. and if so why? just use

[R] tav files

Dear colleagues: I am using the spotfinder 2.2.3 and the .tav files generated there have 20 columns. When I exclude the 3 last columns, the .tav file can not be recognized by aroma package in R platform. What I have to do to generate .tav files with 17 columns only? Thank you in advance.

[R] Invitation to New York, Spain, and Italy; c/ba

Dear potential Speaker: On behalf of the organizing committee, I would like to extend a cordial invitation for you to submit a paper to the IPSI Transactions journal, or to attend one of the upcoming IPSI BgD multidisciplinary, interdisciplinary, and transdisciplinary conferences.

Re: [R] Question about 'text' (add lm summary to a plot)

Use bquote instead of expression, e.g. trees.lm - lm(Volume ~ Girth, trees) trees.sm - summary(trees.lm) trees.co - round(trees.sm$coefficients,2) trees.rsq - round(trees.sm$r.squared,2) plot(Volume ~ Girth, trees) text(10,60, bquote(Intercept : .(trees.co[1,1])%+-%.(trees.co[1,2])), pos = 4)

Re: [R] Question about 'text' (add lm summary to a plot)

[Note: the initial posts have been re-arranged to attempt to maintain the flow from top to bottom] Dan Bolser writes: I would like to annotate my plot with a little box containing the slope, intercept and R^2 of a lm on the data. I would like it to look like...

Re: [R] RandomForest question

From: Uwe Ligges [EMAIL PROTECTED] wrote: Hello, I'm trying to find out the optimal number of splits (mtry parameter) for a randomForest classification. The classification is binary and there are 32 explanatory variables (mostly factors with each up to 4 levels but also some

[R] A Question About Inverse Gamma

Hi R users, I am having a little problem finding the the solution to this problem in R: 1. I need to generate normal distribution of sample size 30, mean = 50, sd = 5. 2. From the statistics obtained in step 1, I need to generate the Inverse Gamma distribution. Your views and help will be

Re: [R] A Question About Inverse Gamma

[EMAIL PROTECTED] wrote: Hi R users, I am having a little problem finding the the solution to this problem in R: 1. I need to generate normal distribution of sample size 30, mean = 50, sd = 5. 2. From the statistics obtained in step 1, I need to generate the Inverse Gamma

Re: [R] R:plot and dots

Here is an example : x - rnorm(20) y - rnorm(20) plot(x, y, type=n) text(x, y, labels=as.character(1:20)) Also look into help(identify) if you want to point out specific points. Regards, Adai On Thu, 2005-07-21 at 16:18 +0200, Clark Allan wrote: hi all a very simple question.

Re: [R] The steps of building library in R 2.1.1

On 22 Jul 2005 00:01:18 +0200, Peter Dalgaard [EMAIL PROTECTED] wrote: Duncan Murdoch [EMAIL PROTECTED] writes: On 7/21/2005 9:43 AM, Gabor Grothendieck wrote: I think you have been using R too long. Something like this is very much needed. There are two problems: 1. the process

Re: [R] principal component analysis in affy

1) Please learn to wrap your emails to 72 characters per line. See http://expita.com/nomime.html 2) You might have better luck with the BioConductor folks. Their mailing list is https://stat.ethz.ch/mailman/listinfo/bioconductor 3) The affy package has many functions including some algorithms

[R] vectorising ifelse()

Hi All, is there any chance of vectorising the two ifelse() statements in the following code: for(i in gp){ new[i,1] = ifelse(srow[i]0, new[srow[i],zippo[i]], sample(1:100, 1, prob =Y1, rep = T)) new[i,2] = ifelse(drow[i]0, new[drow[i]0,zappo[i]], sample(1:100, 1, prob =Y1, rep = T)) }

[R] Re Randomization test for interaction effect

Dear Pedro, How to test for an interaction--or, even, how to pose the question of an interaction--in randomization-based inference is not at all obvious. And, in the permutation test context reliance has been placed on the exchangeability of (estimated) residuals under an additive, homoscedastic

[R] help: pls package

Hello, I have a data set with 15 variables (first one is the response) and 1200 observations. Now I use pls package to do the plsr with cross validation as below. trainSet = as.data.frame(scale(trainSet, center = T, scale = T)) trainSet.plsr = mvr(formula, ncomp = 14, data = trainSet, method =

[R] help: pls package

Hello, I have a data set with 15 variables (first one is the response) and 1200 observations. Now I use pls package to do the plsr with cross validation as below. trainSet = as.data.frame(scale(trainSet, center = T, scale = T)) trainSet.plsr = mvr(formula, ncomp = 14, data = trainSet, method =

Re: [R] :)

olá - Original Message - From: [EMAIL PROTECTED] To: [EMAIL PROTECTED] Sent: Thursday, July 21, 2005 3:04 PM Subject: :) pq nao me liga?? __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

[R] find confounder in covariates

Hi, I was wondering if there is a way, or function in R to find confounders. For istance, a = sample( c(1:3), size=10,replace=T) X1 = factor( c('A','B','C')[a] ) X2 = factor( c('Aa','Bb','Cc')[a] ) Xmat = data.frame(X1,X2,rnorm(10),rnorm(10)) dimnames(Xmat)[[2]] = c('z1','z2','z3','y')

Re: [R] A Question About Inverse Gamma

On Thu, 21 Jul 2005, Sundar Dorai-Raj wrote: [EMAIL PROTECTED] wrote: Hi R users, I am having a little problem finding the the solution to this problem in R: 1. I need to generate normal distribution of sample size 30, mean = 50, sd = 5. 2. From the statistics obtained in step 1, I