[R] xls.getshlib() error in xlsReadWrite package under W7

2010-07-06 Thread milton ruser
Dear all, I installed xlsReadWrite package under Windows 7, and got error when trying to run xls.getshlib(). Please, see the errors below: > require(xlsReadWrite) Loading required package: xlsReadWrite xlsReadWrite version (cran shlib) Copyright (C) 2010 Hans-Peter Suter, Treetron, Switzerland.

Re: [R] avoiding reinstall already installed *package*

2010-05-18 Thread milton ruser
lto:r-help-boun...@r-project.org] On Behalf Of Martin Maechler > > Sent: Tuesday, May 18, 2010 1:25 PM > > To: milton ruser > > Cc: r-help@r-project.org > > Subject: Re: [R] avoiding reinstall already installed *package* > > > > On Tue, May 18, 2010 at 18:06, m

Re: [R] avoiding reinstall already installed *package*

2010-05-18 Thread milton ruser
y installed packages will not be reinstalled. cheers milton On Tue, May 18, 2010 at 12:49 PM, Martin Maechler < maech...@stat.math.ethz.ch> wrote: > { I've modified the subject; I can't stand it hitting square into > my face ... } > > >>>>> "mr&quo

Re: [R] "Re: Change class factor to numeric"

2010-05-18 Thread milton ruser
Try arcilla<-as.numeric(as.character(clay)) best milton On Tue, May 18, 2010 at 12:36 PM, Arantzazu Blanco Bernardeau < aramu...@hotmail.com> wrote: > > > sorry I had a mistake sending my question without a subject. I do resend > again. Please excuse me. > > Hello > > I have a data array with

[R] avoiding reinstall already installed library

2010-05-18 Thread milton ruser
Dear R-experts, I am installing new libraries using install.packages("ggplot2",dependencies=T). But I perceive that many dependencies are already installed. As I am using a low-band internet, how can avoid reinstall installed libraries? cheers milton [[alternative HTML version deleted]]

Re: [R] rcolorbrewer

2010-05-03 Thread milton ruser
?rev On Mon, May 3, 2010 at 4:06 PM, Steve Hempell wrote: > Can you reverse the color scheme order in rcolorbrewer. If so, how? > >[[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/l

Re: [R] sensitivity analysis, input factors

2010-04-13 Thread milton ruser
Hi Szalai I had used only "src" function, and on that case you need to have a vector with your Y variable, and a data-frame with all your X (i.e. explanatory) variables. I have interest on stay in touch with others that have been using sensitivity package! bests milton On Tue, Apr 13, 2010 at

[R] barplot with error bar in lattice

2010-04-01 Thread milton ruser
Dear all, I have a data.frame like below and I need to plot horizonal with error bar only for upper limit. On the code below I am able to plot the bars within groups, but I need (1) change from vertical to horizonal plot and (2) add the error bar. any hint are welcome. milton mydf<-read.ta

Re: [R] reading excel into R

2010-04-01 Thread milton ruser
Dear Cheba Please, install the package "xlsReadWrite" I suppose that read.csv is to read csv files, not xls ones. cheers milton On Thu, Apr 1, 2010 at 6:19 AM, cheba meier wrote: > Dear all, > > I am new R user and I am sure that this question has been asked quite often > and I have also googl

[R] qplot(ggplo2): axis label sizes/colors

2010-03-23 Thread milton ruser
Dear all, How can I change the size and the color of axis? I would like that the xlab to be larger, xvalues to be larger, and xvalues in black instead of grey. x=runif(10) y=runif(10) require(ggplo2) qplot(x,y) cheers milton [[alternative HTML version deleted]] ___

Re: [R] labels gone

2010-03-15 Thread milton ruser
Hi Zhongyi, I must confess I not understood completely what you need, but... Tau<-seq(0.05,0.95,0.05); Pi <- seq(0.19,0.01,-0.01); par(cex.axis=0.8,ps=9,mar=c(1.5,1,0.5,1), oma=c(1,1,0.2,1) ,tck=-0.01); plot(Tau,Pi, type='l', xlab="Tau",ylab="Pi",col=4, xaxt="n", yaxt="n"); axis(1,labels=F) axis(

Re: [R] barplot with a probability density curve

2010-03-14 Thread milton ruser
Hi Olga, may be you can work around with the suggestions below, while other provide the best solution for your: barplot(c(1,2,3,5,2,3,1),names.arg=c("100","200","300","400","500","600","700"), xlab="diameter",ylab="flow",main="some title",space=0.1) lines(lowess(c(1,2,3,5,2,3,1),f=1), type="b",

Re: [R] modify the levels of a variable

2010-03-10 Thread milton ruser
Try data2$X<-factor(data2$X) By the way, avoid call data as "data", matrix as "matrix", vector as "vector". cheers milton On Wed, Mar 10, 2010 at 4:35 PM, clee wrote: > > If anyone can help with this, I would greatly appreciate it. > > When I subset a dataframe, and then create a table of on

Re: [R] how to convert character variables into numeric variables directly

2010-03-08 Thread milton ruser
Try: as.numeric(as.character( typec)) milton On Mon, Mar 8, 2010 at 12:55 PM, Xumin Zeng wrote: > Here is the example. > > > age=18:29 > > height=c(76.1,77,78.1,78.2,78.8,79.7,79.9,81.1,81.2,81.8,82.8,83.5) > > type=c("A", "B", "C", "D","A", "B", "C", "D","A", "B", "C", "D") > > typec=c("0","4",

[R] parsing strings between [ ] in columns

2010-02-18 Thread milton ruser
Dear all, I have a data.frame with a column like the x shown below myDF<-data.frame(cbind(x=c("[[1, 0, 0], [0, 1]]", "[[1, 1, 0], [0, 1]]","[[1, 0, 0], [1, 1]]", "[[0, 0, 1], [0, 1]]"))) > myDF x 1 [[1, 0, 0], [0, 1]] 2 [[1, 1, 0], [0, 1]] 3 [[1, 0, 0], [1, 1]] 4 [[0, 0,

Re: [R] Separating columns, and sorting by rows

2010-02-14 Thread milton ruser
Hi Raging Jim may be this is a starting point. myDF<-read.table(stdin(),head=T,sep=",") mm,Rainfall 1977-02,17.4 1977-03,34.0 1977-04,26.2 1977-05,42.6 1977-06,58.6 1977-07,23.2 1977-08,26.8 1977-09,48.4 1977-10,47.0 1977-11,37.2 1977-12,15.0 1978-01,2.6 1978-02,6.8 1978-03,9.0 1978-04,46.6

Re: [R] release memory

2010-02-08 Thread milton ruser
Try gc() best milton On Mon, Feb 8, 2010 at 5:16 PM, dys0022 wrote: > > I am a new R user with the latest Ubuntu release. > My executions consume large amount of memory (up to 1.6 GB). When I try to > release the memory using "rm(list=ls())", R still occupies 1.6GB. (I also > tried to use rm on

Re: [R] How to choose "seed" in set.seed(seed, ...)?

2010-02-03 Thread milton ruser
Hi Michael, Experts will answer you better than me, but IMHO there is no "best practices". In general you deal with set.seed() when you want that others (or yourself) reproduce the same sequence of random numbers generation as you run on a moment. Best regards milton On Wed, Feb 3, 2010 at 1

Re: [R] R Output and ArcGIS

2010-01-25 Thread milton ruser
Hi Steve, As I have no problem with this because I run windows, I prefere to use .txt files. But may be on your case you install xlsReadWrite (or WriteRead), because .xls file works fine for ArcGis. Remember not use complex (like having dots) or large names because some times you get error on Arc

Re: [R] Logistic regression

2010-01-21 Thread milton ruser
Dear someone, Take a look at drc and Hmisc/Design packages. Next time could be polite you identify yourself bests milton On Thu, Jan 21, 2010 at 4:17 PM, DispersionMap wrote: > > can you do Logistic regression in R, if so how do you do it and how do you > test the fit of a model? > -- > View t

Re: [R] How to write '"' to a csv with the default setting of write.csv?

2010-01-21 Thread milton ruser
Peng, Are you having doble quotes? It is what you want? Did you tried the suggestion? bests milton On Thu, Jan 21, 2010 at 6:38 PM, Peng Yu wrote: > On Thu, Jan 21, 2010 at 5:32 PM, milton ruser > wrote: > > Hi Peng, > > > > May be write.table(...,quote=F) >

Re: [R] How to write '"' to a csv with the default setting of write.csv?

2010-01-21 Thread milton ruser
Hi Peng, May be write.table(...,quote=F) bests milton On Thu, Jan 21, 2010 at 6:29 PM, Peng Yu wrote: > Please see the following example. I can not write '"' to a csv file > successfully. Could somebody let me if it is possible to write '"' to > a csv file with the default setting of write.cs

Re: [R] What is the difference between S and R?

2010-01-21 Thread milton ruser
Hi There, If you are starting with S or R, don't be worried. If you just started (I suppose with R), don't be worried too. Just enjoy and be happy :-) bests milton On Thu, Jan 21, 2010 at 10:44 AM, Peng Yu wrote: > On Thu, Jan 21, 2010 at 9:07 AM, Giovanni Petris wrote: > > > > I have found

Re: [R] 3d trend line

2010-01-21 Thread milton ruser
Hi mlcarte3 (it is you?) ?contour Try algo google R graph gallery bests milton On Thu, Jan 21, 2010 at 2:43 PM, mlcarte3 wrote: > > I have a 3d scatter plot and I need to add a trend line. I can add a > regression plane, but I need a line, not a plane. > > Thanks. > -- > View this message in

Re: [R] concatenate 2 column vectors

2010-01-20 Thread milton ruser
Hi Magali, I confess that I still not figured out what you want, but try this: datatry2=transform(datatry,DIS=as.numeric(as.character(paste(datatry[,1],datatry[,2],sep="" bests milton On Wed, Jan 20, 2010 at 11:14 PM, teurlai magali wrote: > Sorry I forgot the subject in my previous post

Re: [R] how to open excel 2007 (.xlsx) file in R

2010-01-19 Thread milton ruser
Dear V.K., 1. http://www.R-project.org/posting-guide.html 2. it is - in general - good someone that write to the list identify your self. It is a polite way of participate of the list. bests milton On Tue, Jan 19, 2010 at 12:15 AM, vikrant wrote:

Re: [R] How to detect and exclude outliers in R?

2010-01-19 Thread milton ruser
Hi V.S., Did you search first on r-repositories about this issue prior to ask? May be not. RSiteSearch("outliers") bests milton On Tue, Jan 19, 2010 at 1:08 AM, vikrant wrote: > > Suppose I am reading data from a file and the data contains some outliers. > I > want to know if it is possible

Re: [R] (no subject)

2010-01-15 Thread milton ruser
Saleem, I have no idea about this, but may be you can use other data.frame name instead of "data" because data is a name of a function (see ?data). good luck milton On Fri, Jan 15, 2010 at 3:33 PM, saleem mukhtar wrote: > Hello, > > I am running R on cygwin for windows. > > File R1234 contains

Re: [R] Linear regression - R^2

2010-01-09 Thread milton ruser
Dear Csanad, PLEASE do read the posting guide http://www.R-project.org/posting-guide.html Reproducible code/problem are welcome. By the way, who you believe is right? bests milton On Sat, Jan 9, 2010 at 5:40 PM, Csanád Bertók wrote: > Hello, > > I wa

Re: [R] February 2010***New R Courses*** by XLSolutions Corp at 9 USA Cities: San Francisco, New York, Washington DC, Houston, Boston, Las Vegas, Seattle, etc

2010-01-07 Thread milton ruser
Dear Sue, Happy new your to you too. You will offer the course on R or S-Plus? If S-plus is the case, how deal with license? bests milton On Wed, Jan 6, 2010 at 3:24 PM, Sue Turner wrote: > Happy New Year ! > > XLSolutions February 2010 R courses schedule is now available online at > 9 > USA

Re: [R] Fwd: string

2010-01-06 Thread milton ruser
Hi Robert, I use something like this with grass-R: myvar<-paste("",sep="\"") myvar May be it works for you. bests milton On Wed, Jan 6, 2010 at 6:22 AM, wrote: > I would like to assign a variable y the string: style="color: #FF; font-weight: bold; font-size: 14pt"> > How can I do it - i

Re: [R] how to draw abline correctly?

2010-01-04 Thread milton ruser
You can also try: axis(1) axis(2, las=1) axis(2,at=c(0.1361), las=1) bests milton On Sun, Jan 3, 2010 at 9:26 PM, aegea wrote: > > Hello, I am frastruated with this graph, just cannot get what I need. Thank > you for any suggestions or help. I really appreciate it. I wrote the > following

Re: [R] how to draw abline correctly?

2010-01-03 Thread milton ruser
Hi there, May be part of your problem is solved. size=c(80, 40, 20, 16, 10, 8, 4, 2) estimate=c(0.1353, 0.1337, 0.1237, 0.1164, 0.1058, 0.0983, 0.0759 , 0.0607) plot(size, estimate, type="o", lwd=1.5, col="blue", ylim=c(0, 0.2), xlim=c(0,80), pch=19, axes=F) abline(h=0.1361, col="red", lwd=1.5)

Re: [R] function in aggregate applied to specific columns only

2010-01-03 Thread milton ruser
You want this? > basicSubMean <- aggregate(basicSub[c("score")], by=list(basicSub$student), FUN=mean, na.rm=TRUE) > basicSubMean Group.1 score 1 1 55.0 2 2 60.0 3 3 67.5 bests milton On Sun, Jan 3, 2010 at 10:46 PM, david hilton shanabrook < dhsha...@acad.umass.edu> wrote:

Re: [R] Artificial Neural Networks

2010-01-03 Thread milton ruser
Hi Alex, I am not a NN user, but I would suggestion you give more information about what are you planning to do. May be it helps the helper. Happy New Year. miltinho On Sun, Jan 3, 2010 at 3:54 PM, Alex Olafson wrote: > Hi! I am studying to use some R libraries which are applied for working >

Re: [R] Question on Reduce + rollmean

2010-01-02 Thread milton ruser
Dear M.Rahiz, Unfortunatelly I can't reproduce your example. PLEASE do read the posting guide http://www.R-project.org/posting-guide.html bests milton On Sat, Jan 2, 2010 at 9:26 AM, Muhammad Rahiz wrote: > Hello useRs, > > I'd like to perform a mov

Re: [R] Boxplot- input the median point and the median value

2010-01-01 Thread milton ruser
Dear Elaine, The posting guide is http://www.R-project.org/posting-guide.html bests milton On Fri, Jan 1, 2010 at 8:53 PM, elaine kuo wrote: > Dear, > > I am a newbie to R. > > Now I am learning to draw boxplot using graphics(). > > I want to high

Re: [R] How to calculate density function of Bivariate binomial distribution

2010-01-01 Thread milton ruser
Hi Nykee, I checked out ??bivariate on my R installed libraries, and found about a hundred of occurrences of "bivariate". If nobody reply you with a streigth answer, give a look at: VGAM; splancs; sm; prada; geoR; fUtilities; fBasics; akima; adehabitat; vegan; ade4 and spatstat packages. bests m

Re: [R] How to not to terminate read.table if the input file is empty?

2010-01-01 Thread milton ruser
?try() On Fri, Jan 1, 2010 at 3:41 PM, Peng Yu wrote: > read.table terminates the program if the input file is empty. Is there > way to let the program continue and return me a NULL instead of > terminating the program? > > $ Rscript read_empty.R > > read.table("empty_data.txt") > Error in re

Re: [R] How to use read.table with Hebrew column names ?

2010-01-01 Thread milton ruser
Hi Tal, You want the colnames in hebrew? Case not, may be you can use colnames() to attrib new colnames. bests milton On Fri, Jan 1, 2010 at 4:47 PM, Tal Galili wrote: > Hello dear R help group, > > I am trying to read a .txt file, with Hebrew column names, while keeping > the > column names

Re: [R] How x[, 'colname1'] is implemented?

2009-12-31 Thread milton ruser
Hi Peng, If I undertood your point, try this: x<-runif(10) y<-runif(10) z<-runif(10) w<-runif(10) myDF<-data.frame(cbind(x,y,z,w)) myDF myDF[,c("w","z")] Happy new year miltinho On Thu, Dec 31, 2009 at 5:15 PM, Peng Yu wrote: > I don't see where describes the implementation of '[]'. > >

Re: [R] joining two lists...

2009-12-31 Thread milton ruser
Dear Steve, If I undertood you have a list of "data.frame"?. Have you tried: NewData<-rbind(data.frame(myData[1]),data.frame(myData[2])) Please, send us a str(object) so the experts can help you bests milton On Thu, Dec 31, 2009 at 1:05 AM, stevemew wrote: > > My "R" problem.. > > I want

[R] Fwd: Negbin Error Warnings

2009-12-30 Thread milton ruser
: milton ruser Dear Milton, Thanks for your email. I ran the NBR-model and then tried to compute the first differences of some expected values: m3.x1 <-setx(m3, advisory=1) m3.x2 <-setx(m3, advisory=0) m3.s1 <-sim (m3,x=m3.x1, x1=m3.x2) However, I then get 50 error warnings, and there

Re: [R] Negbin Error Warnings

2009-12-30 Thread milton ruser
Dear Clara, could you, please, share the summary() of your data. bests milton 2009/12/30 Clara Brück > Hi, > > I ran a negative binomial regression (NBR) using the Zelig-package and the > negbin model. > When I then try to use the simumlation approach using the setx () and sim() > functions t

Re: [R] pass functions and arguments to function

2009-12-29 Thread milton ruser
Hi Hao, Ok. Sorry for my last post. bests milton On Tue, Dec 29, 2009 at 3:08 PM, Hao Cen wrote: > Hi, > > I wonder how to pass several functions and their arguments as arguments to > a function. For example, the main function is > > f = function(X ) { > process1(X) > ... > > > process

Re: [R] (no subject)

2009-12-29 Thread milton ruser
Hi Hao, I suggest you try again, starting by read posting guide at footnote of this email. How about a title for the message? How about identify yourself? bests milton On Tue, Dec 29, 2009 at 2:59 PM, Hao Cen wrote: > Hi, > > I wonder how to pass several functions and their arguments as argum

Re: [R] 2D array of strings

2009-12-28 Thread milton ruser
Hi Francesco, Be carefull with create a object named str, because you crash str() function. I don't know if it have implications on any package or functions. bests milton On Mon, Dec 28, 2009 at 4:32 PM, Francesco Napolitano wrote: > Sorry for the dumb question, but I couldn't figure this out

Re: [R] help: creating a unified histogram

2009-12-27 Thread milton ruser
Hi Yael, Have you tried to find it on google "R graph gallery" ? bests milton On Sun, Dec 27, 2009 at 2:45 PM, Yael Yohai wrote: > Good evening, > > I would like to put the histograms of several matrices on the same graph, > together on a unified histogram (instead of having one bar per value

Re: [R] help in merging

2009-12-24 Thread milton ruser
Hi there, You can add a order column on x or y and after use this field to order z. Like z<-z[order(z$orderfield),] To generate a order on x or y you can do something like x$xorder<-1:nrow(x) cheers milton On Thu, Dec 24, 2009 at 2:26 PM, wrote: > > Hi All, > I want to "merge" two datase

Re: [R] aggregate binary response data

2009-12-24 Thread milton ruser
may be ?table ? On Thu, Dec 24, 2009 at 7:56 AM, Graham Leask wrote: > Dear list > > I have a response variable coded 0/1 i.e. a binary response. There are > 20,000 individual responses that I would like to aggregate into numbers of > each category (i.e. 0/1) by > group called dn (350 different g

Re: [R] aggregate binary response data

2009-12-24 Thread milton ruser
?aggregate ? On Thu, Dec 24, 2009 at 7:56 AM, Graham Leask wrote: > Dear list > > I have a response variable coded 0/1 i.e. a binary response. There are > 20,000 individual responses that I would like to aggregate into numbers of > each category (i.e. 0/1) by > group called dn (350 different grou

[R] prcomp : plotting only explanatory axis arrows

2009-12-23 Thread milton ruser
Dear all, I have a very large dataset (1712351 , 20) and would like to plot only the arrows that represent the contribution of each variables. On the sample below I woild like to plot only the explanatory variables (Murder, Assault..) and not the sites. prcomp(USArrests) # inappropriate prcomp(U

[R] removing rows with NAs in anywere

2009-12-22 Thread milton ruser
Dear all, I have a data.frame with some NAs that can happens anywere, and I would like to keep only rows without NAs. Thanks in advance for your help, and Merry Xmas. myvect<-runif(100) myvect[sample(1:100)[1:5]]<-NA myDF<-data.frame(matrix(myvect,ncol=5)) myDF miltinho [[alternative HTM

Re: [R] Help,Suggest me some methods to identify tr aining set and test set!!!

2009-12-21 Thread milton ruser
Not ellegant.. but... MyDF<-data.frame(cbind(x=runif(10), y=rnorm(10))) TrainingSize=5 TrainingSize_list<-sample(1:nrow(MyDF))[1:TrainingSize] TrainingSize_list MyDF.training<-MyDF[(1:nrow(MyDF) %in% TrainingSize_list),] MyDF.training MyDF.test<-MyDF[ ! (1:nrow(MyDF) %in% TrainingSize_list),]

Re: [R] write.csv and col.names=F

2009-12-17 Thread milton ruser
?write.table On Thu, Dec 17, 2009 at 6:52 PM, kayj wrote: > > Hi All, > > I always have a problem with write.csv when I want the column names to be > ignored, when I specify col.names=F, I get a header of V1 V2 V3 V4 etc. > > for example I tried > > write.csv(mydata, file="data.csv", quote=FALSE

Re: [R] how to import data from excel to R

2009-12-17 Thread milton ruser
Hi Sarath, try "C:\\Documents and Settings\\admin\\Desktop\\data.txt" or "C:/Documents and Settings/admin/Desktop/data.txt" or yet setwd("your path") mydata<-read.table("data.txt") cheers miltinho On Thu, Dec 17, 2009 at 12:24 PM, wrote: > Hi, > I am using R and I want to know how data can

Re: [R] help with graphing -- Points in my graph are not apparent, always displayed in steps

2009-12-12 Thread milton ruser
Hi Philip, I must confess that I not understood what is the problem. Could you clarify it a little bit more? Cheers miltinho brazil=toronto On Sat, Dec 12, 2009 at 12:39 AM, philip robinson wrote: > > I am trying to graphically represent a large set of data who's result is > not > strictly un

Re: [R] how to creat a matrix

2009-12-11 Thread milton ruser
Hi Him. Did you read the Help/Manual in PDF/A Instroduction to R.PDF? MyVect<-scan() 0 0 0 0 1 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 MyVect MyMat<-matrix(MyVect, ncol=5, byrow=T) MyMat bests milton On Fri, Dec 11, 2009 at 10:47 AM, Moohwan Kim wrote: > Dear R family > > I am attempting

Re: [R] Population Histogram

2009-12-09 Thread milton ruser
1. Read the posting guide: http://www.R-project.org/posting-guide.html 2. Did you installed R? 2.1. Case yes, go to the help, click on "Manual in PDF" bests miltinho On Wed, Dec 9, 2009 at 1:47 PM, terry johnson wrote: > How would I make a popula

Re: [R] hi there!

2009-12-06 Thread milton ruser
Hi Dave, May be a good idea you (1) use subjects avoiding "help", "Hi there", "problem", etc. The list is archived for future search, and nobody with problems with MCMC, Metropholis Hasting algorithm will search "Hi there"; (2) how about minimum reproducible problem? If you are needing use this ki

Re: [R] 4D plotting

2009-12-06 Thread milton ruser
Hi Jared, google R graph library 4d or http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=90 bests milton On Sun, Dec 6, 2009 at 7:31 AM, Jared Nance wrote: > Hello list, > Thanks in advance for whatever help you can give me. I have a data set > that I want to visualize that has fo

Re: [R] length of a density curve (or any curve)

2009-12-04 Thread milton ruser
hi Sylvain, did you try ?density regards milton On Fri, Dec 4, 2009 at 7:19 AM, sylvain willart wrote: > Hello R users, > > When I type > > d <- density(MyData$x) > > I obtain a density object I can plot, > > But I wonder if there is a way to easily compute the length of the > density curve ?

Re: [R] how to seperate a matrix

2009-12-03 Thread milton ruser
Hi there, mymat<-matrix(c(6,8,1, 5,9,2, 20,10,3, 7,11,4, 8,12,5, 25,13,6, 14,14,7), byrow=T, ncol=3) mymat1<-mymat[mymat[,1]=mymat[,2],] mymat2 bests milton On Thu, Dec 3, 2009 at 10:12 PM, aegea wrote: > > Hello, > > I am working on seperate the matrix to two matrices but got trouble on > do

Re: [R] sort a data frame by a vector

2009-12-01 Thread milton ruser
not ellegant.. but... vecDF = data.frame(A1=c("C", "A", "B")) vecDF$A1.order=1:dim(vecDF) vecDF dataDF = data.frame(A1 = c("B", "A", "C"), A2 = c(1,2,3)) dataDF2<-merge(vecDF, dataDF, by=intersect(colnames(vecDF),colnames(dataDF))) dataDF2 dataDF2.ord<-dataDF2[order(dataDF2$A1.order),] dataDF2.ord

Re: [R] Normal tests disagree?

2009-12-01 Thread milton ruser
Hi there, Could you provide a minimum reproducible code, please. Bests milton On Tue, Dec 1, 2009 at 6:11 PM, wrote: > If I have data that I feed into shapio.test and jarque.bera.test yet they > seem to disagree. What do I use for a decision? > > For my data set I have p.value of 0.05496421 re

Re: [R] Ggplot2: Elegant Graphics for Data Analysis (Use R) (2009 Paperback)

2009-11-30 Thread milton ruser
I agree with Ista. The book is very interesting. I put my 2cents too. milton On Mon, Nov 30, 2009 at 6:20 PM, Ista Zahn wrote: > Depends on what you want. It's possible to find most of the > information contained in the book in other places (the ggplto2 > website, Hadley's dissertation, the gg

Re: [R] Arrhenius plot with lattice

2009-11-25 Thread milton ruser
Hi Markus, PLEASE do read the posting guide http://www.R-project.org/posting-guide.html Your example is not reproducible. Help your helper :-) bests miltinho On Wed, Nov 25, 2009 at 5:54 PM, Markus Häge wrote: > hello there, > > I like to print a t

Re: [R] 11 distinguishable colors

2009-11-24 Thread milton ruser
Hi there, I picked it up from r.colors' grass manual: http://colorbrewer2.org/ The code bellow cold be usefull for you explore simbols and colors. x11(900,500) par(mfrow=c(1,2)) x<-runif(20) y<-runif(20) plot(y~x, type="n") Number.of.symbols<-20 for (i in 1:Number.of.symbols) { points(y[i]~x[

[R] (ecological) model sensitivity analysis book

2009-11-20 Thread milton ruser
Dear all, I know that my question may be so general, but I would like to know if is there a good book for model sensitivity analysis (applied to ecology?) that use R? Is there any package dedicated to model sensitivity analysis on R? Cheers milton [[alternative HTML version deleted]] _

Re: [R] Presentation of data in Graphical format

2009-11-15 Thread milton ruser
Google "R graph grallery" Google "R ggplot2" Google "R lattice" and good luck milton On Sun, Nov 15, 2009 at 7:48 AM, Sunita22 wrote: > > Hello > > My data contains following columns: > > 1st column: Posts (GM, Secretary, AM, Office Boy) > 2nd Column: Dept (Finance, HR, ...) > 3rd column: Tasks

Re: [R] refactoring in R

2009-11-13 Thread milton ruser
Hi Peng, If that information is preliminary, so I guess you have a more clear problem and may be you are able to state a minimally reproducible code/example with what you really need. Bests milton On Fri, Nov 13, 2009 at 9:30 PM, Peng Yu wrote: > I'm wondering if there are some tips for ref

Re: [R] LINEAR MIXED EFFECT

2009-11-13 Thread milton ruser
all your data "data" you will crach data() function. cheers milton On Fri, Nov 13, 2009 at 8:30 AM, ANARPCG wrote: > > > > Douglas Bates-2 wrote: > > > > On Thu, Nov 12, 2009 at 10:14 AM, milton ruser > > wrote: > >> Hi Ana, > >> >

Re: [R] LINEAR MIXED EFFECT

2009-11-12 Thread milton ruser
know what function you are using (see first point). > Regardless, sounds like you have missing data and na.action is set to > na.fail (ie, "fail" if any missing data). > > cheers, Dave > > > milton ruser wrote: > > > > Dear Ana Golveia, > > > &

Re: [R] LINEAR MIXED EFFECT

2009-11-11 Thread milton ruser
Dear Ana Golveia, It is completelly impossible someone realise what kind or help you need or what is happening. I suggest you give a look on the posting guide, mainly that part about a minimum reproducible code with self explaining information, etc. Cheers milton On Wed, Nov 11, 2009 at 7:22 AM

[R] trend thresholds on time series

2009-11-07 Thread milton ruser
Dear all, I have several response variables estimated from some simulations, and I would like to identify the thresholds for trend changes. Fro example, below I forced two different response behaviours and on x is time unit. x<-1:1500 y<-x/exp(x^0.2) smaller15<-y[y<15] y<-ifelse(y<15,y+rnorm(len

Re: [R] counting frequencies across two columns

2009-10-31 Thread milton ruser
Hi Jason, As your example is not reproducible, may be something like: myFreq<-data.frame(table(articleID, author)) if you want to know only those articles with 1 author, you can try subset(myFreq, Freq==1) or something like. bests milton On Sun, Nov 1, 2009 at 2:20 AM, Jason Priem wrote:

Re: [R] x-y plot as an histogram

2009-10-29 Thread milton ruser
Hi Federico, Have you checked google R graph library? bests milton On Thu, Oct 29, 2009 at 3:12 PM, Federico Abascal wrote: > Hi, > > I am investigating a problem for which I found no solution. > I have a matrix with two columns. I have plotted it as an x-y plot > (plot(matrix[,1],matrix[,2])

Re: [R] rotate levelplot

2009-10-24 Thread milton ruser
Hi Kang, Could you send a reproducible sample-code? Bests miltinho On Sat, Oct 24, 2009 at 11:32 PM, Kang Min wrote: > Hi all, > > I have a matrix with 5 rows and 10 columns, which represent the grids > on a rectangular map. > I used the code below to plot, but it gives me the map with the 1

Re: [R] updating columns using other column as reference

2009-10-19 Thread milton ruser
Amazing. Thanks a lot milton On Mon, Oct 19, 2009 at 12:55 PM, Henrique Dallazuanna wrote: > Milton, > > Try this: > > myDF[grep("k", names(myDF))] <- sapply(myDF[grep("k", names(myDF))], > function(idx)myDF$namesp[idx]) > > > On Mon, Oct 19,

[R] updating columns using other column as reference

2009-10-19 Thread milton ruser
Dear R-gurus, Just supose I have a dara.frame that looks like myDF<-read.table(stdin(),head=T,sep=",") codID,namesp,k1,k2,k3,k4 1,spA,2,5,6,3 2,spB,4,5,4,6 3,spC,2,1,5,6 4,spC,5,4,3,2 5,spD,1,2,3,4 6,spE,2,4,3,1 I need to update the columns k1-k4 with the namesp, but considering the math between

Re: [R] Filtering on a dataframe- newbie question

2009-10-18 Thread milton ruser
Hi there, Try ?subset ?str may also be useful.. bests milton On Sun, Oct 18, 2009 at 11:10 PM, ANJAN PURKAYASTHA < anjan.purkayas...@gmail.com> wrote: > Hi, > newbie question. I have a data-frame with 3 named columns: Name, Obs1, > Obs2. > The Name column members are made of alphanumeric chara

Re: [R] is that possible to graph 4 dimention plot

2009-10-04 Thread milton ruser
How about : google "r graph gallery 4d" ? abs On Sun, Oct 4, 2009 at 3:31 PM, Duncan Murdoch wrote: > On 04/10/2009 3:14 PM, gcheer3 wrote: > >> Suppose there are 4 variables >> d is a function of a , b and c >> I want to know how a, b and c change will make d change >> It will be straightforwa

Re: [R] Can You Recommend Books for Linear Mixed Models in R

2009-10-01 Thread milton ruser
Hi Pat, I am happy with the book "Mixed Effects Models and Extensions in Ecology with R" by Zuur, Ieno, Walker, Saveliev and Smith. In fact I am reading chapter by chapter, and the reading is very digestible, as well as the examples are quite easy to understand and to be reference for our "real wo

Re: [R] Deleting a column in a dataframe by name

2009-09-29 Thread milton ruser
Peace on your heart... it is still Tuesday, and Friday is away... :-) On Tue, Sep 29, 2009 at 4:54 PM, Rolf Turner wrote: > > On 30/09/2009, at 9:32 AM, milton ruser wrote: > > x=runif(12) >> y=runif(12) >> w=runif(12) >> >> mydf<-data.frame(cbind(x,y,w)

Re: [R] Deleting a column in a dataframe by name

2009-09-29 Thread milton ruser
x=runif(12) y=runif(12) w=runif(12) mydf<-data.frame(cbind(x,y,w)) head(mydf) mydf<-subset(mydf, select=c(-x,-w)) head(mydf) On Tue, Sep 29, 2009 at 4:27 PM, Rolf Turner wrote: > > On 30/09/2009, at 9:15 AM, milton ruser wrote: > > May be this: >> >> FRAME &l

Re: [R] Deleting a column in a dataframe by name

2009-09-29 Thread milton ruser
May be this: FRAME <- FRAME[-c(NAME1, NAME2)] or FRAME<- subset(FRAME, select=(-NAME1, -NAME2)) bests milton On Tue, Sep 29, 2009 at 4:12 PM, Dennis Fisher wrote: > Colleagues, > > Hopefully a simple problem: I want to delete a column with a known name > from a dataframe. I could write: > >

Re: [R] Count number of zeros in a collumn

2009-09-27 Thread milton ruser
Hi Marcio, How about if (min(dart[,1977])==0)... cheers milton On Sun, Sep 27, 2009 at 4:54 PM, Marcio Resende wrote: > > I have a matrix 700x2000 which is sampled in each cycle from another matrix > 788x2000 with the numbers 0,1 and 9 > > There is one specific collumn of this matrix, dart[,197

Re: [R] Teach me how to transpose in R

2009-09-26 Thread milton ruser
*Hum* bbb=t(as.matrix(data2)) ? good luck milton On Sun, Sep 27, 2009 at 12:39 AM, Hyo Lee wrote: > Hi guys, > I need your help!! > > My goal is to make a csv file from ncdf file. > This is the code i've used : > > > hyo=open.ncdf("C:/CRUTEM3.nc") > > hyo > [1] "file C:/CRUTEM3.nc has 4 di

Re: [R] Correlate two time series

2009-09-22 Thread milton ruser
Have you looked at ts() ? On Tue, Sep 22, 2009 at 5:49 PM, Maggie wrote: > I was wondering how to correlate two time series in R? I have to > plotted waver files I need to correlate to one another to see how well > they align.. > > Any guidance would be very much appreciated! > > ___

Re: [R] Superimposing (aligning) two graphs

2009-09-22 Thread milton ruser
Hi Maggie: 1) plot() combined with lines() may do the job. 2) google "R graph gallery" bests milton On Tue, Sep 22, 2009 at 5:23 PM, Maggie wrote: > Hello, I need to superimpose two graphs (plotted representations of > waver files) so I can see how much they align. Is there a quick and > di

Re: [R] Create directory and copy files in R

2009-09-22 Thread milton ruser
Tammy Try avoid directories with spaces... may be this solve your problem. Case not, write us again. bests milton On Mon, Sep 21, 2009 at 10:05 AM, Tammy Ma wrote: > > > HI, All R users, > > My problem is: > > > fn > [1] "C:/Documents and > Settings/lma/Desktop/FamilyAEntrepreneurs/Entrepreneu

Re: [R] plotting least-squares regression against x-axis

2009-09-19 Thread milton ruser
By the way, how about: plot(y,mylm$residuals) So you can have an idea of your error distribution across your response variable domain. :-) milton toronto=brazil On Sat, Sep 19, 2009 at 9:02 PM, Sunil Suchindran wrote: > x <- seq(50) > > y <- 10 + x * 2 + rnorm(50,0,10) > > plot(y~x) > > mylm =

Re: [R] (no subject)

2009-09-13 Thread milton ruser
Hi Karen, As you are just starting (I guess), I suggest you take a shortcut: 1. copy and paste the data from Excel to the notepad, and save it as .TXT 2. If all is ok, you will have a tab-delimited file. 3. Be sure that your variables have its names on first line 4. read your data using setwd("c:

Re: [R] Joining Characters in R {issue with paste}

2009-09-09 Thread milton ruser
You not tryed ?paste :-) paste(a,b,sep="") bests milton On Wed, Sep 9, 2009 at 5:08 PM, Abhishek Pratap wrote: > Hi Guys > I am want to join to strings in R. I am using paste but not getting > desirable result. > > For the sake of clarity, a quick example: > > > a="Bio" > > b="iology" > > pas

Re: [R] fitting nonlinear model

2009-09-09 Thread milton ruser
, data=mydf.new, lty=2) cheers milton On Wed, Sep 9, 2009 at 3:01 PM, milton ruser wrote: > Hi Bill, > > I am not sure what you want, but... > > mydf<-read.table(stdin(), head=T, sep=",") > Np,year > 96,2 > 91,5 > 89,7 > 85,10 > > plot(Np~year,

Re: [R] fitting nonlinear model

2009-09-09 Thread milton ruser
Hi Bill, I am not sure what you want, but... mydf<-read.table(stdin(), head=T, sep=",") Np,year 96,2 91,5 89,7 85,10 plot(Np~year, data=mydf) mymodel<-lm(Np~year, data=mydf) abline(mymodel, col="red") bests milton On Wed, Sep 9, 2009 at 2:27 PM, Bill Hyman wrote: > My data look like: > > Np

Re: [R] (no subject)

2009-09-06 Thread milton ruser
May be you also want: x<-runif(1000,min=0, max=600) x3<-x[x<=300] summary(x3) hist(x3, ,breaks = c(0,50,100,150,200,250,300)) bests milton On Mon, Sep 7, 2009 at 1:01 AM, milton ruser wrote: > Hi Roslina, > > x<-runif(1000,min=0, max=600) > hist(x) > > x2&l

Re: [R] (no subject)

2009-09-06 Thread milton ruser
Hi Roslina, x<-runif(1000,min=0, max=600) hist(x) x2<-ifelse(x<=300,x,300) summary(x2) hist(x2, ,breaks = c(0,50,100,150,200,250,300)) good luck milton On Mon, Sep 7, 2009 at 12:51 AM, Roslina Zakaria wrote: > Hi r-users, > > I would like to know how to put all the data that is greater than

Re: [R] Easy way to get top 2 items from vector

2009-09-03 Thread milton ruser
Hi Noah, Next time try, please, send a short reproducible code/example. May be this can be (not so elegant but) helpfull myDF<-data.frame(cbind(a=runif(5),b=runif(5))) myDF N=2 a.order<-rev(order(myDF$a))[1:N] b.order<-rev(order(myDF$b))[1:N] myDF.max2a<-myDF[a.order,] myDF.max2a myDF.max2b<-m

Re: [R] variable selection in logistic

2009-09-03 Thread milton ruser
Hi Annie, What kind of data (response and explanatory) you have? As a ecological modeller, I always think first what is my (our) mind models. After that I analyze a set of concorrent - *but with ecological meaning* - models. It helps me to test hypothesis as well as to discuss the results. If you

  1   2   3   >