Re: [R] How to add a new row in to an existing data set in R Language?

2012-10-04 Thread PIKAL Petr
Hi ?merge or ?rbind Regards Petr > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > project.org] On Behalf Of killerkarthick > Sent: Friday, October 05, 2012 6:39 AM > To: r-help@r-project.org > Subject: [R] How to add a new row in to an existing data s

Re: [R] problem with convergence in mle2/optim function

2012-10-04 Thread Bert Gunter
Presumably you checked out the CRAN Optimization task view to see if you could find any joy there, right? (I doubt there is, but ...) -- Bert On Thu, Oct 4, 2012 at 10:12 PM, Adam Zeilinger wrote: > Hello R Help, > > I am trying solve an MLE convergence problem: I would like to estimate four > p

Re: [R] Rainflow, range pair counting

2012-10-04 Thread Dandy Ehlert
That is right. It is a part of materials engineering. Could say it’s a standard method in it. Better description is always on Wikipedia :) http://en.wikipedia.org/wiki/Rainflow-counting_algorithm Have found packs for all other big statistic software like Matlab or Diadem but none for R. And I

Re: [R] How to add a new row in to an existing data set in R Language?

2012-10-04 Thread killerkarthick
Thanks for your answer Arun, That above image shows the merge in SAS but its possible as same in R please reply me let me know? -- View this message in context: http://r.789695.n4.nabble.com/How-to-merge-data-set-in-R-Language-tp4644855p4645119.html Sent from the R help mailing list archive at

[R] How to add a new row in to an existing data set in R Language?

2012-10-04 Thread killerkarthick
Hi. I have one dataset with 3 variables with 4 observations. I want to add the some more observation into the existing dataset. Data;Mydata center drugchange 101 P 18 102 P 14 103 D 12 104 D 18 105 P 10 This is the data set. So

Re: [R] How to add a new row in to an existing data set in R Language?

2012-10-04 Thread killerkarthick
Hi Ista Zahn, Thanks for your advice. Please see the following Image. i am expecting the result should same in the image. -- View this message in context: http://r.78969

[R] problem with convergence in mle2/optim function

2012-10-04 Thread Adam Zeilinger
Hello R Help, I am trying solve an MLE convergence problem: I would like to estimate four parameters, p1, p2, mu1, mu2, which relate to the probabilities, P1, P2, P3, of a multinomial (trinomial) distribution. I am using the mle2() function and feeding it a time series dataset composed of fou

Re: [R] Subsetting a group of data

2012-10-04 Thread arun
HI, Another way to get the result: dat1 <- read.table(text=" Cola   Colb 1  1  1 2  1  1 3  1  1 4  -1    1 5  -1    -1 6  -1    -1 7  1  -

Re: [R] barplot with some 0 frequencies

2012-10-04 Thread Richard M. Heiberger
barchart(-3:3 ~ letters[1:7], origin=0) On Fri, Oct 5, 2012 at 12:05 AM, David Winsemius wrote: > > On Oct 4, 2012, at 4:49 PM, Guillaume2883 wrote: > > > Hi all, > > > > I am back with a new question ! > > I recorded the occurence of 4 differents event on 20 places for a given > time > > period

Re: [R] barplot with some 0 frequencies

2012-10-04 Thread David Winsemius
On Oct 4, 2012, at 4:49 PM, Guillaume2883 wrote: > Hi all, > > I am back with a new question ! > I recorded the occurence of 4 differents event on 20 places for a given time > period. > Now, I want to do some barplot of the frequency of theses events for each > place, so it should be easy. My pr

Re: [R] Text to Speech In R

2012-10-04 Thread C.H.
Pizza can be ordered from R (fortune 66), why not tts? Under Debian, I can tts by apt-get install the tts engine espeak and then use the simple system call under R > system("espeak", input="hello world, howdy") On Thu, Oct 4, 2012 at 8:31 PM, Alex Zhang wrote: > Dear all, > > Excuse me for my

Re: [R] Error not found - for file retrieval with head()

2012-10-04 Thread jim holtman
Couple of things: 1. 'head' it to look at objects, not files 2. you specified 'turtlehatch.csv' without quotes which told R to look for an object by that name 3. you probably want to use 'file.show' On Thu, Oct 4, 2012 at 5:47 PM, Jhope wrote: > Hi R-listers, > > I am in plyr and using the scrip

Re: [R] Calculating all possible ratios

2012-10-04 Thread genome1976
Hi Rui, A while ago you helped me with calculaing all possible ratios from a dataset. This is the code I am using as suggested by you. data <- read.table("new_data.txt", header=T, row.names=1, sep="\t") pairwise.ratios <- function(x, prefix="probeset", char=":"){ n <- ncol(x)

[R] barplot with some 0 frequencies

2012-10-04 Thread Guillaume2883
Hi all, I am back with a new question ! I recorded the occurence of 4 differents event on 20 places for a given time period. Now, I want to do some barplot of the frequency of theses events for each place, so it should be easy. My problem is that I want to see the frequencies of the 4 events on my

[R] Error not found - for file retrieval with head()

2012-10-04 Thread Jhope
Hi R-listers, I am in plyr and using the script below and when I try to pull up the .csv file with head() it does not show up. I am not sure where I am going wrong. I have also checked the file box for turtlehatch.csv and it still does not show up when I head(turtlehatch.csv). Am I missing a comma

Re: [R] Subsetting a group of data

2012-10-04 Thread arun
Hi, You can also try this: dat1 <- read.table(text=" Cola   Colb 1  1  1 2  1  1 3  1  1 4  -1    1 5  -1    -1 6  -1    -1 7  1  -1 8   

[R] converting matlab code including C/mex to R

2012-10-04 Thread Kjetil Halvorsen
The title says it all! I am starting upon such a project, but now next to nothing about matlab/mex and how matlab-to-C communication is done. Any hints? links? (I am starting now to read up on mex) kjetil __ R-help@r-project.org mailing list https://st

Re: [R] [R-SIG-Mac] problem with the installation of r commander on a mac

2012-10-04 Thread Simon Urbanek
In fun(libname, pkgname) : couldn't connect to display ":0" suggests that X11 is not running. Note that X11 is no longer part of OS X since Mountain Lion. Cheers, Simon On Oct 4, 2012, at 4:32 PM, John Fox wrote: > Dear Marco, > > To provide some background, you and I have corresponded about

Re: [R] smoothScatter plot

2012-10-04 Thread JiangZhengyu
Hi John, Thanks for your email. Your way works good. However, I was wondering if you can help with a smoothed scatter plot that has shadows with different darker blue color representing higher density of points. Zhengyu Date: Thu, 4 Oct 2012 05:46:46 -0800 From: jrkrid...@inbox.com Subject: R

Re: [R] convert multi dimensional array to list

2012-10-04 Thread David Winsemius
On Oct 4, 2012, at 8:57 AM, anto.r wrote: > Hi Michael > > thanks! That was the option if I kept it an array. The list format with $ > sign since it leaves me feeling that the names are there and can be easily > accessed. Why would you rather not use the $ sign? It would be better in programing

Re: [R] help with making figures

2012-10-04 Thread Bert Gunter
First convert your data from short fat to long thin format: Column headings for the transformed data frame are: Year Site Result This can be done in a variety of ways, but check out the reshape2 package or the reshape function ins the stats package. Once you have the data in this form,

[R] extract fit values from geom_smooth

2012-10-04 Thread Joseph Tannhuber
Dear all, I have a plot with two gaussian distributions: ggplot( tofdf, aes( x = tof, y = counts ) ) + geom_histogram( stat = "identity", position = "stack", fill = 2 ) + geom_smooth( method = "nls", formula = y ~ (a/b)*exp(-(x-c)^2/(2*b^2)) + (d/e)*exp(-(x-f)^2/(2*e^2)), se=F,

[R] help with making figures

2012-10-04 Thread megalops
I need to make about 30 figures and I am trying to create a program in R that will make my life a lot easier. First I will tell you how my data is setup. I have 30 sites and then data for each year at the site. I have 10 years of data for each site. Below is a small chunk of my data to show th

Re: [R] Creating vegetation distance groups from one column

2012-10-04 Thread Jean V Adams
Try this data.to.analyze$VegIndex <- cut(data.to.analyze$Veg, breaks=seq(0, 70, 5), include.lowest=TRUE) plot(data.to.analyze$VegIndex) Jhope wrote on 10/04/2012 02:25:09 PM: > > Hi, > > Allow me to recap my question. In plyr I am trying to group distances of > nests to the vegetati

Re: [R] Coda, HPDinterval and multiple chains

2012-10-04 Thread ilai
How about require(coda) data(line) str(line) cline <- as.mcmc(do.call(rbind,line)) str(cline) # Thus HPDinterval(cline) # (or any FUN.mcmc) sum(cline[,'alpha'] >= median(cline[,'alpha']))/nrow(cline) Cheers On Thu, Oct 4, 2012 at 12:22 PM, NORRIS Paul wrote: > Dear all, > > I'm not 100% sure

Re: [R] R help - Adding a column in a data frame with multiple conditions

2012-10-04 Thread arun
Hi, By extending Sarah's solution to the whole dataset: dat1<-read.table(text="  A B C  D  1 X 90 88  1 Y 72 70  1 Z 67 41  2 X 74 49  2 Y 42 50  2 Z 81 56  3 X 92 59  3 Y 94 80  3 Z 80 82 ",sep="",header=TRUE,stringsAsFactors=FALSE) dat1$E<-unlist(lapply(split(dat1,dat1$A),function(x) ifelse(x$

Re: [R] problem with the installation of r commander on a mac

2012-10-04 Thread John Fox
Dear Marco, To provide some background, you and I have corresponded about your problem, and now at least you appear to have the tcltk package installed. Did you reinstall R as I suggested? Have you checked that tcltk is working, following the instructions I sent you in an earlier email? I don't

Re: [R] can stepAIC be customized to exclude coefficients with p-value less than certain values?

2012-10-04 Thread Marc Schwartz
On Oct 4, 2012, at 3:22 PM, David Winsemius wrote: > > On Oct 4, 2012, at 11:26 AM, liang@us.pwc.com wrote: > >> For example, if coefficient's p-value is less than 0.1 I want the stepwise >> to automatically drop that variable. Can the stepAIC be customized to do >> that? SAS seems to b

Re: [R] can stepAIC be customized to exclude coefficients with p-value less than certain values?

2012-10-04 Thread David Winsemius
On Oct 4, 2012, at 11:26 AM, liang@us.pwc.com wrote: > For example, if coefficient's p-value is less than 0.1 I want the stepwise > to automatically drop that variable. Can the stepAIC be customized to do > that? SAS seems to be able to customized stepwise function with p-value > or cook

Re: [R] R help - Adding a column in a data frame with multiple conditions

2012-10-04 Thread Sarah Goslee
Hi Libby, You had an accumulation of small errors, from an extra ) to an unclear understanding of how indexing works in R. Also, you shouldn't call your dataframe df, or use square brackets in column names. That said, what about: sampledata <- structure(list(A = c(1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L

Re: [R] (no subject)

2012-10-04 Thread R. Michael Weylandt
On Thu, Oct 4, 2012 at 8:55 PM, Sarah Goslee wrote: > On Thu, Oct 4, 2012 at 1:40 PM, Chia-Yi Chiu > wrote: >> Hi >> >> I would like to learn how the R function "hclust" deals with ties. It is >> written in Fortran, so I cannot access the code. Thanks!! > > Sure you can. Just download the source

[R] geoRglm with factor variable as covariable

2012-10-04 Thread Filoche
Dear R users. I'm trying to fit a generalised linear spatial mode using the geoRglm package. To do so, I'm preparing my data (geodata) as follow: geoData9093 = as.geodata(data9093, coords.col= 17:18, data.col=15,* covar.col=16*) where covar.col is a factor variable (years in this case 90-91-92-

Re: [R] (no subject)

2012-10-04 Thread Jeff Newmiller
Please use a relevant subject line. You can easily download source code for all CRAN packages and read it yourself. --- Jeff NewmillerThe . . Go Live... DCN:Basics: ##.#.

Re: [R] Creating vegetation distance groups from one column

2012-10-04 Thread Jhope
Hi, Allow me to recap my question. In plyr I am trying to group distances of nests to the vegetation into classes that are defined by (0-5m, 6-10m, 11-15m, 16-20m, 21-25m, 26-30m, 31-35m). Each row is a nest and all the distances to the vegetation is in a column. As well I would like to create c

Re: [R] (no subject)

2012-10-04 Thread Sarah Goslee
On Thu, Oct 4, 2012 at 1:40 PM, Chia-Yi Chiu wrote: > Hi > > I would like to learn how the R function "hclust" deals with ties. It is > written in Fortran, so I cannot access the code. Thanks!! Sure you can. Just download the source from r-project.org. -- Sarah Goslee http://www.functionaldiv

[R] R help - Adding a column in a data frame with multiple conditions

2012-10-04 Thread Libby M Gertken
Hi, I am trying to add a column of numbers to a data frame in R with multiple conditions. Here is a simplified example df: [A] [B] [C] [D] [E] [1] 1 X 90 88 [2] 1 Y 72 70 [3] 1 Z 67 41 [4] 2 X 74 49 [5] 2 Y 42 50 [6] 2 Z 81 56 [7] 3 X 92 59 [8] 3 Y 94 80 [9] 3 Z 80 82 I would like column [E] t

[R] P value discrepancy plots

2012-10-04 Thread bilelsan
Hello R-users, I am trying to plot P value discrepancy plots as describe by Russell Davidson and James G. MacKinnon in "Graphical Methods for Investigating the Size and Power of Hypothesis Tests (1998)". I know ecdf () to obtain the empirical cumulative distribution frequencies and knots() to extra

[R] Coda, HPDinterval and multiple chains

2012-10-04 Thread NORRIS Paul
Dear all, I'm not 100% sure if this question is best directed at the r-list, or a mailing list concerned with Bayesian analysis, so please accept my apologies if another audience may be more appropriate. I have been using the rjags package to run Jags models with multiple chains and store the

[R] data structure for plsr

2012-10-04 Thread Emma Jones
I am having a similar problem understanding the data structure of the "yarn" dataset described in the "[R] data structure for plsr" posts. I have spectroscopic data I'd like to run through a PLSR and have read the tutorial series, but still do not understand the data format required for the code to

[R] (no subject)

2012-10-04 Thread Chia-Yi Chiu
Hi I would like to learn how the R function "hclust" deals with ties. It is written in Fortran, so I cannot access the code. Thanks!! [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinf

Re: [R] R combining vectors into a data frame but without a continuous common variable

2012-10-04 Thread lucy88
Oh my word, you're a genius!! That is absolutely perfect, thank you so much!! I've no idea how you've learnt these things but I would never ever have been able to do that. You've just made my day so much better after the horror of confusion before. Thank you!! -- View this message in conte

[R] problem with the installation of r commander on a mac

2012-10-04 Thread Marco Mello
Dear list members, I’m trying to install R Commander under Mac OSX Mountain Lion (10.8.2). After following all the steps described in the installation notes (http://socserv.socsci.mcmaster.ca/jfox/Misc/Rcmdr/installation-notes.html), I got this error message: = Loading required package: tc

Re: [R] Broken Links on http://www.r-project.org

2012-10-04 Thread Hasan Diwan
The R Graphics Gallery has moved to http://gallery.r-enthusiasts.com/ and there's another R Graphics Manual at http://rgm2.lab.nig.ac.jp/RGM2 -- H On 26 September 2012 04:56, Viechtbauer Wolfgang (STAT) < wolfgang.viechtba...@maastrichtuniversity.nl> wrote: > I was not sure who I should contact a

Re: [R] R combining vectors into a data frame but without a continuous common variable

2012-10-04 Thread arun
Hi Lucy, No problem. Just a correction to my earlier email. dat1<-read.table("Landeck_vec.txt",sep="",header=TRUE,stringsAsFactors=FALSE) dat2<-read.table("Kaurnetal_vec.txt",sep="",header=TRUE,stringsAsFactors=FALSE) colnames(dat1)[1]<-"Date" (Rui: #dat2 Date format is inconsistent.) dat2$Date

Re: [R] Broken Links on http://www.r-project.org

2012-10-04 Thread Uwe Ligges
On 26.09.2012 10:56, Viechtbauer Wolfgang (STAT) wrote: I was not sure who I should contact about this, so I am posting this here. There are a few broken links on the R website. 1) http://www.r-project.org/search.html -> link to the Nabble R Forum. I belive the correct/new URL should be: htt

[R] can stepAIC be customized to exclude coefficients with p-value less than certain values?

2012-10-04 Thread liang . che
For example, if coefficient's p-value is less than 0.1 I want the stepwise to automatically drop that variable. Can the stepAIC be customized to do that? SAS seems to be able to customized stepwise function with p-value or cooks'd. thanks!

Re: [R] Failure of sas7bdat package

2012-10-04 Thread Nordlund, Dan (DSHS/RDA)
> -Original Message- > From: b.rowling...@gmail.com [mailto:b.rowling...@gmail.com] On Behalf > Of Barry Rowlingson > Sent: Thursday, October 04, 2012 10:12 AM > To: Nordlund, Dan (DSHS/RDA) > Cc: Fisher Dennis; r-h...@stat.math.ethz.ch > Subject: Re: [R] Failure of sas7bdat package > > On

Re: [R] convert multi dimensional array to list

2012-10-04 Thread arun
Hi, You can try this: list2<-lapply(1:dim(my.df.cast)[3],function(i) my.df.cast[,,i,])  names(list2)<-unlist(dimnames(my.df.cast)[3])  list2$May #, , variable = tmax # #  name #year   Alf Que  # 2011  18  16  # 2012  18  16 #, , variable = tmin #  #     name #year   Alf Que  # 2011  12  11  #

Re: [R] Rainflow, range pair counting

2012-10-04 Thread Bert Gunter
I'm going to make a wild guess that this might have something to do with "bump hunting" =local conditional extrema. If so, the prim package might be applicable. Otherwise ("...cycles in time series") time series methods might be appropriate. But I basically agree with others that the OP has not c

Re: [R] Error in if (any(ch)) { : missing value where TRUE/FALSE needed

2012-10-04 Thread Uwe Ligges
On 03.10.2012 19:47, liang@us.pwc.com wrote: Can someone please help with the error message below? Not really, the message says it all: You cannot improve the AIC since it is numerically -Inf already! Uwe Ligges thanks! Start: AIC=-Inf value ~ 1 + Core_CPI__ + GDP_change + Un

Re: [R] Errors when saving output from WinBUGS to R

2012-10-04 Thread Uwe Ligges
So have you closed WinBUGS regularly (and not via Esc in R) now? Uwe Ligges On 03.10.2012 21:29, yiyiwang wrote: Sorry for any cross messaging. Here are what I've tried and more info. The error sticks :( I installed the latest R2.15.1 (64bit) on my machine, per Uwe's comments. The same error

Re: [R] R combining vectors into a data frame but without a continuous common variable

2012-10-04 Thread Rui Barradas
Hello, Try the following. url1 <- "http://r.789695.n4.nabble.com/file/n4644986/Landeck_vec.txt"; url2 <- "http://r.789695.n4.nabble.com/file/n4644986/Kaurnetal_vec.txt"; dat1 <- read.table(url1, header = TRUE) dat2 <- read.table(url2, header = TRUE) str(dat1) str(dat2) # Precip is a factor,

Re: [R] Rainflow, range pair counting

2012-10-04 Thread Marc Schwartz
A quick Google search suggests that the methods are used in material fatigue studies where varying stress loads are placed upon the material, which puts this into a materials engineering context. Using rseek.org does not give any joy. Perhaps someone with experience in that domain will chime in

Re: [R] Legend Truncated Using filled.contour

2012-10-04 Thread David Winsemius
On Oct 3, 2012, at 12:58 PM, Kirsten wrote: > Hey everyone, > > I'm working on a contour plot depicting asymptomatic prevalence at varying > durations of infectiousness and force of infection. I've been able to work > everything out except for this one - my legend title keeps getting cut off.

Re: [R] Failure of sas7bdat package

2012-10-04 Thread Barry Rowlingson
On Thu, Oct 4, 2012 at 6:07 PM, Nordlund, Dan (DSHS/RDA) wrote: > Looking at the source code on github, Matt has already added W32_7PRO to the > list of known hosts. So, you could download the source tarball and build the > package yourself, or just extract the SAS7bdat.R file from the tarball

Re: [R] Failure of sas7bdat package

2012-10-04 Thread Nordlund, Dan (DSHS/RDA)
> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > project.org] On Behalf Of Barry Rowlingson > Sent: Thursday, October 04, 2012 9:52 AM > To: Fisher Dennis > Cc: r-h...@stat.math.ethz.ch > Subject: Re: [R] Failure of sas7bdat package > > On Thu, Oct 4, 2

Re: [R] Failure of sas7bdat package

2012-10-04 Thread Barry Rowlingson
On Thu, Oct 4, 2012 at 4:55 PM, Fisher Dennis wrote: > R 2.15.1 > OS 10.7 > > Colleagues > > I have been an enthusiastic user of the sas7bdat package in R. However, > several recent sas7bdat files sent to me from different sources cannot be > read by the package. The error message is: >

Re: [R] "Explore" SPSS function in R

2012-10-04 Thread Ista Zahn
Hi Maria, I suggest starting with the by() function. Here is an example to get you started: by(mtcars$mpg, INDICES = mtcars[c("cyl", "am")], FUN = function(x) { list(summary = summary(x), mean.CI95 = confint(lm(x ~ 1))), }) Best, Ista On Thu, Oct 4, 2012 at 10:28 AM,

Re: [R] Find the average line from several lines

2012-10-04 Thread David Winsemius
On Oct 3, 2012, at 10:31 PM, Greg Lollback wrote: > Hi people, > > I'm looking at the influence of sample size on a dataset, so I randomly > sampled from the data with different sample sizes multiple times using > the "sample" and "replicate" function. For example: > > /replicate(1000, sample

Re: [R] convert multi dimensional array to list

2012-10-04 Thread anto.r
Hi Michael thanks! That was the option if I kept it an array. The list format with $ sign since it leaves me feeling that the names are there and can be easily accessed. Why would you rather not use the $ sign? I use R-Studio and there names can be selected from a drop-down list, I have found it

[R] Analysis of Means in R

2012-10-04 Thread Data Analytics Corp.
Hi, Does anyone know of a package for doing Analysis of Means in R? Or have any functions they would like to share? Thanks, Walt Walter R. Paczkowski, Ph.D. Data Analytics Corp. 44 Hamilton Lane Plainsboro, NJ 08536 (V) 609-936-8999 (F) 609-

Re: [R] "Explore" SPSS function in R

2012-10-04 Thread William Dunlap
Can you tell us what "Explore" gives you for the following dataset? > d<-data.frame(Quant=round(log2(1:24), 2), Categ=LETTERS[rep(1:4,c(8,7,6,3))]) > write.table(d, row.names=FALSE, sep="\t") "Quant" "Categ" 0 "A" 1 "A" 1.58"A" 2 "A" 2.32"A" 2.58"A" 2.81"A" 3

Re: [R] Failure of sas7bdat package

2012-10-04 Thread Anthony Damico
I can't help troubleshoot the read.sas7bdat function, but if you have the fixed-width data and a SAS importation script, the SAScii package might help you work around this issue. :) On Thu, Oct 4, 2012 at 11:55 AM, Fisher Dennis wrote: > R 2.15.1 > OS 10.7 > > Colleagues > > I have been an ent

Re: [R] Subsetting a group of data

2012-10-04 Thread Rui Barradas
Hello, Try the following. dat <- read.table(text=" Cola Colb 1 1 1 2 1 1 3 1 1 4 -11 5 -1-1 6 -1-1 7 1 -1 8

[R] Failure of sas7bdat package

2012-10-04 Thread Fisher Dennis
R 2.15.1 OS 10.7 Colleagues I have been an enthusiastic user of the sas7bdat package in R. However, several recent sas7bdat files sent to me from different sources cannot be read by the package. The error message is: Error in read.sas7bdat(FILENAME) : unknown host W32_7PRO please r

Re: [R] Rainflow, range pair counting

2012-10-04 Thread R. Michael Weylandt
R is used by many many folks, and most of us don't have the domain knowledge to make heads or tails of what you're mentioning here. You'll need to greatly clarify and perhaps ask on the appropriate R-SIG-* list. (Ecology maybe?) Perhaps rseek.org is also of help. Cheers, Michael On Thu, Oct 4, 20

Re: [R] "Explore" SPSS function in R

2012-10-04 Thread R. Michael Weylandt
On Thu, Oct 4, 2012 at 3:28 PM, MARIA RODRIGUEZ wrote: > Hi everyone, > > Does anybody knows if there is an equivalent R function that gives the same > outcome as in "Explore" function in SPSS ? > (Analize->Descriptive Statistics->Explore) > > It does a categorical vs quantitative variables analy

[R] problems with plotting and printing aareg

2012-10-04 Thread Asa Johannesen
Hi all, I've ventured into the world of nonparametric survival and I would like to use the "maxtime" option for printing and plotting my aareg fit. However, my fit does not have "test.var2" and this stops the print and plot when adding a maxtime. My code is as follows: Response<-Surv(Time,Eve

Re: [R] Help with R Fitting an inverse Gamma

2012-10-04 Thread R. Michael Weylandt
On Thu, Oct 4, 2012 at 11:42 AM, kmammasis wrote: > Dear all, > > I am new in R and would like to ask for someone's help in understanding > where I go wrong with the following code: > # > rm(list=ls()) # Please don't include this line into your postings: it's quite rude to clobber the data of fol

Re: [R] running a regression

2012-10-04 Thread R. Michael Weylandt
You'll need to supply a reproducible example: http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example Cheers, Michael On Thu, Oct 4, 2012 at 2:19 PM, Dereje Bacha wrote: > Hi all, > > When I run a regression, I got the following response. your help is > appreciat

[R] "Explore" SPSS function in R

2012-10-04 Thread MARIA RODRIGUEZ
Hi everyone, Does anybody knows if there is an equivalent R function that gives the same outcome as in "Explore" function in SPSS ? (Analize->Descriptive Statistics->Explore) It does a categorical vs quantitative variables analysis. ( But not linear regression) I need to compare intragroup (cat

Re: [R] Subsetting a group of data

2012-10-04 Thread R. Michael Weylandt
On Thu, Oct 4, 2012 at 3:20 PM, Ian Arvin wrote: > Hi, > > > I am making my way down the learning curve of R, and find it a great > language with so many helpful users! Up it hopefully, no? > > Below is an example of what I'm trying to do, but can't quite figure out the > right path to go down.

Re: [R] Problem with colors in contour plot

2012-10-04 Thread Loukia Spineli
The contpour plot is very impressive!!! The contour plot I want to display should have only 2 colors: white a grey. I have created this plot for only one trial (see, ppts) and now I am attempting to do the same for many trials. On Thu, Oct 4, 2012 at 3:51 PM, Rui Barradas wrote: > Hello, > > So

Re: [R] How to convert array dim names to expressions for plotting?

2012-10-04 Thread Hofert Jan Marius
Dear Bert, thanks, I was confused (of course). So the right question would have been: How to convert strings to expressions (plotmath expressions). But that I could easily find, thanks. Cheers, Marius Here is the working example (showing plotmath and just-character labels): require(grid) ta

[R] R combining vectors into a data frame but without a continuous common variable

2012-10-04 Thread lucy88
Hello, I have two different files which I'd like to combine to make one data frame but I've no idea how to do it! The first file has two columns; one is the date, the following is a binary code for debris flow events. Then my other file has also two columns; the date and then precipitation data.

Re: [R] Error in installing package wordcloud in R

2012-10-04 Thread R. Michael Weylandt
Try using a different repo: something seems off about whichever one you've tried. M On Thu, Oct 4, 2012 at 9:33 AM, M Bindu/TVM/TCS wrote: > Dear Sir, > > I faced problem during the installation of wordcloud package in R -2.15 . I > have mentioned the error below. After selecting the cran mirro

[R] RPostgreSQL Date format problem

2012-10-04 Thread scottyjock
Hi, I've just started using RPostgreSQL within R to interrogate a PostgreSQL database. The problem I have is that a column containing dates shows up differently in R than in the database: here is the output direct from the psql database in Linux: db=> select pat_dob_dt from patients where pat_dob

[R] Subsetting a group of data

2012-10-04 Thread Ian Arvin
Hi, I am making my way down the learning curve of R, and find it a great language with so many helpful users! Below is an example of what I'm trying to do, but can't quite figure out the right path to go down. Here's what I have: Main is a time series of data with columns Cola and Colb

Re: [R] convert multi dimensional array to list

2012-10-04 Thread R. Michael Weylandt
On Thu, Oct 4, 2012 at 2:54 PM, anto.r wrote: > I have an 4D named array (the result of recast() on a dataframe) and would > like to convert it to a list with the names retained. > > Example: >>my.df<-expand.grid(name=c("Alf","Que"), month=c("May","Jun"), > year=c("2011","2012")) >>my.df$tmax<-sam

[R] Error in installing package wordcloud in R

2012-10-04 Thread M Bindu/TVM/TCS
Dear Sir, I faced problem during the installation of wordcloud package in R -2.15 . I have mentioned the error below. After selecting the cran mirror it shows the error. I will be thankful for your response to this solution. Error Selection: 34 Warning: unable to access index for repository

[R] (no subject)

2012-10-04 Thread dick liang
producing a multi-figure plot, i am try to add beta_1, beta_2,.. beta_9 to ylab using expression or substitution, but cannot work out like for (i in 1:9){ plot(seq(1/m, 1-1/m, 1/m), beta.q[,i], type="l", col=1, ylim=range(beta.q), xlab="quantile", ylab=expresion(beta[i])) } any suggestio

[R] Rainflow, range pair counting

2012-10-04 Thread Dandy Ehlert
Hello I got some question about R. Im searching for a function or R-code which makes me some collectives. In the literature they are called: > level crossing counting > rainflow counting > range pair counting Cant believe that’s not available for R. May I just loo

[R] running a regression

2012-10-04 Thread Dereje Bacha
Hi all, When I run a regression, I got the following response.  your help is appreciated.  Error in constants (lhs, cnames_symb) :   The hypothesis "Cost_PKL-QK.X_(Intercept)=0" is not well formed: contains bad coefficient/variable names.  In addition: Warning message: In constants(lhs, cnames

[R] convert multi dimensional array to list

2012-10-04 Thread anto.r
I have an 4D named array (the result of recast() on a dataframe) and would like to convert it to a list with the names retained. Example: >my.df<-expand.grid(name=c("Alf","Que"), month=c("May","Jun"), year=c("2011","2012")) >my.df$tmax<-sample(15:20,4) >my.df$tmin<-sample(10:15,4) >my.df.melt=mel

Re: [R] t-test

2012-10-04 Thread Jose Iparraguirre
Hi, You can also have a look at this paper (no subscription needed): Herberich, E.; Sikorski, J.; and Hothorn, T. (2010). "A Robust Procedure for Comparing Multiple Means under Heteroscedasticity in Unbalanced Designs". PLoSONE, Vol. 5, Issue 3, e9788. doi:10.1371/journal.pone.0009788 In this

Re: [R] (minor) R syntax error in help page to the function makeCluster of library(snow)

2012-10-04 Thread R. Michael Weylandt
On Thu, Oct 4, 2012 at 2:30 PM, Wolfgang Raffelsberger wrote: > Dear list, > > I just realized that one of the examples given in the help page to the > function makeCluster of the library(snow) has a small syntax error : > > ## to get started > library(snow) > ?makeCluster > > .. will open a halp

Re: [R] How to convert array dim names to expressions for plotting?

2012-10-04 Thread Bert Gunter
It is not clear to me what you mean by "expressions" in this context. Dimension names but be character strings -- see ?dimnames, which says: "The dimnames of a matrix or array can be NULL or a list of the same length as dim(x). If a list, its components are either NULL or a character vector with p

[R] How to convert array dim names to expressions for plotting?

2012-10-04 Thread Hofert Jan Marius
Hi, I would like to use the dim names of an array as labels in a plot. In the case of character labels, this is no problem. However, if I would like to pass expressions through array dim names, they only appear as character labels in the plot and are not evaluated as expressions (see the minima

Re: [R] smoothScatter plot

2012-10-04 Thread John Kane
Hi, Do you mean something like this? = scatter.smooth(x,y)scatter.smooth(x,y) = It looks like invoking that dcols <- densCols(x,y) is callling in some package that is masking the basic:

[R] (minor) R syntax error in help page to the function makeCluster of library(snow)

2012-10-04 Thread Wolfgang Raffelsberger
Dear list, I just realized that one of the examples given in the help page to the function makeCluster of the library(snow) has a small syntax error : ## to get started library(snow) ?makeCluster .. will open a halp page containing the command towards the end of the examples : cl <- makeCluste

Re: [R] Text to Speech In R

2012-10-04 Thread Bhupendrasinh Thakre
Integrating question and my 2 cents. Use R2WD to print your output in word file Use TTS facility of windows to do the speech Never implemented myself but may work Best Regards, Bhupendrasinh Thakre Sent from my iPhone On Oct 4, 2012, at 7:55 AM, Alex Zhang wrote: > Thanks, Mike! I am on Wi

Re: [R] Text to Speech In R

2012-10-04 Thread Jeff Newmiller
Please don't shout. Speech synthesis is associated in my mind with sound generation, which is rather platform-specific. R is a platform agnostic research tool. R can be interfaced with many kinds of external libraries, so many things are possible. However, they aren't necessarily simple, and hi

[R] Is there any package for Vector Auto-regressive with exogenous variable other than fastVAR?

2012-10-04 Thread Paul Gilbert
Package dse also does vector auto-regression with exogenous variables (and also vector ARMA and state-space models with exogenous variables). But I don't understand what you mean by "not taking the base in the model" so I don't know if it will solve your problem. Paul On 12-10-04 06:00 AM, r

Re: [R] t-test

2012-10-04 Thread arun
Hi John, You are right.  With more than two groups, the type 1 error rate should be a problem. A.K. - Original Message - From: John Kane To: arun Cc: Sent: Thursday, October 4, 2012 8:40 AM Subject: Re: [R] t-test My stats are lousy but isnt Nico doing some multiple t-tests when an

Re: [R] Text to Speech In R

2012-10-04 Thread Alex Zhang
Thanks, Mike! I am on Windows for this project tho... - Alex From: R. Michael Weylandt To: Alex Zhang Cc: "r-help@r-project.org" Sent: Thursday, October 4, 2012 8:50 AM Subject: Re: [R] Text to Speech In R On Thu, Oct 4, 2012 at 1:31 PM, Alex Zhang wrote:

Re: [R] Text to Speech In R

2012-10-04 Thread Alex Zhang
Jeff - so, do you mean that you are sure it can NOT AND SHOULD NOT? Thanks, - Alex From: Jeff Newmiller To: Alex Zhang ; Alex Zhang ; "r-help@r-project.org" Sent: Thursday, October 4, 2012 8:46 AM Subject: Re: [R] Text to Speech In R Why do you think this i

Re: [R] Problem with colors in contour plot

2012-10-04 Thread Rui Barradas
Hello, Something like this? myRamp <- colorRampPalette(c("lightblue", "darkblue")) mypal <- myRamp(20) mypal[1] <- rgb(1, 1, 1) #?filled.contour ## Persian Rug Art: x <- y <- seq(-4*pi, 4*pi, len = 27) r <- sqrt(outer(x^2, y^2, "+")) filled.contour(cos(r^2)^2, frame.plot = FALSE,

Re: [R] Text to Speech In R

2012-10-04 Thread R. Michael Weylandt
On Thu, Oct 4, 2012 at 1:31 PM, Alex Zhang wrote: > Dear all, > > Excuse me for my searching skills: I just couldn't figure out how to do any > simple text to speech in R console. > > What I want to do is, out put simple English words or sentences to speaker as > voice. No fanciness. For example

Re: [R] Class for time series

2012-10-04 Thread arun
Hi, library(xts) xts(dat1[,-1:2],order.by=as.POSIXct(paste(dat1[,1],dat1[,2],sep=" "),format="%Y-%m-%d %H:%M:%S")) #Error in .subset(x, j) : only 0's may be mixed with negative subscripts xts(dat1[,3:7],order.by=as.POSIXct(paste(dat1[,1],dat1[,2],sep=" "),format="%Y-%m-%d %H:%M:%S")) #

Re: [R] Text to Speech In R

2012-10-04 Thread Jeff Newmiller
Why do you think this is something R should be able to do? --- Jeff NewmillerThe . . Go Live... DCN:Basics: ##.#. ##.#. Live Go...

  1   2   >