[R] RWeka prediction

2009-04-26 Thread Dajiang Liu
Dear All,I encountered a problem when I use RWeka for prediction. Specifically, I use the following: res=J48(X1~.,data=mydata); predict(res), #it worked fine but when I tried to use a different data set, i.e. predict(res,newdata=mynewdata); all the predictions I get is 0, which apparently is

[R] RWeka prediction

2009-04-26 Thread Dajiang Liu
Dear All,I encountered a problem when I use RWeka for prediction. Specifically, I use the following: res=J48(X1~.,data=mydata); predict(res), #it worked fine but when I tried to use a different data set, i.e. predict(res,newdata=mynewdata); all the predictions I get is 0, which apparently is

[R] help with plotting results of lda

2009-04-26 Thread pgseye
Hi, I've performed an lda and obtained a classification table for some of my data: efa.dfa-lda(groups~.,efa.scores.8,CV=T) str(efa.dfa) List of 5 $ class: Factor w/ 2 levels 1,2: 1 2 1 2 1 1 2 2 1 2 ... $ posterior: num [1:160, 1:2] 0.99083 0.00852 0.93983 0.23186 0.85931 ... ..-

Re: [R] argument 'exclude' in function xtabs

2009-04-26 Thread Matthieu Lesnoff
Ok thanks Erik Indeed, exclude in xtabs seems working whith character vectors: x - c(rep(c(A,B,C), 2)) x [1] A B C A B C xtabs(~ x, exclude = B) x A C 2 2 and not directly with factors: x - factor(rep(c(A,B,C), 2)) x [1] A B C A B C Levels: A B C levels(x) [1] A B C xtabs(~ x,

Re: [R] omit empty cells in crosstab?

2009-04-26 Thread Bernardo Rangel Tura
On Fri, 2009-04-24 at 13:12 -0700, sjaffe wrote: small example: a-c(1.1, 2.1, 9.1) b-cut(a,0:10) c-data.frame(b,b) d-table(c) dim(d) ##result: c(10, 10) But only 9 of the 100 cells are non-zero. If there were 10 columns, the table have 10 dimensions each of length 10, so have 10^10

Re: [R] RWeka prediction

2009-04-26 Thread Achim Zeileis
On Sun, 26 Apr 2009, Dajiang Liu wrote: Dear All,I encountered a problem when I use RWeka for prediction. Specifically, I use the following: res=J48(X1~.,data=mydata); predict(res), #it worked fine but when I tried to use a different data set, i.e. predict(res,newdata=mynewdata); all the

Re: [R] Changing gird marks in ggplot2

2009-04-26 Thread ONKELINX, Thierry
Dear Chris, Changing coord_cartesian(ylim = c(0, 5)) into coord_cartesian() + scale_y_continuous(limits = c(0, 5)) That should solve your problem. HTH, Thierry ir. Thierry Onkelinx

Re: [R] argument 'exclude' in function xtabs

2009-04-26 Thread Matthieu Lesnoff
Ok thanks Erik Indeed, exclude in xtabs seems working whith character vectors: x - c(rep(c(A,B,C), 2)) x [1] A B C A B C xtabs(~ x, exclude = B) x A C 2 2 and not directly with factors: x - factor(rep(c(A,B,C), 2)) x [1] A B C A B C Levels: A B C levels(x) [1] A B C xtabs(~ x, exclude =

[R] Conditional plot labels

2009-04-26 Thread Christian Bustamante
Hi all, I'm trying to do multiple graphs in a window like this: ___ ___ ___ ylab |__| |__| |__| ___ ___ ___ ylab |__| |__| |__| ___ ___ ___ ylab |__| |__| |__| xl xl xl If I try to put the labels manually, some graphs become smaller

[R] Question of Quantile Regression for Longitudinal Data

2009-04-26 Thread Helen Chen
Hi, I am trying to estimate a quantile regression using panel data. I am trying to use the model that is described in Dr. Koenker's article. So I use the code the that is posted in the following link: http://www.econ.uiuc.edu/~roger/research/panel/rq.fit.panel.R How to estimate the panel

[R] Scatterplot of two groups side-by-side?

2009-04-26 Thread nonu...@yahoo.de
Dear all I'm realy new to R, so I hope you can help me, as I didn't find any solution in the common books. Since some days I'm trying to create the following plot: A scatterplott showning two different groups side-by-side with according regression lines. Both datasets only have the same five

[R] THE EQUIVALENT OF SQL INNER TABLE JOIN IN R

2009-04-26 Thread Nigel Birney
Hello all, Apologize for the newbie question. What's the easiest way to do a SQL inner table join in R? Say I have a table containing column names A, B, C and another which has columns named C, D, E. I would like to do an inner table join on C and produce a table A, B, C, D, E. thanks a lot,

[R] simulate arima model

2009-04-26 Thread Rebecca1117
I am new in R. I can simulate Arma, using Arima.sim However, I want to simulate an Arima Model. Say (1-B)Zt=5+(1-B)at. I do not know how to deal with 5 in this model. Can any one could help me? Thank you very much! Regards, -- View this message in context:

Re: [R] THE EQUIVALENT OF SQL INNER TABLE JOIN IN R

2009-04-26 Thread Gabor Grothendieck
1. ?merge 2. sqldf package whose home page is at: http://sqldf.googlecode.com On Sat, Apr 25, 2009 at 9:15 PM, Nigel Birney na...@cam.ac.uk wrote: Hello all, Apologize for the newbie question. What's the easiest way to do a SQL inner table join in R? Say I have a table containing column

[R] Stochastic Gradient Ascent for logistic regression

2009-04-26 Thread Tim LIU
Hi. guys, I am trying to write my own Stochastic Gradient Ascent for logistic regression in R. But it seems that I am having convergence problem. Am I doing anything wrong, or just the data is off? Here is my code in R - lbw -

Re: [R] THE EQUIVALENT OF SQL INNER TABLE JOIN IN R

2009-04-26 Thread Peter Dalgaard
Nigel Birney wrote: Hello all, Apologize for the newbie question. What's the easiest way to do a SQL inner table join in R? Say I have a table containing column names A, B, C and another which has columns named C, D, E. I would like to do an inner table join on C and produce a table A, B, C,

Re: [R] Nomogram with stratified cph in Design package

2009-04-26 Thread reneepark
I'm sorry - I meant a median survival estimate, not a median risk. I see - I didn't realize that by stratifying it would pool the levels of the stratified variable. Hm, that is unfortunate, considering the stratified variable is one that I would like to keep in the nomogram. Thank you for your

Re: [R] Question of Quantile Regression for Longitudinal Data

2009-04-26 Thread Tirthankar Chakravarty
This is a nontrivial problem. This comes up often on the Statalist (-qreg- is for cross-section quantile regression): You want to fit a plane through the origin using the L-1 norm. This is not as easy as with L-2 norm (LS), as it is more than a matter of dropping a constant predictor yet

Re: [R] Conditional plot labels

2009-04-26 Thread baptiste auguie
Hi, Have you considered using high-level plotting functions provided by the ggplot2 or lattice package? Here's a dummy example, x - seq(0, 10, length=100) y1 - sin(x) y2 - cos(x) y3 - x^2/100 y4 - 1/x d - data.frame(x, y1, y2, y3, y4) library(reshape) dm - melt(d, id=x) dm$type1 -

Re: [R] Scatterplot of two groups side-by-side?

2009-04-26 Thread baptiste auguie
Hi, You could do this very easily using ggplot2, #install.packages(ggplot2, dep=TRUE) library(ggplot2) c - ggplot(mtcars, aes(y=wt, x=mpg)) + facet_grid(. ~ cyl) c + stat_smooth(method=lm) + geom_point() See more examples on Hadley's website: http://had.co.nz/ggplot2/ Hope this helps,

Re: [R] 3 questions regarding matrix copy/shuffle/compares

2009-04-26 Thread Esmail
Hello David, Let me try again, I don't think this was the best post ever I've made :-) Hopefully this is clearer, or otherwise I may break this up into three separate simple queries as this may be too long. == is not an assignment operator in R, so the answer is that it would do neither. -

[R] Help to select the raw in a data.frame with the max value

2009-04-26 Thread Alessandro
Dear User, thank for the attention. I have a data.frame with 5 columns (ex:ID, a1,a2,a3,a4) and 1000 rows. I wish to find the absolute max value for all data.frame and save a new data.frame with the row where is that value. Ex: ID: 1,2,3,4,5,6,7,8,9,10 a1:1,2,3,4,5,6,7,8,9,10

[R] re moving entries from one vector that are in another

2009-04-26 Thread onyourmark
I have various objects defined but I am trying to remove a set of elements in one vector from another and using the loops at the end of this post but I am getting the error at the very end of this post. str(x) num [1:923, 1:923] 1 -0.00371 -0.00102 -0.00204 -0.00102 ... - attr(*,

[R] Is their any function can generate orthogonal tables(e.g. L_8(2^7)

2009-04-26 Thread dzuswxbylw
I want to generate some orthogonal tables in my experiment design. I searched http://search.r-project.org/ (RSiteSearch) whith key words such as orthogonal table, orthogonal design, latin square etc and get no useful result. the same result get by searching via google's insite search in r-cran

Re: [R] dotplot: labeling coordinates for each point

2009-04-26 Thread Qifei Zhu
Hi David, Thanks! It looks much better now. but is there any way to add (x,y) coordinates as labels to all the points in the graph? Best case if I can enforce some conditions saying if (y10,000) label, else no label. Any advice is appreciated. Best, Tony -Original Message- From: David

Re: [R] Help to select the raw in a data.frame with the max value

2009-04-26 Thread Jorge Ivan Velez
Dear Alessandro, Here is one way: DF - data.frame(ID,a1,a2,a3,a4) Row - which( DF == max(DF[,-1]), arr.ind = TRUE)[1] DF[Row,] # ID a1 a2 a3 a4 # 10 10 10 20 30 40 See ?which and ?max for more details. HTH, Jorge On Sun, Apr 26, 2009 at 8:02 AM, Alessandro

Re: [R] Nomogram with stratified cph in Design package

2009-04-26 Thread Frank E Harrell Jr
David Winsemius wrote: On Apr 25, 2009, at 6:57 PM, reneepark wrote: Hello, I am using Dr. Harrell's design package to make a nomogram. I was able to make a beautiful one without stratifying, however, I will need to stratify to meet PH assumptions. This is where I go wrong, but I'm not

Re: [R] 3 questions regarding matrix copy/shuffle/compares

2009-04-26 Thread David Winsemius
On Apr 26, 2009, at 7:48 AM, Esmail wrote: Hello David, Let me try again, I don't think this was the best post ever I've made :-) Hopefully this is clearer, or otherwise I may break this up into three separate simple queries as this may be too long. == is not an assignment operator in

[R] Problem installing packages

2009-04-26 Thread Jarek Jasiewicz
since 2.9.0 version I have a problem with installing packages: install.packages(sp) --- Please select a CRAN mirror for use in this session --- Loading Tcl/Tk interface ... done Warning: unable to access index for repository http://piotrkosoft.net/pub/mirrors/CRAN/src/contrib Warning messages:

Re: [R] issue building my own package... moving from Apple OS to Windows

2009-04-26 Thread Uwe Ligges
Daryl Morris wrote: Thanks for the various responses. It was easier than I thought to get all the tools together and setup Windows paths for the build process. I've been successful now on Windows. BUT... how do I make a package which DOES NOT require R CMD INSTALL to install? Obviously,

Re: [R] Scatterplot of two groups side-by-side?

2009-04-26 Thread Stefan Grosse
On Sun, 26 Apr 2009 09:29:39 + (GMT) nonu...@yahoo.de nonu...@yahoo.de wrote: ND I'm realy new to R, so I hope you can help me, as I didn't find any ND solution in the common books. Have a look at the online resources at: http://cran.r-project.org/other-docs.html There is also stuff on

Re: [R] 3 questions regarding matrix copy/shuffle/compares

2009-04-26 Thread Esmail
David Winsemius wrote: Yes. As I said before I am going to refrain from posting speculation until you provide valid R code that will create an object that can be the subject of operations. The code I have provided works, here is a run that may prove helpful: POP_SIZE = 6 LEN = 8

Re: [R] Help to select the raw in a data.frame with the max value

2009-04-26 Thread David Winsemius
On Apr 26, 2009, at 8:02 AM, Alessandro wrote: Dear User, thank for the attention. I have a data.frame with 5 columns (ex:ID, a1,a2,a3,a4) and 1000 rows. I wish to find the absolute max value for all data.frame and save a new data.frame with the row where is that value. Ex: ID:

[R] R: constrained optimization

2009-04-26 Thread mauede
Thank you. Unluckily what makes the problem only apparenttly simple (for me) is that we have not differentiable functions and the parameter space is not continuous ... which reduces dramatically the number of choices. I would be grateful to chat with anyone who has tackled a similar problem.

Re: [R] Scatterplot of two groups side-by-side?

2009-04-26 Thread John Fox
Dear Karin, If I understand correctly what you want, the scatterplot function in the car package isn't designed to produce it, but there are many ways to draw side-by-side scatterplots. Here is one, using basic R graphics: par(mfrow=c(1,2)) by(Data, Data$group, function(x) {

Re: [R] Problem installing packages

2009-04-26 Thread Uwe Ligges
Jarek Jasiewicz wrote: since 2.9.0 version I have a problem with installing packages: install.packages(sp) --- Please select a CRAN mirror for use in this session --- Loading Tcl/Tk interface ... done Warning: unable to access index for repository

Re: [R] fclustindex, e1071 package

2009-04-26 Thread Uwe Ligges
Caroline Wallis wrote: Hi, I'm using e1071 package to do fuzzy cluster analysis. My dataset (ra) has 5237 observations and 2 variables - depth and velocity. I used fuzzy cmeans to create 6 fuzzy classes. ra.flcust6-cmeans(ra,6,iter.max=100,verbose=F,dist=euclidean,method=cmea

Re: [R] Overlapping parameters k in different functions in ipred

2009-04-26 Thread Uwe Ligges
WilDsc0p wrote: Dear List, I have a question regarding ipred package. Under 10-fold cv, for different knn ( = 1,3,...25), I am getting same misclassification errors: # library(ipred) data(iris) cv.k = 10 ## 10-fold cross-validation bwpredict.knn

Re: [R] help with plotting results of lda

2009-04-26 Thread Uwe Ligges
Works for me with library(MASS) plot(lda(Species~., data=iris)) hence you may want to profide the data to enable us to reproduce your problem... Uwe Ligges pgseye wrote: Hi, I've performed an lda and obtained a classification table for some of my data:

Re: [R] How to get rid of loop?

2009-04-26 Thread Uwe Ligges
Ken-JP wrote: set.seed(1) x - runif(100) # I want to calculate y such that: # # 1. if x0.75, y - 1 # 2. else if x0.25, y - -1 # 3. else if y_prev==1 x0.5, y - 0 # 4. else if y_prev==-1 x0.5, y - 0 # 5. else y - y_prev # # 1. and 2. are directly doable without looping. # # How do I do 3.-5.

Re: [R] Problem installing packages

2009-04-26 Thread Jarek Jasiewicz
I tried several mirrors Uwe Ligges pisze: Jarek Jasiewicz wrote: since 2.9.0 version I have a problem with installing packages: install.packages(sp) --- Please select a CRAN mirror for use in this session --- Loading Tcl/Tk interface ... done Warning: unable to access index for repository

Re: [R] 3 questions regarding matrix copy/shuffle/compares

2009-04-26 Thread David Winsemius
On Apr 26, 2009, at 9:43 AM, Esmail wrote: David Winsemius wrote: Yes. As I said before I am going to refrain from posting speculation until you provide valid R code that will create an object that can be the subject of operations. The code I have provided works, here is a run that may

Re: [R] Problem installing packages

2009-04-26 Thread Jarek Jasiewicz
I tried several mirrors But, what may be more important. This error was mentioned earlier on VISTA and WindowsXP, I use Ubuntu 8.04 Uwe Ligges pisze: Jarek Jasiewicz wrote: since 2.9.0 version I have a problem with installing packages: install.packages(sp) --- Please select a CRAN mirror

[R] Memory issues in R

2009-04-26 Thread Neotropical bat risk assessments
How do people deal with R and memory issues? I have tried using gc() to see how much memory is used at each step. Scanned Crawley R-Book and all other R books I have available and the FAQ on-line but no help really found. Running WinXP Pro (32 bit) with 4 GB RAM. One SATA drive

Re: [R] 3 questions regarding matrix copy/shuffle/compares

2009-04-26 Thread hadley wickham
I want to (1) create a deep copy of pop, I have already said *I* do not know how to create a deep copy in R. Creating a deep copy is easy, because all copies are deep copies. You need to try very hard to create a reference in R. Hadley -- http://had.co.nz/

Re: [R] 3 questions regarding matrix copy/shuffle/compares

2009-04-26 Thread David Winsemius
My understanding of the OP's request was for some sort of copy which did change when entries in the original were changed; the sort of behavior that might be seen in a spreadsheet that had a copy by reference. On Apr 26, 2009, at 11:28 AM, hadley wickham wrote: I want to (1) create a

Re: [R] 3 questions regarding matrix copy/shuffle/compares

2009-04-26 Thread Esmail
David, Good news! It seems that R has deep copy by default. I ran this simplified test and it seems I can change 'pop' without changing the saved version. POP_SIZE = 4 LEN = 8 pop=create_pop_2(POP_SIZE, LEN) cat('printing original pop\n') print(pop) keep_pop = pop pop[1,1] = 99 cat('printing

Re: [R] 3 questions regarding matrix copy/shuffle/compares

2009-04-26 Thread Esmail
hadley wickham wrote: I want to (1) create a deep copy of pop, I have already said *I* do not know how to create a deep copy in R. Creating a deep copy is easy, because all copies are deep copies. You need to try very hard to create a reference in R. Hi Hadley Right you are .. I discovered

Re: [R] 3 questions regarding matrix copy/shuffle/compares

2009-04-26 Thread hadley wickham
In that case, you would want a shallow copy, and you'd need to jump through a lot of hoops to do that in R. Hadley On Sun, Apr 26, 2009 at 10:35 AM, David Winsemius dwinsem...@comcast.net wrote: My understanding of the OP's request was for some sort of copy which did change when entries in the

Re: [R] 3 questions regarding matrix copy/shuffle/compares

2009-04-26 Thread Esmail
David Winsemius dwinsem...@comcast.net wrote: My understanding of the OP's request was for some sort of copy which did change when entries in the original were changed; the sort of behavior that might be seen in a spreadsheet that had a copy by reference. You misunderstood (my phrasing

Re: [R] Stochastic Gradient Ascent for logistic regression

2009-04-26 Thread Ravi Varadhan
Hi Tim, There are two main problems with your implementation of Stochastic gradient algorithm: 1. You are only implementing one cycle of the algorithm, i.e. it cycles over each data point only once. You need to do this several time, until convergence of parameters is obtained. 2.

Re: [R] Problem installing packages

2009-04-26 Thread Uwe Ligges
Jarek Jasiewicz wrote: I tried several mirrors But, what may be more important. This error was mentioned earlier on VISTA and WindowsXP, I use Ubuntu 8.04 Unfortunately, I cannot reproduce this for sp on Windows XP. Uwe Uwe Ligges pisze: Jarek Jasiewicz wrote: since 2.9.0 version I

Re: [R] Memory issues in R

2009-04-26 Thread David Winsemius
On Apr 26, 2009, at 11:20 AM, Neotropical bat risk assessments wrote: How do people deal with R and memory issues? They should read the R-FAQ and the Windows FAQ as you say you have.

Re: [R] Memory issues in R

2009-04-26 Thread Stefan Grosse
On Sun, 26 Apr 2009 09:20:12 -0600 Neotropical bat risk assessments neotropical.b...@gmail.com wrote: NBRA NBRAHow do people deal with R and memory issues? NBRAI have tried using gc() to see how much memory is used at each NBRA step. Scanned Crawley R-Book and all other R books I have

[R] Matching in R

2009-04-26 Thread dirk567
Dear R users, I am trying to do exact matching on a large dataset (500.000 obs), about equal size of treatment and controll group, with replacement: As for the moment I use the Match function of the Matching library. I match on 2 covariates and all observations in the treatment group have at

Re: [R] help with plotting results of lda

2009-04-26 Thread Prof Brian Ripley
On Sun, 26 Apr 2009, Uwe Ligges wrote: Works for me with library(MASS) plot(lda(Species~., data=iris)) hence you may want to profide the data to enable us to reproduce your problem... He is trying to plot the results from a cross-validation. As the help page clearly states, that is a

[R] Install packages not working in latest version?

2009-04-26 Thread Neotropical bat risk assessments
Seems that the latest version R version 2.9.0 (2009-04-17) has a glitch and will not install packages. Issue with unzipping? Works fine with R version 2.8.1 (2008-12-22) install.packages(ff) --- Please select a CRAN mirror for use in this session --- trying URL

Re: [R] Install packages not working in latest version?

2009-04-26 Thread Uwe Ligges
Neotropical bat risk assessments wrote: Seems that the latest version R version 2.9.0 (2009-04-17) has a glitch and will not install packages. Issue with unzipping? Works fine with R version 2.8.1 (2008-12-22) install.packages(ff) --- Please select a CRAN mirror for use in this session

[R] comparing matrices

2009-04-26 Thread Esmail
I'm trying to compare two matrices made up of bits. doing a simple comparison of matA == matB yields this sort of output. [,1] [,2] [,3] [,4] [,5] [,6] [1,] FALSE TRUE FALSE TRUE TRUE FALSE [2,] TRUE TRUE TRUE TRUE TRUE TRUE [3,] FALSE TRUE FALSE FALSE FALSE TRUE

[R] figure layout

2009-04-26 Thread hesicaia
Hello, I have a specific question regarding figure layout. I am tryng to make a 2 by 1 figure but I would like to make the bottom figure slightly larger than the top figure. I have read through the help posts and have tried to use fig=c(),new=T and have also tried to use split.screen and

[R] doubt in vglm output

2009-04-26 Thread priyabrata panigrahi
i am using vglm for multiple logistic regression. i have 1 response variable (total 4 category) and 5 predictor. Call: vglm(formula = class ~ PC1 + PC2 + PC3 + PC4 + PC5, family = multinomial(), na.action = na.pass) Coefficients: (Intercept):1 (Intercept):2 PC1:1 PC1:2

[R] Problem with create a tree

2009-04-26 Thread Grześ
I have a small database (file csv): 'District';'HouseType';'Income';'PreviousCustomer';'Outcome' 'Suburban';'Detached';'High';'No';'Nothing' 'Suburban';'Detached';'High';'Yes';'Nothing' 'Rural';'Detached';'High';'No';'Responded' ... itd. After instruction str() as a result I get:

Re: [R] ANOVA/statistics question

2009-04-26 Thread drmh
Hello again, In my situation, I have three variables: pretest, posttest, and cohesion. I want to work out the correlation between postest and cohesion. I looked at multiple sets of data and created ANOVA tables of them. However, as pretest and postest are sometimes correlated (with a

Re: [R] Memory issues in R

2009-04-26 Thread David Winsemius
Then post the material that would make sense for Windows. What _does_ memory.limits() return? This _was_ asked and you did not answer. How many other objects do you have in your workspace? How big are they? Jim Holtman offered this function that displays memory occupation by object and

Re: [R] comparing matrices

2009-04-26 Thread ONKELINX, Thierry
Have a look at all.equal matA - matrix(1:4, ncol = 2) matB - matA all.equal(matA, matB) matB[1,1] - -10 all.equal(matA, matB) HTH, Thierry ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research

Re: [R] THE EQUIVALENT OF SQL INNER TABLE JOIN IN R

2009-04-26 Thread Wacek Kusnierczyk
Peter Dalgaard wrote: Nigel Birney wrote: Hello all, Apologize for the newbie question. What's the easiest way to do a SQL inner table join in R? Say I have a table containing column names A, B, C and another which has columns named C, D, E. I would like to do an inner table join on C and

Re: [R] comparing matrices

2009-04-26 Thread baptiste auguie
I'm not sure I'm following you but have you tried, identical(matrix(c(1,1,1,1),ncol=2), matrix(c(1,1,1,1),ncol=2)) ?all.equal ?isTRUE ?identical and possibly the compare package, compare(matrix(c(1,1,1,1),ncol=2),matrix(c(1,1,1,1),ncol=2)) HTH, baptiste On 26 Apr 2009, at 18:02, Esmail

Re: [R] Problem installing packages

2009-04-26 Thread Jarek Jasiewicz
unfortunately this problem is difficult do reproduce. If no, it would be mentioned and removed earlier. It probably heppens in some specific cirumstances. I I look for that Jarek Uwe Ligges pisze: Jarek Jasiewicz wrote: I tried several mirrors But, what may be more important. This error

Re: [R] comparing matrices

2009-04-26 Thread David Winsemius
On Apr 26, 2009, at 1:02 PM, Esmail wrote: I'm trying to compare two matrices made up of bits. doing a simple comparison of matA == matB identical( matrix((1:4), ncol=2), matrix((1:4), nrow=2)) [1] TRUE identical( matrix((1:4), ncol=2), matrix((2:5), nrow=2)) [1] FALSE yields

Re: [R] comparing matrices

2009-04-26 Thread Esmail
baptiste auguie wrote: I'm not sure I'm following you but have you tried, identical(matrix(c(1,1,1,1),ncol=2), matrix(c(1,1,1,1),ncol=2)) ?all.equal ?isTRUE ?identical and possibly the compare package, compare(matrix(c(1,1,1,1),ncol=2),matrix(c(1,1,1,1),ncol=2)) HTH, baptiste Hi

Re: [R] comparing matrices

2009-04-26 Thread Esmail
ONKELINX, Thierry wrote: Have a look at all.equal matA - matrix(1:4, ncol = 2) matB - matA all.equal(matA, matB) matB[1,1] - -10 all.equal(matA, matB) Hi Thierry, Thanks, all.equal does indicate if it's all equal so that works great! Much nicer than my hack - thanks, Esmail

Re: [R] figure layout

2009-04-26 Thread Cuvelier Etienne
hesicaia a écrit : Hello, I have a specific question regarding figure layout. I am tryng to make a 2 by 1 figure but I would like to make the bottom figure slightly larger than the top figure. I have read through the help posts and have tried to use fig=c(),new=T and have also tried to

[R] eager to learn how to use sapply, lapply, ...

2009-04-26 Thread mauede
After a year my R programming style is still very C like. I am still writing a lot of for loops and finding it difficult to recognize where, in place of loops, I could just do the same with one line of code, using sapply, lapply, or the like. On-line examples for such high level function do not

Re: [R] eager to learn how to use sapply, lapply, ...

2009-04-26 Thread jim holtman
I think you can replace your 'for' loop with vectorized operations: if(NW 0){ Rnorm - Rnorm + sum(abs(Y[Nstart:(Nstart + NW - 1)]) ^ 2) } On Sun, Apr 26, 2009 at 1:42 PM, mau...@alice.it wrote: After a year my R programming style is still very C like. I am still writing a lot of

Re: [R] ANOVA/statistics question

2009-04-26 Thread Peter Flom
drmh douglasrmhol...@googlemail.com wrote Hello again, In my situation, I have three variables: pretest, posttest, and cohesion. I want to work out the correlation between postest and cohesion. cor(cohesion, posttest) gives you this. I looked at multiple sets of data and created ANOVA

Re: [R] eager to learn how to use sapply, lapply, ...

2009-04-26 Thread Gabor Grothendieck
Also you don't need the abs since you are squaring it anyways and seq with length argument is a bit cleaner: ix - seq(Nstart, length = NW) sum(Y[ix]^ 2) Also please read the last line on every message to r-help. The code in your post lacks all 4 of the asked for ingredients: 1. there are no

Re: [R] dotplot: labeling coordinates for each point

2009-04-26 Thread Deepayan Sarkar
On 4/26/09, Qifei Zhu zhu_qi...@yahoo.com.sg wrote: Hi David, Thanks! It looks much better now. but is there any way to add (x,y) coordinates as labels to all the points in the graph? Best case if I can enforce some conditions saying if (y10,000) label, else no label. Any advice is

[R] Generate ramified structures

2009-04-26 Thread Talita Perciano
Hello, I would like to generate ramified structures like plant root systems, river networks or trees and save the generated structure as an image. Does anyone knows if there is a way to do that with R? Thank you in advance, Talita [[alternative HTML version deleted]]

[R] Rmpi failing to install with all latest MPI packages and config arguments

2009-04-26 Thread Vince Fulco
On FC10 with openmpi and mpich2 installed, the command R CMD INSTALL Rmpi_0.5-7.tar.gz --configure-args=--with-mpi=/usr/lib64/openmpi (or /usr/lib64/mpich2) fails with the error ''cannot find mpi.h. Doing a (s)locate indicates no header file labeled as such. Would appreciate any trailheads. TIA,

[R] RWeka: How to access AttributeEvaluators

2009-04-26 Thread Michael Olschimke
Hi, I'm trying to use Information Gain for feature selection. There is a InfoGain implementation in Weka: *weka.attributeSelection.InfoGainAttributeEval* Is it possible to use this function with RWeka? If yes how? list_Weka_interfaces doesn't show it and there is no make function for

Re: [R] simulate arima model

2009-04-26 Thread Rolf Turner
On 26/04/2009, at 3:56 PM, Rebecca1117 wrote: I am new in R. I can simulate Arma, using Arima.sim However, I want to simulate an Arima Model. Say (1-B)Zt=5+(1-B)at. I do not know how to deal with 5 in this model. Can any one could help me? Thank you very much! If this is a homework

[R] lme - nlminb problem, convergence error code = 1

2009-04-26 Thread Katrin Wolf
Hi there, I have one problem calculating an linear mixed model. I have a repeated measurement project for several gases. These were measured in 3 different altitudes (sites) and three different positions within these sites (plot). each I now wanted to find out if the gasfluxes rather depend

Re: [R] Question of Quantile Regression for Longitudinal Data

2009-04-26 Thread roger koenker
I was trying to resist responding to this question since the original questioner had already been admonished twice last october about asking questions on R-help about posted code that was not only not a part of R-base, but not even a part of an R package. But the quoted comment about Stata is

[R] R 64-bit for Ubuntu 9.04 64-bit

2009-04-26 Thread Tom La Bone
I just installed Ubuntu 9.04 but there does not seem to be repository for binaries for this version. Are there going to be such repositories set up in the near future? Tom -- View this message in context: http://www.nabble.com/R-64-bit-for-Ubuntu-9.04-64-bit-tp23246229p23246229.html Sent

[R] how to install R really *locally*?

2009-04-26 Thread Oliver Kullmann
Hello, my first attempt at installing version 2.9.0 failed because I got an error Error in library(pspline) : there is no package called 'pspline' Later I realised that this comes from HOME/.RProfil, and removing that files solves that problem. However, I'm actually glad that this error

Re: [R] How to get rid of loop?

2009-04-26 Thread Peter Alspach
Tena koe Ken Would something along the following lines do what you require: set.seed(1) x - runif(100) y - rep(NA, length(x)) y[x0.25] - -1 y[x0.75] - 1 y[-1][y[-length(y)]%in%1 (x[-1]=0.25 x[-1]0.5)] - 0 # etc HTH Peter Alspach -Original Message- From:

Re: [R] R 64-bit for Ubuntu 9.04 64-bit

2009-04-26 Thread Dirk Eddelbuettel
On 26 April 2009 at 13:25, Tom La Bone wrote: | I just installed Ubuntu 9.04 but there does not seem to be repository for | binaries for this version. Are there going to be such repositories set up in | the near future? As I understand it Michael and Vincent are working on it right now. This

Re: [R] Rmpi failing to install with all latest MPI packages and config arguments

2009-04-26 Thread Dirk Eddelbuettel
On 26 April 2009 at 14:50, Vince Fulco wrote: | On FC10 with openmpi and mpich2 installed, the command R CMD INSTALL | Rmpi_0.5-7.tar.gz --configure-args=--with-mpi=/usr/lib64/openmpi (or | /usr/lib64/mpich2) fails with the error ''cannot find mpi.h. Doing a | (s)locate indicates no header file

Re: [R] Rmpi failing to install with all latest MPI packages and config arguments

2009-04-26 Thread Vince Fulco
Agreed on the redirect to SIG. mpich2-devel and -libs are installed. V. On Sun, Apr 26, 2009 at 4:51 PM, Dirk Eddelbuettel e...@debian.org wrote: On 26 April 2009 at 14:50, Vince Fulco wrote: | On FC10 with openmpi and mpich2 installed, the command R CMD INSTALL | Rmpi_0.5-7.tar.gz

Re: [R] function returns R object with name based on input

2009-04-26 Thread Jennifer Brea
Thanks for your replies. I ended up using the following: df = data.frame(year =c(1991,1991,1992,1992,1993,1993,1992,1991),x=rnorm(8),y=rnorm(8)) df year x y 1 1991 0.5565083 -1.31364232 2 1991 0.1686598 -0.20344656 3 1992 -0.1010090 -0.65681852 4 1992 0.6130324

[R] Bumps chart in R

2009-04-26 Thread Andreas Christoffersen
Hi there, I would like to make a 'bumps chart' like the ones described e.g. here: http://junkcharts.typepad.com/junk_charts/bumps_chart/ Purpose: I'd like to plot the proportion of people in select countries living for less then one USD pr day in 1994 and 2004 respectively. I have already

Re: [R] eager to learn how to use sapply, lapply, ...

2009-04-26 Thread hadley wickham
Have a look at the plyr package and associated documentation - http://had.co.nz/plyr Hadley On Sun, Apr 26, 2009 at 12:42 PM, mau...@alice.it wrote: After a year my R programming style is still very C like. I am still writing a lot of for loops and finding it difficult to recognize where,

Re: [R] Bumps chart in R

2009-04-26 Thread hadley wickham
In statistics, a bumps chart is more commonly called a parallel coordinates plot. Hadley On Sun, Apr 26, 2009 at 5:45 PM, Andreas Christoffersen achristoffer...@gmail.com wrote: Hi there, I would like to make a 'bumps chart' like the ones described e.g. here:

Re: [R] Bumps chart in R

2009-04-26 Thread Gabor Grothendieck
Have a look at plotweb in the bipartite package. On Sun, Apr 26, 2009 at 6:45 PM, Andreas Christoffersen achristoffer...@gmail.com wrote: Hi there, I would like to make a 'bumps chart' like the ones described e.g. here: http://junkcharts.typepad.com/junk_charts/bumps_chart/ Purpose: I'd

[R] RODBC - XLSX files - dropping/clearing sheets

2009-04-26 Thread Daniel Bradley
Hi! I'm manipulating XLSX data using RODBC, however a limitation which appears to be driver based is that you can't clear or drop sheets from the XLSX files, as per the following example: library(RODBC) xlsx-odbcDriverConnect(DRIVER=Microsoft Excel Driver (*.xls, *.xlsx, *.xlsm,

Re: [R] How to get rid of loop?

2009-04-26 Thread Ken-JP
The code below shows what I'm trying to get rid of. If there is no way to get rid of the loop, I will try to use package( inline ). I'm just curious as to whether there is a vector way of doing this algorithm. #

Re: [R] Bumps chart in R

2009-04-26 Thread Mike Lawrence
Here's a ggplot2 based solution: #load the ggplot2 library library(ggplot2) #here's the data provided by Andreas countries - c(U-lande, Afrika syd for sahara, Europa og Centralasien, Lantinamerika og Caribien,Mellemøstenog Nordafrika,Sydasien,ØStasien og stillehaveet, Kina, Brasilien) poor_1990

[R] How to create a graph layout?

2009-04-26 Thread Christian Bustamante
I all, I want to create a graph layout in a 3x3 matrix like this: ylab |__| |__| |__| ___ ___ ___ ylab |__| |__| |__| ___ ___ ___ ylab |__| |__| |__| xl xl xl With this layout, then I'll insert the 9 plots. How ca I create it? -- CdeB

[R] Plotting polynomial fit

2009-04-26 Thread Ronnen Levinson
Hi. Is there an analog to abline() that can be used to plot a polynomial fit? For example, I can draw the straight-line fit fit - lm(y ~ x) via abline(coef=fit$coef) but I'm not sure how to draw the polynomial fit fit - lm(y ~ poly(x,2)) I do see the function curve(), but not

Re: [R] RODBC - XLSX files - dropping/clearing sheets

2009-04-26 Thread Gabor Grothendieck
Try using RDCOMClient or rcom: library(RDCOMClient) xl - COMCreate(Excel.Application) # next line optional xl[[Visible]] - TRUE wb - xl[[Workbooks]]$Open(/mydir/sample.xlsx) sheet - wb$Sheets(Sheet2) xl[[DisplayAlerts]] - FALSE

  1   2   >