[R] Cuetion abut dunn test

Hello, I wanted to know which are the commands to perform a Dunn test, and how must the data be entered. Thanks, Nicolás Tarjeta de crédito Yahoo! de Banco Supervielle. Solicitá tu nueva Tarjeta de crédito. De tu PC directo a tu casa. www.tuprimeratarjeta.com.ar

[R] Creating Matrix

Hello, I have two data. x-c(1, 2, 1, 3, 2) y-c(3, 1, 2, 3, 5) How do i create matrix from this two. what i want is this x y 1 1 3 2 2 1 3 1 2 4 3 3 5 2 5 thanks Claire -- View this message in context: http://www.nabble.com/Creating-Matrix-tp17168173p17168173.html Sent

Re: [R] [RsR] function in nls argument -- robust estimation

Martin, Kate, Fernando, et al. Be careful bootstrapping robust estimators. The trouble is that when resampling is done with replacement, the outliers can be selected too many times which would ruin the standard error estimates. Martin is right that bootstrapping would be fine if there are not

[R] Monty Hall simulation

Is it possible to simulate the Monty Hall problem using R? If so, could someone please show me how? Thanks for any help rendered. -- View this message in context: http://www.nabble.com/Monty-Hall-simulation-tp17169235p17169235.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] Creating Matrix

See ?cbind and ?matrix. Gabor On Sat, May 10, 2008 at 03:21:26PM -0700, Claire_6700 wrote: Hello, I have two data. x-c(1, 2, 1, 3, 2) y-c(3, 1, 2, 3, 5) How do i create matrix from this two. what i want is this x y 1 1 3 2 2 1 3 1 2 4 3 3 5 2 5

Re: [R] Random number generation

Dennis, I assume that there is a set.seed() somewhere in your script, possibly in something you source()d (hopefully not in anything library()d). Have you tried successively removing/commenting parts of the script before the sample() command until the problem goes away? That way you should

Re: [R] Monty Hall simulation

cirrus74 [EMAIL PROTECTED] [Sun, May 11, 2008 at 03:44:46AM CEST]: Is it possible to simulate the Monty Hall problem using R? If so, could someone please show me how? Thanks for any help rendered. The kind of simulation, as any thinking about this seeingly paradoxical situation, depends on

[R] loess and locpoly

Dear list, I've got a question concerning difference between loess and locpoly. I have to use a plug-in method to chose a bandwith so I take locpoly method to fit a curve. My problem is:I know how to get predicted values in loess: m=loess(y~x) y_fitted=predict(m). But how to get the

Re: [R] Is there in R a function equivalent to the mround, as found in most spreadsheets?

Hello Luca. Dr. Ottorino-Luca Pantani wrote: ... c(1803.02, 193.51, 3.47) Each solution is to be taken with 3 different pipettes (5000, 250 and 10 µL Volume max) and each of those delivers volumes in steps of 50 µL, 5 µL or 1µL, respectively Since the above values would

[R] Problems with blanks and umlauts in filenames with Rserve

Dear R developers, I'm trying to load files by means of Rserve. If the files have blanks in it names or german umlauts the loading failes(for example when using the rgdal lib with readGDAL() ). In the R application without Rserve this works. Is there a general encoding recipe or switch that

Re: [R] Bar Charts

hoogeebear wrote: Hi, I created some bar charts. My first one is concerned with males, and my second concerned with females. Is there a way I can put the charts into one chart? There is 2 different columns in each file. Here is my new file containing males and females: gender,familar

Re: [R] Random number generation

Stephan Kolassa wrote: Have you tried successively removing/commenting parts of the script before the sample() command until the problem goes away? That way you should be able to pinpoint the offending script command. Hi, This brings up a question I have .. is there a way to do *block*

Re: [R] Random number generation

On Sun, 11 May 2008, Esmail Bonakdarian wrote: Stephan Kolassa wrote: Have you tried successively removing/commenting parts of the script before the sample() command until the problem goes away? That way you should be able to pinpoint the offending script command. Hi, This brings up a

Re: [R] substitute in graphics - Uwe's help desk

Hello, I still have my problem. I couldn't make Uwe Ligges's example work. I wonder why :-( The following doesn't work : param.list=list(mu1=0,mu2=0,s1=3,s2=2,s3=2,s4=4) plot(1:8,type=n) text(8,3,adj=1,labels=substitute(with * mu == bgroup((,atop(mu1,mu2),)) * , * Sigma[x] ==

Re: [R] Random number generation

Hello there, Prof Brian Ripley wrote: On Sun, 11 May 2008, Esmail Bonakdarian wrote: Stephan Kolassa wrote: Have you tried successively removing/commenting parts of the script before the sample() command until the problem goes away? That way you should be able to pinpoint the offending

Re: [R] Creating Matrix

Here is how you can apply the mat function mentioned by Cassardi, x-c(1, 2, 1, 3, 2) y-c(3, 1, 2, 3, 5) mat-matrix(c(x,y),5,2) ##the first parameter gives the data vector wich is filled columnwise in the matrix, then comes the row and column dimensions) colnames(mat)-c(x,y) you can

Re: [R] Creating Matrix

how do i calculate the p-value of trend test power oquote author=Claire_6700 Hello, I have two data. x-c(1, 2, 1, 3, 2) y-c(3, 1, 2, 3, 5) 1. How do i create matrix from this two. what i want is this x y 1 1 3 2 2 1 3 1 2 4 3 3 5 2 5 2. what is the best way to use

Re: [R] Compact Indicator Matrices

On Sat, May 10, 2008 at 5:27 AM, amarkos [EMAIL PROTECTED] wrote: An indicator matrix is a binary matrix with orthogonal columns whose rows sum to 1. A row of this matrix could be [0 1 0 0]. My problem is to group the similar rows (profiles) so that to create a compact form of the matrix. I'm

[R] superscript text on graph legend

Is that possible to create superscript text on the graph legend, for example to put cm2 (centimeter square) on the legend. Please show me how to do it. Thanks. -- Agus Susanto [[alternative HTML version deleted]] __ R-help@r-project.org

Re: [R] subset rows in two dataframes

Not exactly. I need something to subset ONLY rows common to both dataframes. In the provided example, dat1 and dat2 have no common rows so I would expect: [1] v1 v2 0 rows (or 0-length row.names) But I can´t do it... On Sun, 11 May 2008 10:07:25 -0400, Zhuanshi He [EMAIL PROTECTED] said:

[R] combining bar and column graphs?

Hi, I'm hoping to find out whether R, or an R add-on, can generate a particular type of graph. And, more basically, whether such a type of graph even makes sense. I'm looking for something resembling both a column chart and a bar chart, where the basic visual unit is a solid rectangle of

Re: [R] subset rows in two dataframes

It is giving you exactly what you are asking for. You asked first which of dat1$v1 were in dat2$v; you got a TRUE on the second value (2006-05-03): dat1$v1 %in% dat2$v1 [1] FALSE TRUE TRUE FALSE TRUE FALSE TRUE TRUE FALSE FALSE FALSE FALSE FALSE You then asked for which of dat1$v2 were in

Re: [R] subset rows in two dataframes

Here is one way to compare the entire rows between the data frames: x1 - do.call(paste, dat1) x2 - do.call(paste, dat2) dat1[x1 %in% x2,] [1] v1 v2 0 rows (or 0-length row.names) x1 [1] 2006-01-03 7312.5 2006-05-03 3352.5 2006-05-04 4252.5 2006-05-11 3825 [5] 2006-05-12 2700 2006-05-16

Re: [R] substitute in graphics - Uwe's help desk

Gustave Lefou wrote: Hello, I still have my problem. I couldn't make Uwe Ligges's example work. I wonder why :-( The following doesn't work : param.list=list(mu1=0,mu2=0,s1=3,s2=2,s3=2,s4=4) plot(1:8,type=n) text(8,3,adj=1,labels=substitute(with * mu == bgroup((,atop(mu1,mu2),)) * , *

[R] problems installing rJava

Hello, I'm having trouble installing rJava on R-2.6.1 / Suse10.3. What could be wrong? Regards, Juan Pablo This is the output of install.packages(rJava, dependencies=T) --- Please select a CRAN mirror for use in this session --- Loading Tcl/Tk interface ... done probando la URL

Re: [R] subset rows in two dataframes

Dear Jim, Maybe u want this, subset(dat2, time1 %in% dat2$v1 time2 %in% dat2$v1) v1 v2 2 2006-05-09 7065.0 3 2006-05-04 3622.5 5 2006-07-14 3532.5 7 2006-05-12 6480.0 8 2006-05-17 4612.5 15 2006-07-05 4837.5 16 2006-07-06 3352.5 18 2006-07-24 6772.5 20 2006-07-18 5625.0

Re: [R] Compact Indicator Matrices

On May 11, 4:47 pm, Douglas Bates [EMAIL PROTECTED] wrote: Do you mean that you want to collapse similar rows into a single row and perhaps a count of the number of times that this row occurs? Let me rephrase the problem by providing an example. Input: A = [,1] [,2] [1,]11

Re: [R] subset rows in two dataframes

Dear Jim, The following codes maybe helps. for (i in 1:length(dat1[,1])) { for (j in 1:length(dat2[,1])) { if (dat1[i,1] == dat2[j,1] dat1[i,2] == dat2[j,2]) print (j) } } time1- as.Date(c(2006-01-03, 2006-05-03, 2006-05-04, 2006-05-11, 2006-05-12, 2006-05-16, 2006-05-19,

[R] Robust Stepwise Regression

Hello, I am interested in performing a stepwise regression using the robust regression technique to estimate the models at each stage. At the moment I am using the lm code and the stepAIC code to select the best OLS model. This final model is then re-estimated using the rlm code. I am hoping

Re: [R] lme nesting/interaction advice

On 10 May 2008, at 07:36, Kingsford Jones wrote: Federico, I think you'll be more likely to receive the type of response you're looking for if you formulate your question more clearly. The inclusion of commented, minimal, self-contained, reproducible code (as is requested at the bottom of

[R] Fundamental formula and dataframe question.

There is a very useful and apparently fundamental feature of R (or of the package pls) which I don't understand. For datasets with many independent (X) variables such as chemometric datasets there is a convenient formula and dataframe construction that allows one to access the entire X matrix

[R] positioning of color key in levelplot

Is there a way of positioning the color key in levelplot when the axes are on a categorical (rather than numerical) scale? I've put some sample code below. I need to add a secondary y axis to the right side of my plot but then the labels interfere with the color key (which is currently on the

Re: [R] Fundamental formula and dataframe question.

On 11-May-08 18:58:45, Myers, Brent wrote: There is a very useful and apparently fundamental feature of R (or of the package pls) which I don't understand. For datasets with many independent (X) variables such as chemometric datasets there is a convenient formula and dataframe construction

Re: [R] positioning of color key in levelplot

On 5/11/08, E C [EMAIL PROTECTED] wrote: Is there a way of positioning the color key in levelplot when the axes are on a categorical (rather than numerical) scale? I've put some sample code below. I need to add a secondary y axis to the right side of my plot but then the labels interfere

Re: [R] howto import .xls and .ods

On Fri, May 02, 2008 at 07:35:37AM +0100, Prof Brian Ripley wrote: There is a *manual* on R Data Import/Export, not just an FAQ. This is the first request I have seen for .ods (whatever that is -- The most well-known application that uses this file format is the Calc (Spreadsheet) part of

[R] Finding unmatched data between two dataframe using several factors

Hi R users I am trying to find unmatched data from two dataframes. I would like to find unmatched data based on several factors. For the following data: dat1 - data.frame(x = paste(A, 1:6, sep=), y = c(andy,bob,ciaran,dan, eion, fred)) dat1 dat2

Re: [R] lme nesting/interaction advice

On Sun, May 11, 2008 at 07:52:50PM +0100, Federico Calboli wrote: The main point of my question is, having a 3 way anova (or ancova, if you prefer), with *no* nesting, 2 fixed effects and 1 random effect, why is it so boneheaded difficult to specify a bog standard fully crossed model?

Re: [R] Not to draw the xaxis ticks in ggplot2

Hi ronggui wrote: Thanks,Hadley. Another question is how can I know there is a grob path called xaxis::ticks? Is there any easy way to figure out? grid.ls() ? Paul Thanks in advance. On Wed, May 7, 2008 at 10:06 PM, hadley wickham [EMAIL PROTECTED] wrote: On Tue, May 6, 2008 at

Re: [R] lme nesting/interaction advice

On 12/05/2008, at 9:45 AM, Andrew Robinson wrote: On Sun, May 11, 2008 at 07:52:50PM +0100, Federico Calboli wrote: The main point of my question is, having a 3 way anova (or ancova, if you prefer), with *no* nesting, 2 fixed effects and 1 random effect, why is it so boneheaded difficult to

Re: [R] Finding unmatched data between two dataframe using severalfactors

Andrew Have you tried: dat2[!paste(dat2[,1], dat2[,2])%in%paste(dat1[,1], dat1[,2]),] ? HTH Peter Alspach -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Andrew McFadden Sent: Monday, 12 May 2008 9:40 a.m. To: r-help@r-project.org

[R] Running external python script

Hi R users, I have a python script. Assume that I would like to run this external python script when executing a R script or command line R. I don't know how to this using R code. I really appreciate if you could help me. Thanks for your help Rostam [[alternative HTML version deleted]]

Re: [R] Running external python script

On 12/05/2008, at 11:17 AM, rostam shahname wrote: Hi R users, I have a python script. Assume that I would like to run this external python script when executing a R script or command line R. I don't know how to this using R code. I really appreciate if you could help me. Thanks for your

Re: [R] Running external python script

mac os X On Sun, May 11, 2008 at 8:24 PM, Rolf Turner [EMAIL PROTECTED] wrote: On 12/05/2008, at 11:17 AM, rostam shahname wrote: Hi R users, I have a python script. Assume that I would like to run this external python script when executing a R script or command line R. I don't know

Re: [R] Running external python script

On 12/05/2008, at 11:37 AM, rostam shahname wrote: mac os X No, no, no!!! That's what ***you*** type (in R)! :-) cheers, Rolf On Sun, May 11, 2008 at 8:24 PM, Rolf Turner [EMAIL PROTECTED] wrote: On 12/05/2008, at 11:17 AM, rostam

Re: [R] lme nesting/interaction advice

On Mon, May 12, 2008 at 10:34:40AM +1200, Rolf Turner wrote: On 12/05/2008, at 9:45 AM, Andrew Robinson wrote: On Sun, May 11, 2008 at 07:52:50PM +0100, Federico Calboli wrote: The main point of my question is, having a 3 way anova (or ancova, if you prefer), with *no* nesting, 2 fixed

Re: [R] [Partial Summary] A fundamental formula and dataframe question.

Two very good responses to this question, but I wonder, Is there some more complete documentation on using this form of model and dataframe construction? I've been using R for ~5 years now and wasn't aware of it. Response 1: Insert a matrix as a column of the dataframe using I().

[R] How to eliminate perticular date

Hi R-expert, I try to eliminate 29 Feb but I got an error message below: feb_data1 - Pooraka_data[Pooraka_data$Month==2,] feb_28days - feb_data1 [feb_data1$Day==28,] feb_29days - feb_data1 [feb_data1$Day==29,] ## delete 29 Feb feb_no_29 - feb_data1 [-(feb_29days),] feb_no_29 - feb_data1

Re: [R] How to eliminate perticular date

What you are getting back in 'feb_no_29' is a data frame which can not be used for indexing. Do 'str(feb_no_29)' and see what you set. What you probably want to do is to use 'which' to find out which rows match the criteria and then delete them: feb_29 - which((Pooraka_data$Month == 2)

[R] problem configuring package udunits

Hi R Users, I am new to running R on a Linux platform (I'm used to Windows) - I'm running R 2.7.0 on Ubuntu 7.10 (Gutsy) as sudo (without Emacs). My architecture is Pentium D (x86_64). I am having problems successsfully configuring the downloaded package 'udunits'. When I execute

Re: [R] Not to draw the xaxis ticks in ggplot2

Thanks. It is helpful. Best On Mon, May 12, 2008 at 4:34 AM, Paul Murrell [EMAIL PROTECTED] wrote: Hi ronggui wrote: Thanks,Hadley. Another question is how can I know there is a grob path called xaxis::ticks? Is there any easy way to figure out? grid.ls() ? Paul Thanks in

Re: [R] Creating data.frames dynamically

Hi Worik: Is this what you want? Names - letters[1:5]; Dates- 1:20; d- data.frame(dates=Dates, a = vector(mode=numeric, length=length(Dates))); for(i in 2:5){d[,paste(sep=,Names[i])]-with(d,a)}; d dates a b c d e 1 1 0 0 0 0 0 2 2 0 0 0 0 0 3 3 0 0 0 0 0 4 4 0 0 0 0

Re: [R] How to eliminate perticular date

Why dont you use the rattle GUI for data selection ? Regards, Ajay On Mon, May 12, 2008 at 6:24 AM, Roslina Zakaria [EMAIL PROTECTED] wrote: Hi R-expert, I try to eliminate 29 Feb but I got an error message below: feb_data1 - Pooraka_data[Pooraka_data$Month==2,] feb_28days - feb_data1

[R] Insert a recorde into a table using SQL

Dear list, I want to insert a recorde into a SQLite table by dbGetQuery(), but there is a problem when the value contains quotation mark. dd-data.frame(txt=c(having both ' and \ in character.,OK)) library(RSQLite) Loading required package: DBI con-dbConnect(dbDriver(SQLite),:memory:)