[R] Read and show Bitmap Images

2012-01-15 Thread Alaios
Dear all, I am looking for a function that can plot bitmap images and by plotting I mean a function that can read an image's matrix structure with integers and assign colors. Do you please suggest me what else I can do for plotting these images? I would like to thank you for your reply B.R Alex

[R] rho stat from a fasta sequence file

2012-01-15 Thread utpal
Hi all, I have a sequence file (fasta format) and want to calculate the rho statistics for dinucleotide abundance value on my data.. the code which I use is (using seqinr library and current working directory) seq_info<-read.fasta("gene.txt") rho(seq_info[1],2) but it yields only the dinucleoti

Re: [R] problem with table.CAPM in PerformanceAnalytics

2012-01-15 Thread Dom Pazzula
Pete, That seems to work.  I wonder what the problem was. Thanks, Dom   Dominic J. Pazzula +-+-+-+-+-+-+-+-+-+-+- domp...@yahoo.com From: Pete Brecknock To: r-help@r-project.org Sent: Sunday, January 15, 2012 8:13 PM Subject: Re: [R] problem with table.CAPM in

Re: [R] A recent Post

2012-01-15 Thread Roy Mendelssohn
I understand your point., and is why I have hesitated to say anything in the past My point is the following last line of that part of the guide: > Ad hominem comments are absolutely out of place I guess with a Bayesian approach you can give "absolutely" a fuzzy meaning. Why I am not a Bayes

Re: [R] A recent Post

2012-01-15 Thread Joshua Wiley
Hi Roy, This still needs refinement and may benefit from moving to a Bayesian approach, but this is roughly my approach: status <- function(help, harshness, OP, contributions, prior, duration) { params <- c(help, harshness, OP, contributions, prior, duration) stopifnot(all(params > 0 & param

Re: [R] ggplot2 stacked bar - sum of values rather than count

2012-01-15 Thread J Toll
On Sun, Jan 15, 2012 at 9:15 PM, Richard M. Heiberger wrote: > This is based on lattice, not ggplot, and I hope I understand what you are > trying to do. > > require(lattice) > > x <- read.table(text=" >     Value Bar Segment >  1.10020075   1   1 > -1.37734577   2   1 >  2.50702876  

Re: [R] A recent Post

2012-01-15 Thread Roy Mendelssohn
>From the r-help posting guidelines: > Good manners: Remember that customs differ. Some people are very direct. > Others surround everything they say with hedges and apologies. Be tolerant. > Rudeness is never warranted, but sometimes `read the manual' is the > appropriate response. Don't waste

Re: [R] ggplot2 stacked bar - sum of values rather than count

2012-01-15 Thread Richard M. Heiberger
This is based on lattice, not ggplot, and I hope I understand what you are trying to do. require(lattice) x <- read.table(text=" Value Bar Segment 1.10020075 1 1 -1.37734577 2 1 2.50702876 3 1 0.58737028 3 2 0.21106851 3 3 -2.50119261 4

Re: [R] problem with table.CAPM in PerformanceAnalytics

2012-01-15 Thread Pete Brecknock
Dom Pazzula wrote > > All, > I'm attempting to run this: > table.CAPM(series[,"Strat.Return",drop=FALSE],series[,"spy.Return",drop=FALSE]) > > > and getting this error > Error in as.vector(data[, i]) : subscript out of bounds > > > I've searched around and cannot find a solution to the proble

Re: [R] rpart question

2012-01-15 Thread Joshua Wiley
Hi Amanda, Sorry for the bit of a slow response (classes and research have been chaotic). Below are details on what I looked at and a few suggestions at the end for what you can do. To the general R community: summary.rpart() makes explicit the default dropping behavior of `[` which makes me thi

[R] ggplot2 stacked bar - sum of values rather than count

2012-01-15 Thread J Toll
Hi, I'm trying to create a stacked bar plot using ggplot2. Rather than plotting the count of each of the 13 "Bar" factors on the Y axis, I would like to represent the sum of the Values associated with each of the 13 "Bar" factors. Is there a way to do that? Given the following data, that would ob

Re: [R] howto test a package without installation

2012-01-15 Thread Ken Hutchison
Did not fully read the without installing it part. Mea Culpa, Ken On Jan 15, 2012, at 8:56 PM, Jeff Newmiller wrote: > I don't believe you can. However, you need not install it into a system-wide > library directory... your personal library (e.g. > /home/jonas/R/x86_64-pc-linux-gnu-li

Re: [R] howto test a package without installation

2012-01-15 Thread Jeff Newmiller
I don't believe you can. However, you need not install it into a system-wide library directory... your personal library (e.g. /home/jonas/R/x86_64-pc-linux-gnu-library/2.14) should be sufficient. --- Jeff Newmiller

[R] Fwd: Trouble installing packages on R2.14.1

2012-01-15 Thread Ken Hutchison
Begin forwarded message: > From: Ken Hutchison > Date: January 15, 2012 8:54:49 PM EST > To: Ben Bolker > Subject: Re: [R] Trouble installing packages on R2.14.1 > > Check browser proxy settings and run R.exe with the proper flags to use them > from cmd. > Hope that helps, > Ken > > >

Re: [R] howto test a package without installation

2012-01-15 Thread Ken Hutchison
install.packages('pathtotargz', repos=FALSE) I believe that will get you going. Hope that's helpful, Ken On Jan 15, 2012, at 8:33 PM, Joshua Wiley wrote: > Hi Jonas, > > Look at Hadley Wickham's devtools package. It is designed with this > sort of thing. That said, it really is not too

Re: [R] Trouble installing packages on R2.14.1

2012-01-15 Thread Ben Bolker
Crawford Mark derm.qld.gov.au> writes: > > Hi, > Have just recently downloaded R version 2.14.1 and can't install > packages (Code below), 32bit version on windows. From the looks of > things access is the problem but I don't know where to start, would it > be my system that is the problem or i

Re: [R] Need help interpreting the logit regression function

2012-01-15 Thread cberry
Rohit Pandey writes: > Hello R community, > > I have a question about the logistic regression function. > Specifically, when the predictor variable has not just 0's and 1's, > but also fractional values (between zero and one). I get a warning > when I use the "glm(formula = ... , family = binomia

Re: [R] R-help Digest, Vol 107, Issue 14

2012-01-15 Thread Paul Gilbert
On 12-01-14 06:00 AM, r-help-requ...@r-project.org wrote: Date: Fri, 13 Jan 2012 02:09:09 -0800 (PST) From: statquant2 To:r-help@r-project.org Subject: Re: [R] simulating stable VAR process Message-ID:<1326449349804-4291835.p...@n4.nabble.com> Content-Type: text/plain; charset=us-ascii Hello Pau

[R] Trouble installing packages on R2.14.1

2012-01-15 Thread Crawford Mark
Hi, Have just recently downloaded R version 2.14.1 and can't install packages (Code below), 32bit version on windows. From the looks of things access is the problem but I don't know where to start, would it be my system that is the problem or is the repository not set up for the latest version yet

Re: [R] Troubles with stemming (tm + Snowball packages) under MacOS

2012-01-15 Thread Julien Velcin
I use the version 1.6: $ java -version java version "1.6.0_26" Julien On Jan 15, 2012, at 8:55 PM, Milan Bouchet-Valat wrote: Le dimanche 15 janvier 2012 à 16:32 +0100, Julien Velcin a écrit : Unfortunately, it doesn't work. I've installed JGR and launched my script. I still obtain an error:

[R] problem with table.CAPM in PerformanceAnalytics

2012-01-15 Thread Dom Pazzula
All, I'm attempting to run this: table.CAPM(series[,"Strat.Return",drop=FALSE],series[,"spy.Return",drop=FALSE]) and getting this error Error in as.vector(data[, i]) : subscript out of bounds I've searched around and cannot find a solution to the problem.  I've used this in the past without pr

Re: [R] howto test a package without installation

2012-01-15 Thread Joshua Wiley
Hi Jonas, Look at Hadley Wickham's devtools package. It is designed with this sort of thing. That said, it really is not too difficult to install as long as you have a working tool chain (which you will need to test it anyway). R CMD INSTALL /tmp/sitools R require(sitools) and away you go testi

[R] howto test a package without installation

2012-01-15 Thread Jonas Stein
Hi, how can i play around with my first selfwritten package [*] without to install it to my debian system? I think of something like doing this: /tmp/$ R R version 2.11.1 (2010-05-31) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 > library(/tmp/sitools) 3 * k

Re: [R] Calling Fortran from R: Issues with dynamic loading of fortran dll

2012-01-15 Thread Duncan Murdoch
On 12-01-15 5:34 PM, David Stevens wrote: I successfully used .Fortran to load and execute my fortran procedures under WinXP and 32 bit R. Alas, the same isn't true with my next Windows 7/64 machine, R 2.14.1 (64 bit) and the gnu gfortran (64) compiler (mingw64 v. 4.6.1). Though I'm able to compi

[R] Need help interpreting the logit regression function

2012-01-15 Thread Rohit Pandey
Hello R community, I have a question about the logistic regression function. Specifically, when the predictor variable has not just 0's and 1's, but also fractional values (between zero and one). I get a warning when I use the "glm(formula = ... , family = binomial(link = "logit"))" which says: "

Re: [R] Union/Intersect two continuous sets

2012-01-15 Thread William Dunlap
Late last month Hans Borchers suggested the following solution (which I trivially converted to a function): u0 <- function(M) { # Mnew the union of intervals in M, an (nx2)-matrix with n > 1: o <- order(M[, 1], M[, 2]) L <- M[o, 1]; R <- M[o, 2] k <- 1 Mnew <- matrix(c(L[k], R[

Re: [R] ggExtra package installation fails

2012-01-15 Thread Tengfei Yin
Hi Milan, thanks for the response, I guess grid.arrange does't really serve the purpose I want(or maybe my grid using level just limits the usage..), like mentioned by Baptiste. On Sun, Jan 15, 2012 at 5:13 PM, baptiste auguie < baptiste.aug...@googlemail.com> wrote: > ggExtra was not compatible

Re: [R] ggExtra package installation fails

2012-01-15 Thread baptiste auguie
ggExtra was not compatible with the recent changes in ggplot2 and had become partially redundant, so I removed it. grid.arrange and arrangeGrob in gridExtra won't help with the alignment of axes. align.plots, which I can always send you offlist if you want, extracted the size of axes and legends f

Re: [R] deviance and variance - GAM models

2012-01-15 Thread Simon Wood
On 01/13/2012 08:25 PM, collifu wrote: Hi all, This is pretty basic but I am not an expert and I couldn't find anything in the forum or my statistics book about it. I was reading a paper and the authors were using both "explained deviance" and "explained variance" as synonyms. They were describi

[R] Calling Fortran from R: Issues with dynamic loading of fortran dll

2012-01-15 Thread David Stevens
I successfully used .Fortran to load and execute my fortran procedures under WinXP and 32 bit R. Alas, the same isn't true with my next Windows 7/64 machine, R 2.14.1 (64 bit) and the gnu gfortran (64) compiler (mingw64 v. 4.6.1). Though I'm able to compile the routines from the command line us

Re: [R] general question on Spotfire

2012-01-15 Thread Frank Harrell
Thanks for your note Uwe. Yes I think a lot of the needed work was related to implementing R functions that many of us use that are not available in S-Plus, plus what to do about plotmath. It wasn't enough to just be able to load the R package. I don't think implementation of the needed R functi

Re: [R] ggExtra package installation fails

2012-01-15 Thread Milan Bouchet-Valat
Le dimanche 15 janvier 2012 à 15:43 -0600, Tengfei Yin a écrit : > Hi dear all, > > I cannot find ggExtra source code or install it by > > install.packages("ggExtra", repos="http://R-Forge.R-project.org";) This package has been removed, it's been marked as deprecated by its author: http://groups.

[R] ggExtra package installation fails

2012-01-15 Thread Tengfei Yin
Hi dear all, I cannot find ggExtra source code or install it by install.packages("ggExtra", repos="http://R-Forge.R-project.org";) I am really interested in one of the functions inside called "align.plots", but cannot find the package on-line. Do I miss something? or is there any function simil

Re: [R] More issues with apply

2012-01-15 Thread LCOG1
Sorry I wasn't completely clear, I am not 100% sure I know all the terminology but what I want the array to store is both 'Taz' with the 'Bldgtype' count designations like the result from data set one returns. When there isn't more than one type of 'Bldgtype' like in data set 2 the 'Bldgtype' desi

Re: [R] Regression for Rank Order Dependent Variable

2012-01-15 Thread Iasonas Lamprianou
Dear colleagues, I am looking for a package to analyze rank-order data in the form of regression. My outcome (dependent) variable is a rank ordering in the sense that the respondents to a  questionnaire were asked a single question: to assign the values 1 to 4 to four different objects (e.g. t

Re: [R] parLapply within a function

2012-01-15 Thread Uwe Ligges
Katja, since the information given in this message is not sufficient to be able to help (question whether you are forking or parallelization happens through, e.g., MPI or sockets) and no reproducible code is given yet (as the posting guide ask you to provide), let us discuss this off-list and

Re: [R] add column with values found in another data frame

2012-01-15 Thread jdog76
Thanks so much Jorge and Pete ! Pete Brecknock wrote > > > frame2$age = frame1[match(frame2$ID, frame1$ID),"age"] > > This match function was exactly what I needed to know! -- View this message in context: http://r.789695.n4.nabble.com/add-column-with-values-found-in-another-data-frame-tp4

Re: [R] problem in R related to import n export

2012-01-15 Thread jim holtman
My quick guess is that it not an R problem, but with your version of EXCEL. What version are you using, since this used to be a limitation of Excel. Just tried to output an object with 1436 columns and read it fine into 2007 version of EXCEL. On Sun, Jan 15, 2012 at 10:09 AM, mukul purva wrote:

Re: [R] Nabble? Was Re: function to replace values doesn't work on vectors

2012-01-15 Thread Uwe Ligges
On 14.01.2012 20:41, Duncan Murdoch wrote: On 12-01-13 11:48 AM, Sarah Goslee wrote: ... I hope that it was a momentary glitch; greater disagreement between Nabble and the email list will cause all sorts of fun. If the interface, whatever it is, starts stripping out code? I'll have to quit ans

Re: [R] general question on Spotfire

2012-01-15 Thread Uwe Ligges
On 12.01.2012 17:38, Frank Harrell wrote: As a slight aside, Tibco/Spotfire originally planned to provide a capability to load R packages into S-Plus. This always seemed to me to be a hard thing to do, and if my understanding is correct, this proved to be too difficult to do in S-Plus, at lea

Re: [R] Word Wrap

2012-01-15 Thread Martin Morgan
On 01/15/2012 09:10 AM, Barry Rowlingson wrote: On Sun, Jan 15, 2012 at 4:40 PM, Francisco wrote: Hello, I have to write a big sentence with cat() and I would like that R automatically adds a new line when it is needed (when the text arrives at the end of the window), the same as Windows Notepa

Re: [R] Word Wrap

2012-01-15 Thread Florent D.
see the 'fill' argument for cat(). On Sun, Jan 15, 2012 at 12:10 PM, Barry Rowlingson wrote: > On Sun, Jan 15, 2012 at 4:40 PM, Francisco > wrote: >> Hello, >> I have to write a big sentence with cat() and I would like that R >> automatically adds a new line when it is needed (when the text arr

Re: [R] Word Wrap

2012-01-15 Thread Barry Rowlingson
On Sun, Jan 15, 2012 at 4:40 PM, Francisco wrote: > Hello, > I have to write a big sentence with cat() and I would like that R > automatically adds a new line when it is needed (when the text arrives at > the end of the window), the same as Windows Notepad does (Word Wrap). > How could I do? A li

Re: [R] par.plot() for repeated measurements

2012-01-15 Thread JerryLee
Dear Jonas, Thanks for you your suggestion. I tried double quotes (as in following code). It still does not work. Attached is the data. library(gamlss) data=read.csv("D:\\lzg\\jointmodel\\data.csv",header=T) data$ID<-factor(data$ID) par.plot(factpal~timetodeath2,data=data,sub=ID,color=FALSE,ylim=

[R] Word Wrap

2012-01-15 Thread Francisco
Hello, I have to write a big sentence with cat() and I would like that R automatically adds a new line when it is needed (when the text arrives at the end of the window), the same as Windows Notepad does (Word Wrap). How could I do? Thank you F. __

Re: [R] How can I doing Quality adjusted survival analysis in R?

2012-01-15 Thread Pedro Mota Veiga
Thaks Ryusuke, One question, is it possible to have an increased quality adjusted survival during a period of time? I ask that because I saw a paper published in a cientific journal with a adjusted survival curve that increase during a period and dicrease after that. Best regards Pedro Mota V

[R] problem in R related to import n export

2012-01-15 Thread mukul purva
hello, i hav a prob in R lang. i want to do correlation of data frame 22810(gene) rows and 1436 colums(experiment) when i covert this file in csv or xls format it reads only 1024 cloumns. n when i do correlation of data mean 22810 *22810 matrix made in terminal then i export in csv or xls b

[R] A question about cointegration - How can we find the standard deviation in the cointegration relationship ?

2012-01-15 Thread miao
Hello, I am using urca package to run cointegration. I would like to find the standard error in the (normalized, Johansen) cointegration relationship. How can I do it? As far as I know, The function "cajorls" in the "urca" package provides the normalized cointegrating relationships. Never

Re: [R] metafor: weights computation in Mantel-Haenszel method

2012-01-15 Thread Viechtbauer Wolfgang (STAT)
Yes, I forgot to divide the weights by the sum of the weights (and multiply by 100) for weights.rma.mh(). This bug was pointed out to me a few weeks ago and this will be changed in the next version of the package. Best, Wolfgang -- Wolfgang Viechtbauer, Ph.D., Statistician Department of Psychi

[R] Reading MINE output into a matrix

2012-01-15 Thread Matt Considine
I've benefited from this list with input on how to build up a symmetrical matrix. The purpose of that query was to work with the output from the MINE routine posted at www.exploredata.net To the extent it helps others, here is the script that I was working on an which turns a given MINE outpu

Re: [R] The Future of R | API to Public Databases

2012-01-15 Thread Jason Edgecombe
At first, I thought that RDF and SDMX were two competing standards and was disheartened, but it appears that there is collaboration between them, yay! http://groups.google.com/group/publishing-statistical-data/browse_thread/thread/531b1b5a73397c1c?pli=1 On 01/15/2012 08:31 AM, Benjamin Weber w

Re: [R] Troubles with stemming (tm + Snowball packages) under MacOS

2012-01-15 Thread Milan Bouchet-Valat
Le vendredi 13 janvier 2012 à 15:49 +0100, Julien Velcin a écrit : > Dear all, > > I have some troubles using the stemming algorithm provided by the tm > (text mining) + Snowball packages. > Here is my config: > > MacOS 10.5 > R 2.12.0 / R 2.13.1 / R 2.14.1 (I have tried several versions) > >

[R] How can I doing Quality adjusted survival analysis in R?

2012-01-15 Thread Chris
You can prepare only the so-called Partitioned Survival Analysis, but those are without the quality of life weights. I describe Partitioned analysis in my paper in The American Statistician and give citations to description of the method in a  Biometrics paper. There are several variations on q

Re: [R] determining the difference between 2 character strings

2012-01-15 Thread R. Michael Weylandt
On Jan 15, 2012, at 3:28 AM, Philip Robinson wrote: > Hi, > > I am struggling, I have 2 lists with shared elements, one ~600, one > ~1000, and I need to determine the difference between them. > > They are character strings, and to use setdiff(), or unique() I need > vectors. You can have ch

Re: [R] The Future of R | API to Public Databases

2012-01-15 Thread Benjamin Weber
Yes, R-devel would be the right mailing list for this discussion. As some people pointed out, the problem definition is vague. This was to encourage people to share their *different* perceptions about the problem and to get to some extent a consensus. My starting point has come from my mind, conse

Re: [R] Quality adjusted survival analysis

2012-01-15 Thread ryusuke
paper for reference http://www.google.com.hk/url?sa=t&rct=j&q=Package+iSQoL&source=web&cd=1&ved=0CCYQFjAA&url=http%3A%2F%2Fwww.stat.sinica.edu.tw%2Fjshwang%2Fweb%2Fwp-content%2Fplugins%2Fdownload-monitor%2Fdownload.php%3Fid%3D2&ei=fcMST6CgJ4yIrAfb_bT9AQ&usg=AFQjCNG4elXcSmsAElv2M5bHoqTNf8qXGQ&cad=rj

Re: [R] Is it possible to "right align" text in R graphics?

2012-01-15 Thread Majid Einian
Hello Tall, sorry for late reply. I meant just the direction of the closing parentheses. This is the way it should be (and is on my linux system) : http://myup.ir/images/93380129091837074902.png This is the way it is (on windows system) http://myup.ir/images/40188837660458525731.png On 1/3/12, T

Re: [R] Correct Localized Numbers on Plots, related to glibc!

2012-01-15 Thread Prof Brian Ripley
You need to read ?Sys.setlocale (surely part of the homework the R posting guide asked of you). It is vital that you give the 'at a minimum' information requested in the posting guide when asking such questions, and we also absolutely need to know what graphics device you are trying to use.

Re: [R] determining the difference between 2 character string

2012-01-15 Thread Jim Holtman
why cann't you use setdiff or unique? in the sample data how do you tell the difference? do you compare the whole string, or partial? a better definition of your problem and what you have tried would help. it would seem like there is nothing the same between your sample data. Sent from my iP

[R] Correct Localized Numbers on Plots, related to glibc!

2012-01-15 Thread Majid Einian
Dear R Helpers, I want to localize my plots, i.e. the numbers by x & y axis be Persian, using Persian numerals and Persian decimal separator. I change the locale to fa_IR.utf8, but nothing on plots change. I can change the numerals shaping to Persian ones (۱۲۳۴ instead of 1234) using some non-stan

[R] determining the difference between 2 character strings

2012-01-15 Thread Philip Robinson
Hi, I am struggling, I have 2 lists with shared elements, one ~600, one ~1000, and I need to determine the difference between them. They are character strings, and to use setdiff(), or unique() I need vectors. I don't know how to force these character strings into a form where you can use functio