Re: [R] ddply from plyr package - any alternatives?

2011-08-25 Thread Tal Galili
Hi Adam, I don't think there is a faster alternative to plyr, without doing it in nested for loops, with a lot of book-keeping of variables (but if someone here were to correct me, I'd be happy to know). Two things to consider: 1) See if you can optimizing your function. (there is a lot of mater

Re: [R] Little problem with robustbase::covOKG() [?]

2011-08-25 Thread kv
super! works like a charm, best -- View this message in context: http://r.789695.n4.nabble.com/Little-problem-with-robustbase-covOKG-tp3752169p3767510.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list ht

Re: [R] ddply from plyr package - any alternatives?

2011-08-25 Thread Paul Hiemstra
Hi Adam, A recent thread on R-help deals exactly with your problem. In one of the responses I compare ddply to a number of alternative solutions (using ave and data.table) [1]. The test in the e-mail shows that for large amounts of unique categories, ddply is quite slow. Hadley (Wickham, author o

Re: [R] Regression by factor using "sapply"

2011-08-25 Thread Paul Hiemstra
I would take a look at the ddply function from the plyr package. Without a reproducible example (most importantly an example dataset) I cannot give you some example code. cheers, Paul On 08/24/2011 05:51 PM, elh wrote: > Apologies for the elementary nature of the question (yes, I'm another > new

[R] how to read a group of files into one dataset?

2011-08-25 Thread Jie TANG
for example : I have files with the name "ma01.dat","ma02.dat","ma03.dat","ma04.dat",I want to read the data in these files into one data.frame flnm<-paste("obs",101:114,"_err.dat",sep="") newdata<-read.table(flnm,skip=2) data<-(flnm,skip=2) but the data only contains data from the flnm[1] I als

[R] axes start at end of plot

2011-08-25 Thread Sumukh Sathnur
Hello all, I used the following script to generate a plot, but when I try to place axes, the axes start at the end of my plot and continue off the plot; the attached image shows what happens when I try to generate axes. library(fields) as.matrix(read.table("Matrix.txt", sep="\t"))->x #the mat

Re: [R] how to read a group of files into one dataset?

2011-08-25 Thread Eik Vettorazzi
Hi Jie, you have to merge the sequential data.frames, and depending on the structure of your inputs and the way you want your resulting data.frame (which you both didn't specify) either ?merge or ?rbind should help. cheers Am 25.08.2011 10:17, schrieb Jie TANG: > for example : I have files wit

[R] help about copula package

2011-08-25 Thread Deniz SIGIRLI
Hello all, I've got a data set with two variables. I want to learn which copula is suitable for my data set. I tried gofCopula function.I used this code:gofCopula(gumbelCopula(1.5, d=2),x, method = "itau")x is my data set, and "1.5" is my patameter which I found by converting Kendall tau to gu

[R] Odp: Boxplot orders

2011-08-25 Thread Petr PIKAL
Hi > > Hi there, > > I have length data of an organism over the year and I want to make a > boxplot. I get the boxplot just fine but the months are all out of order. In > the data set they are in order from Jan-Dec...how can I fix this problem? Probably months are either character or factor v

[R] Odp: Howto convert Linear Regression data to text

2011-08-25 Thread Petr PIKAL
> > Dear all, > > How can I covert lm data to text in the form of "y=ax+b, r2" and how do I > calculate R-squared(r2)? > > Thanks. > > Code: > x=18:29 > y=c(7.1,7,7.7,8.2,8.8,9.7,9.9,7.1,7.2,8.8,8.7,8.5) > res=lm(y~x) I am not sure what do you mean by text. There are some lm extraction func

[R] Autocorrelation using acf

2011-08-25 Thread Vincy Pyne
Dear R list As suggested by Prof Brian Ripley, I have tried to read acf literature. The main problem is I am not the statistician and hence have some problem in understanding the concepts immediately. I came across one literature (http://www.stat.nus.edu.sg/~staxyc/REG32.pdf) on auto-correlatio

Re: [R] ddply from plyr package - any alternatives?

2011-08-25 Thread AdamMarczak
Thank you for suggestions, apparently data.table is much quicker than ddply and it's fantastic to use. I forgot to mention in my topic I'm looking for alternative in R 2.10 version as on my platform our server runs older version of software which only support R up to version of R-2.10 (upgrade is

[R] Application of results from smooth.spline outside R

2011-08-25 Thread Jesper
Hi, I want to use the result from smooth.spline outside R. I take my data ,which is 180 point stored in x and y s <- smooth(x,y) I can know use to e.g. find the interpolated value at e.g. x=500 predict (s,500) My problem is, that i don't know how to implement the predict function. I have look

[R] Create two uniformly random variables correlated

2011-08-25 Thread Soberon Velez, Alexandra Pilar
Hello, I want to create two random variables (x1,x2) both with uniform distribution bounded by (-1) and (1) that has a correlation of 0.6 between them. Does somebody know how I can do it? For normal random variables I known how to implement it with the rmvnorm command but I don't know how t

[R] How to vary the distance between the pairs of a side-by-side barplot?

2011-08-25 Thread Michael Eisenring
Hello, for my project I wanted to create a side-by-side barplot. For each of 143 treatments there are always two measurement. This two measurements (a pair) are represented by a red resp. a yellow bar. I was able to write the code for the side-by-side barplot. Unfortunately however I am not able

[R] FW: How to count numbers of a vector and use them as index values?

2011-08-25 Thread Freeh Sophia
Hi Paul, I would use something like this: > x <- c(2,2,3,3,4,6) > table(x) x 2 3 4 6 2 2 1 1 > x <- factor(x, levels=1:8) > table(x) x 1 2 3 4 5 6 7 8 0 2 2 1 0 1 0 0 Sarah On Sun, Jul 31, 2011 at 5:41 PM, Paul Menzel wrote: > Dear R folks, > > > I am sorry to ask this simple question, but my

Re: [R] Change color in forest.rma (metafor)

2011-08-25 Thread Paola Tellaroli
Thank you for your attention and help! In this way I get the diamond coloured, but actually I would have the squares representing the values of the individual studies coloured. Is it somehow possible? *Paola* 2011/8/24 Viechtbauer Wolfgang (STAT) < wolfgang.viechtba...@maastrichtuniversity.nl>

Re: [R] unused argument(s) (Header = True) help!

2011-08-25 Thread shardman
Thanks everyone, that's really helpful, Cheers, Sam -- View this message in context: http://r.789695.n4.nabble.com/unused-argument-s-Header-True-help-tp3765651p3767482.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.o

[R] rpart: plot without scientific notation

2011-08-25 Thread Jay
While I'm very pleased with the results I get with rpart and rpart.plot, I would like to change the scientific notation of the dependent variable in the plots into integers. Right now all my 5 or more digit numbers are displayed using scientific notation. I managed to find this: http://tolstoy.new

Re: [R] how to read a group of files into one dataset?

2011-08-25 Thread Mohammed Ouassou
inputDataPath <- "/home/.../bla/"; #Directory containing data files szPattern <- ".dat"; # File extension # Get all files name in the specified directory file2process <- list.files(inputDataPath, pattern=szPattern); # Get number of files to be processed iFileCnt <- len

Re: [R] axes start at end of plot

2011-08-25 Thread Jim Lemon
On 08/25/2011 06:31 PM, Sumukh Sathnur wrote: Hello all, I used the following script to generate a plot, but when I try to place axes, the axes start at the end of my plot and continue off the plot; the attached image shows what happens when I try to generate axes. library(fields) as.matrix(read

[R] apply function to spatial grid data frame to calculate CTI

2011-08-25 Thread Tom Vanwalleghem
Dear R helpers, I'm new to R and struggling with applying functions to a spatial grid data frame, called data for example.data has a variable named "slope" and "uparea" Basically, I need to calculate the tan(data$slope)Next, I need to divide data$uparea by data$slope and take the log() of that a

Re: [R] Function rank() for data frames (or multiple vectors)?

2011-08-25 Thread Sebastian Bauer
Hi! On 08/24/2011 07:46 PM, David Winsemius wrote: I was looking for an elegant solution ;) In the real case I have double values and this would be quite inefficient then. Still no r-code: Then what about rank(order(...) , further-ties.method-argument) ? I think that, as order() always give

Re: [R] ddply from plyr package - any alternatives?

2011-08-25 Thread Hadley Wickham
> z <- ddply(past, c("GEO_CNTRY_NAME","PROD_SEG_NAME"), >  function(x) summary(lm(VAL~fy,x))$r.squared) > > But when ave is not exactly doing what I need. Above code runs under a > minute for my data set where as ave runs over 8 mins. It's hard to know without a reproducible example, but I doubt t

Re: [R] Autocorrelation using acf

2011-08-25 Thread Raphael Saldanha
Hi Vincy, Take a look on the material bellow, maybe they can help you: http://www.statoek.wiso.uni-goettingen.de/veranstaltungen/zeitreihen/sommer03/ts_r_intro.pdf http://www.maths.bris.ac.uk/~mazlc/TSA/r-ts.pdf http://www.stat.pitt.edu/stoffer/tsa2/R_time_series_quick_fix.htm On Thu, Aug 25,

[R] X11 problem

2011-08-25 Thread carol white
Hello, TO update R under Ubuntu on my macchine, I downloaded R-2.13.1. When running ./configure, I got the following error message: checking for X... no configure: error: --with-x=yes (default) and X11 headers/libs are not available But I could use X11() function in R-2-10. So where does the p

Re: [R] X11 problem

2011-08-25 Thread Kevin E. Thorpe
On 08/25/2011 08:27 AM, carol white wrote: Hello, TO update R under Ubuntu on my macchine, I downloaded R-2.13.1. When running ./configure, I got the following error message: checking for X... no configure: error: --with-x=yes (default) and X11 headers/libs are not available But I could use

Re: [R] Application of results from smooth.spline outside R

2011-08-25 Thread Jean V Adams
I'm assuming that you meant to write s <- smooth.spline(x, y) You can look at the code for the predict method for smooth.spline by typing getAnywhere(predict.smooth.spline) If a new value is provided (i.e., the 500 in your example) then the predict method for smooth.spline.fit is applied to t

Re: [R] How to vary the distance between the pairs of a side-by-side barplot?

2011-08-25 Thread Jean V Adams
>From the help for barplot ?barplot The argument "... space may be specified by two numbers, where the first is the space between bars in the same group, and the second the space between the groups. If not given explicitly, it defaults to c(0,1) if height is a matrix and beside is TRUE ..." S

[R] Possible Error in generic function rzigp in ZIGP Package

2011-08-25 Thread Chee Chen
Dear All, I would like to report a possible mistake in the generic function rzigp in ZIGP package, and seek for help. Yesterday I tried the "rzigp" function, in this way: 1. generate poisson means from a pareto distribuion, 2. generate dispersion for each poisson from gamma distribution 3. gene

Re: [R] Create two uniformly random variables correlated

2011-08-25 Thread Jean V Adams
Have you tried searching the R-Help archives? I found some answers to a similar question by Googling r multivariate uniform distribution Jean Soberon Velez, Alexandra Pilar wrote on 08/25/2011 04:15:11 AM: > > Hello, > > > > I want to create two random variables (x1,x2) both with unifo

Re: [R] X11 problem

2011-08-25 Thread carol white
Thanks for your reply. I used the following commands: sudo apt-get update sudo apt-get install r-base - Original Message - From: Kevin E. Thorpe To: carol white Cc: "r-h...@stat.math.ethz.ch" Sent: Thursday, August 25, 2011 2:34 PM Subject: Re: [R] X11 problem On 08/25/2011 08:27 AM, c

Re: [R] Opening package manual from within R

2011-08-25 Thread Uwe Ligges
On 24.08.2011 21:59, Tyler Rinker wrote: Apparently my request to view the help pages Well, although the questions we see on R-help suggest that some useRs do not look at those help pages, most of us do, but in text or html mode rather than the pdf version. Best, Uwe Ligges is not a

Re: [R] X11 problem

2011-08-25 Thread Kevin E. Thorpe
On 08/25/2011 08:49 AM, carol white wrote: Thanks for your reply. I used the following commands: sudo apt-get update sudo apt-get install r-base Is that how you installed R-2-10 or how you tried to install 2.13.1? - Original Message - From: Kevin E. Thorpe To: carol white Cc: "r-h...

Re: [R] how to read a group of files into one dataset?

2011-08-25 Thread Dennis Murphy
Hi: Similar in vein to the other respondents, you could try something like this: On Thu, Aug 25, 2011 at 1:17 AM, Jie TANG wrote: > for example : I have files with the name >  "ma01.dat","ma02.dat","ma03.dat","ma04.dat",I want to read the data in > these files into one data.frame > # Your file n

Re: [R] X11 problem

2011-08-25 Thread carol white
both. But it worked for R-2-10. I added  deb http://cran.ma.imperial.ac.uk/bin/linux/ubuntu lucid/ into  /etc/apt/sources.list and ran apt-get update as root and got the following error message: W: GPG error: http://cran.ma.imperial.ac.uk lucid/ Release: The following signatures couldn't

[R] Synchronizing R libraries on N machines?

2011-08-25 Thread Giovanni Petris
Hello! I am using R on two different machines (under Ubuntu and OS X, but this is probably irrelevant) and I would like to keep the two installations 'synchronized', in particular in terms of installed packages. For example, if I install package xxx on my Linux machine, I would like to find it ins

Re: [R] X11 problem

2011-08-25 Thread Alexander Engelhardt
Am 25.08.2011 15:05, schrieb carol white: both. But it worked for R-2-10. I added deb http://cran.ma.imperial.ac.uk/bin/linux/ubuntu lucid/ into /etc/apt/sources.list and ran apt-get update as root and got the following error message: W: GPG error: http://cran.ma.imperial.ac.uk lucid/ Re

[R] Question about object permanence/marshalling

2011-08-25 Thread Albert-Jan Roskam
Hello,   I am trying to write some code that dumps R objects to the harddisk in a binary format so they can be quickly re-used later. Goal is to save time. The objects may be quite large (e.g. classes for a GUI). I was thinking that save() and load() would be suitable for this (until now I only

Re: [R] X11 problem

2011-08-25 Thread 颜林林
In order to build R-2.13.1 from source with X11, I guess you need to install X11 develop library. Please try this: sudo apt-get install libx11-dev On Thu, 2011-08-25 at 05:49 -0700, carol white wrote: > Thanks for your reply. I used the following commands: > sudo apt-get update > sudo apt-get inst

Re: [R] X11 problem

2011-08-25 Thread carol white
gpg --keyserver keyserver.ubuntu.com --recv-key E084DAB9 gpg: requesting key E084DAB9 from hkp server keyserver.ubuntu.com gpgkeys: HTTP fetch error 7: couldn't connect to host gpg: no valid OpenPGP data found. gpg: Total number processed: 0 How to solve the problem? Cheers, Carol - Origin

Re: [R] Synchronizing R libraries on N machines?

2011-08-25 Thread Rainer M Krug
On Thu, Aug 25, 2011 at 3:25 PM, Giovanni Petris wrote: > Hello! > > I am using R on two different machines (under Ubuntu and OS X, but this > is probably irrelevant) and I would like to keep the two installations > 'synchronized', in particular in terms of installed packages. For > example, if I

Re: [R] apply function to spatial grid data frame to calculate CTI

2011-08-25 Thread Jean V Adams
It's not clear to me what problem you're having. tan(data$slope) should work. As should log(data$uparea/data$slope) I suggest that you provide a small subset of example data along with an example of any code you've tried, and the output that you'd like. Jean Tom Vanwalleghem wrote on

[R] {R} How to extract correctly from vector?

2011-08-25 Thread J. Augusiak
Dear list, I have problem that I cannot solve and would like to ask your opinion. I tried to ask a few days ago already but got no answer and all my attempts to solve it by myself since then failed. Sorry for repeated posting! Here the problem broken down a bit. My problem basically is, tha

[R] Sorting order of reorder with multiple variables

2011-08-25 Thread markm0705
I've been building a ranked dot plot for several days now and am sorting the data using the reorder command. What I don't understand is how reorder works when mutiple varibles are plotted by grouping. In the example below I'm using re-order to sort by a variable name Resv_Prop, but I'm plotting u

Re: [R] X11 problem

2011-08-25 Thread carol white
Thanks I ran your command but still get the same error message when I run ./configure configure: error: --with-x=yes (default) and X11 headers/libs are not available - Original Message - From: Linlin Yan (颜林林) To: carol white Cc: "r-h...@stat.math.ethz.ch" Sent: Thursday, August 25

Re: [R] {R} How to extract correctly from vector?

2011-08-25 Thread Sarah Goslee
I suspect you need a get() in there somewhere, but this isn't a completely reproducible example so I didn't try it. get(filename) will return the R object whose name is given by the string filename. Sarah On Thu, Aug 25, 2011 at 9:06 AM, J. Augusiak wrote: > Dear list, > > > > I have problem th

Re: [R] {R} How to extract correctly from vector?

2011-08-25 Thread jim holtman
Do you want to use 'get' in the last statement: for (i in 1:10){assign(nam[i],as.ltraj(xy=get(char[i]),date=DATE,...) I would also suggest the you use a 'list' instead of creating a lot of variables in the workspace; e.g., result <- lapply(char, function(x) as.ltraj(xy = get(x), date = DATE, ...

Re: [R] {R} How to extract correctly from vector?

2011-08-25 Thread R. Michael Weylandt
It sounds like if you know you can pass the string "A1" where you need it to go, the command get() *might* be what you are looking for. Also, you can make your code much faster/prettier because R is smart about recycling vectors of different lengths: char <- paste("A",1:10, sep="") x <- paste("X

Re: [R] X11 problem

2011-08-25 Thread Rainer M Krug
On Thu, Aug 25, 2011 at 3:51 PM, carol white wrote: > Thanks I ran your command but still get the same error message when I run > ./configure > > configure: error: --with-x=yes (default) and X11 headers/libs are not > available > If you are installing via apt-get, why do you have to run ./confi

Re: [R] X11 problem

2011-08-25 Thread Karl Brand
Hi Carol, When i ran into this problem, using the instructions as guide from here: http://cran.r-project.org/bin/linux/ubuntu/ I did the following: Went to this site: http://keyserver.ubuntu.com:11371/ searched for: 0xE084DAB9 Then clicked the hyperlinked "E084DAB9", copy, pasted and saved a

Re: [R] Question about object permanence/marshalling

2011-08-25 Thread jim holtman
The problem I think is in your unmarshal. 'load' will load the object into the local environment, not the global. You have to explicitly return it, that means you have to know the name that its was 'save'd by; > unmarshal <- function(xdr) { >     object <- strsplit(strsplit(xdr, "\\.")[[1]][[1]]

Re: [R] X11 problem

2011-08-25 Thread 颜林林
Sorry for the wrong package name. It should be: sudo apt-get install xorg-dev On Thu, 2011-08-25 at 06:51 -0700, carol white wrote: > Thanks I ran your command but still get the same error message when I run > ./configure > > configure: error: --with-x=yes (default) and X11 headers/libs are not

[R] Odp: {R} How to extract correctly from vector?

2011-08-25 Thread Petr PIKAL
> > Dear list, > > > > I have problem that I cannot solve and would like to ask your opinion. I > tried to ask a few days ago already but got no answer and all my attempts to > solve it by myself since then failed. Sorry for repeated posting! Here the > problem broken down a bit. > > > > M

Re: [R] X11 problem

2011-08-25 Thread 颜林林
Because the official source is not as new as the latest version yet. On Thu, 2011-08-25 at 15:59 +0200, Rainer M Krug wrote: > > > On Thu, Aug 25, 2011 at 3:51 PM, carol white > wrote: > Thanks I ran your command but still get the same error message > when I run ./configure >

[R] Choropleth in R

2011-08-25 Thread Reza Salimi-Khorshidi
Hi all, I would like to use R in order to generate a Heatmap/Choropleth overlaid on the UK map. Can someone please help with finding the right package/code? Any help is much appreciated. Best, Reza [[alternative HTML version deleted]] __ R-help@

Re: [R] Importing data from MS EXCEL (.xls) to R XXXX

2011-08-25 Thread Kathleen Rollet
Hi Dan, You can save your .xls as .csv, set up your directory under the file where you saved it, then widge= read.csv("WidgeOne.csv", header=T) It should work. Kathleen. > From: jorgeivanve...@gmail.com > Date: Wed, 24 Aug 2011 18:31:21 -0400 > To: dan.abne...@gmail.com > CC: r-help@r

Re: [R] Create two uniformly random variables correlated

2011-08-25 Thread Enrico Schumann
This topic has been discussed (quite recently) on R-help. http://r.789695.n4.nabble.com/Generating-uniformly-distributed-correlated-da ta-td3314905.html Enrico > -Ursprüngliche Nachricht- > Von: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] Im Auftrag von Soberon

Re: [R] how to read a group of files into one dataset?

2011-08-25 Thread Hadley Wickham
> # Method 2: Use the plyr package > > library('plyr') > bdf <- ldply(mlply(files, read.csv, header = TRUE), rbind) Or just bdf <- ldply(files, read.csv, header = TRUE) Hadley -- Assistant Professor / Dobelman Family Junior Chair Department of Statistics / Rice University http://had.co.nz/ _

[R] plotting muliple CI ellipses for lm estimates

2011-08-25 Thread natalia . vizcaino . palomar
Hi, I have read your post about plotting multiple CI ellipses and I would like too, but I don´t know how. I would like to plot an ellipse for each group and I have slope and intercepts, but I don´t know what to put in the arguments shape and radius. CAn you help me please? __

Re: [R] rpart: plot without scientific notation

2011-08-25 Thread Jean V Adams
I don't see a quick solution to this. You could contact the maintainer of the rpart.plot package, Stephen Milborrow maintainer("rpart.plot") or you could try to modify the rpart.plot() function yourself to meet your needs rpart.plot Jean Jay wrote on 08/25/2011 05:30:25 AM: > >

Re: [R] Choropleth in R

2011-08-25 Thread Jean V Adams
Check out the maps package for a choropleth map. Specifically the map() function (and its argument col=). ?map Check out the image() function for a heat map which can be overlaid. ?image Jean Reza Salimi-Khorshidi wrote on 08/25/2011 09:00:04 AM: > > Hi all, > I would like to use

Re: [R] apply function to spatial grid data frame to calculate CTI

2011-08-25 Thread Paul Hiemstra
On 08/25/2011 10:54 AM, Tom Vanwalleghem wrote: > Dear R helpers, > I'm new to R and struggling with applying functions to a spatial grid data > frame, called data for example.data has a variable named "slope" and "uparea" > Basically, I need to calculate the tan(data$slope)Next, I need to divide

Re: [R] Function rank() for data frames (or multiple vectors)?

2011-08-25 Thread David Winsemius
On Aug 25, 2011, at 7:56 AM, Sebastian Bauer wrote: Hi! On 08/24/2011 07:46 PM, David Winsemius wrote: I was looking for an elegant solution ;) In the real case I have double values and this would be quite inefficient then. Still no r-code: Then what about rank(order(...) , further-ties.

[R] Combining a greek letter and a newline on the same label

2011-08-25 Thread Sébastien Vigneau
Hi, I would like to add a two-lines label consisting of: "Delta word1\n word2" (with Delta being the actual greek letter). However, when I try "expression(paste(Delta,"word1\nword2"))", I obtain the following display: "word1\n Delta word2". Do you have any idea on how to obtain the desired displ

Re: [R] X11 problem

2011-08-25 Thread carol white
Thanks I could run configure and make but got warning after running configure: configure: WARNING: you cannot build DVI versions of the R manuals configure: WARNING: you cannot build DVI versions of all the help pages configure: WARNING: you cannot build info or HTML versions of the R manuals conf

Re: [R] Question about object permanence/marshalling

2011-08-25 Thread Albert-Jan Roskam
Hi Jim,   Thanks for your reply. It seems that save and load can only be used for datasets (as the title in ?load suggests). I'd be very glad if I'm mistaken though!  Cheers!! Albert-Jan ~~ All right, but apart from the sanitat

Re: [R] rpart: plot without scientific notation

2011-08-25 Thread David Winsemius
On Aug 25, 2011, at 10:38 AM, Jean V Adams wrote: I don't see a quick solution to this. You could contact the maintainer of the rpart.plot package, Stephen Milborrow maintainer("rpart.plot") or you could try to modify the rpart.plot() function yourself to meet your needs rpart.pl

[R] R hangs after htmlTreeParse

2011-08-25 Thread Simon Kiss
Dear colleagues, I'm trying to parse the html content from this webpage: http://timesofindia.indiatimes.com/searchresult.cms?sortorder=score&searchtype=2&maxrow=10&startdate=2001-01-01&enddate=2011-08-25&article=2&pagenumber=1&isphrase=no&query=IIM&searchfield=§ion=&kdaterange=30&date1mm=01&date1dd

Re: [R] read.table truncated data?

2011-08-25 Thread zhenjiang xu
Thanks for your replies. I looked at those lines and didn't spot anything unusual. > tail(a) test_id gene_id gene locus sample_1 sample_2 status 21418 tY(GUA)J1 - SUP7 chr10:354243-354332 air1rrp6 air2rrp6 OK 21419 tY(GUA)J2 - SUP4 chr10:542955-543044 air1rrp6

Re: [R] X11 problem

2011-08-25 Thread Dirk Eddelbuettel
On 25 August 2011 at 22:24, Linlin Yan (颜林林) wrote: | Because the official source is not as new as the latest version yet. That is not true --- you get the most recent R version prebuilt for Ubuntu (and other Linux distros) via the CRAN repositories. See the R FAQ about where to get to current U

Re: [R] Question about object permanence/marshalling

2011-08-25 Thread jim holtman
You probably need to change 'marshal' also. I see that you are passing in the character string and this is what you are saving. what you probably is something like this where you pass in the object and then get its name: >> marshal <- function(object.value) { object <- deparse(substitut

Re: [R] read.table truncated data?

2011-08-25 Thread jim holtman
But did you try the following: x <- read.table(, comment.char = '', quote = '') Most cases is that there is a missing quote somewhere in your data. use a text editor and search for single and double quotes. On Thu, Aug 25, 2011 at 11:49 AM, zhenjiang xu wrote: > Thanks for your replies. I l

Re: [R] Append a value to a vector

2011-08-25 Thread Claudio Zanettini
Oh thanks, yes now I understand. I did not realized that I did not save it. :) 2011/8/24 Jean V Adams > > Claudio Zanettini wrote on 08/24/2011 > 04:33:50 PM: > > > > > Thank you, this work fine, > > and is not contorted like mine:) > > In this case lastV=LastI but depending on the data that I

Re: [R] Combining a greek letter and a newline on the same label

2011-08-25 Thread David Winsemius
On Aug 25, 2011, at 11:23 AM, Sébastien Vigneau wrote: Hi, I would like to add a two-lines label consisting of: "Delta word1\n word2" (with Delta being the actual greek letter). But you are not saying how word1 and word2 are supposed to be positioned relative to each other and to Delta.

Re: [R] Change color in forest.rma (metafor)

2011-08-25 Thread Viechtbauer Wolfgang (STAT)
The color of the squares is also currently hard coded. The thing is, there are so many different elements to a forest plot (squares, lines, polygons, text, axes, axis labels, etc.), if I would add arguments to set the color of each element, things would really get out of hand (as far as I am c

Re: [R] Split data frame by date (POSIXlt)

2011-08-25 Thread MacQueen, Don
I suspect, but have not tested, that your src$date element has class "POSIXlt", which is internally a list, in which case splitting by it might not work properly, and might be the cause of your out of bounds error message. One of these might do the job: src$date <- format(strptime(tickdata$da

Re: [R] read.table truncated data?

2011-08-25 Thread zhenjiang xu
Thanks, Jim. quote='' works. And then I found a single quote in each of these lines: 3262 10403 17544 24685 31826 38967 None of them near the position the table got truncated. Why is it? And read.table is a great function. Is it possible for it to give a warning message when the data gets truncat

Re: [R] X11 problem

2011-08-25 Thread 颜林林
Thank you for the information! On Thu, 2011-08-25 at 10:53 -0500, Dirk Eddelbuettel wrote: > On 25 August 2011 at 22:24, Linlin Yan (颜林林) wrote: > | Because the official source is not as new as the latest version yet. > > That is not true --- you get the most recent R version prebuilt for Ubuntu

[R] Construct a File Path: File Path Unknown

2011-08-25 Thread Tyler Rinker
I am not a programmer and am self-taught so I may lack the language to ask this appropriately (perhaps why an rseek search was unfruitful). Let's say I saved a file to my desktop called foo.pdf. Then I want R to return the file path of foo.pdf (pretend I don't know the location(path) of

[R] Optim function with multivariate inputs

2011-08-25 Thread Noah Silverman
Hi, I have function that I want to optimize. Am playing with the optim() function in R Two issues: 1) I can't seem to get it to work with a function that takes multiple inputs. Dummy Example: myFunc <- function(A,B,D,D){ # Do stuff return E } > myFunc(1,2,3,4) [1] 12 # wor

[R] Question about FlexMix package

2011-08-25 Thread Andra Isan
Hello All, I am trying to use FlexMix package for my two linear regression models: http://finzi.psych.upenn.edu/R/library/flexmix/doc/regression-examples.pdf y1 = w1 * x1 + w2 *x2 + w3 * x3 +beta1 y2 = w4 * x4 + w4 *x4 + w4 * x4 +beta2 and I would like to combine these two regression models wi

[R] Segment out of the Graph

2011-08-25 Thread Claudio Zanettini
Hello everyone, I have a graph and a segment parallel to the x axis at y=-10, x=0, and bars on it. Now the question is, Is there a way to leave the segment there but let the graph axis start from the origin? In this way the segment will be out of the graph Thanks [[alternative HTML ver

Re: [R] Segment out of the Graph

2011-08-25 Thread R. Michael Weylandt
Look at ylim, as an optional argument to plot. Michael On Thu, Aug 25, 2011 at 1:07 PM, Claudio Zanettini < claudio.zanett...@gmail.com> wrote: > Hello everyone, > I have a graph and a segment parallel to the x axis at y=-10, x=0, and > bars on it. > Now the question is, > Is there a way to le

Re: [R] Segment out of the Graph

2011-08-25 Thread Claudio Zanettini
Yes I tried but if I set the lim to 0 then it will not displayed the line that is at -10, right? 2011/8/25 R. Michael Weylandt > Look at ylim, as an optional argument to plot. > > Michael > > On Thu, Aug 25, 2011 at 1:07 PM, Claudio Zanettini < > claudio.zanett...@gmail.com> wrote: > >> Hello ev

Re: [R] Segment out of the Graph

2011-08-25 Thread R. Michael Weylandt
"lim" is not the argument: "ylim" is. You put in a vector of length 2 comprising the min and max y you wish: consider this: x = -5:5; y = x^2; z = rep(-5,11); layout(1:2) plot(x,y,type="b"); lines(x,z,col=2) plot(x,y,ylim = c(-8,max(y)+3),type="b"); lines(x,z,col=2) For your work, you'd need yli

Re: [R] Segment out of the Graph

2011-08-25 Thread Claudio Zanettini
> I tried setting ylim=(0, 1000) > but the segment that I have draw is at y=-10 > so if I set the y origin to 0 I don t have the segment, > if a l leave it at -10 I have the segment but the axis start from -10. > > I would like to have both. > So to have a graph with ylim=( 0, 1000) > and under it

Re: [R] Segment out of the Graph

2011-08-25 Thread R. Michael Weylandt
I'm not really sure what you are asking for, but I might venture to suggest a few ideas: 1) Do you mean to say you want two plots: one with ylim around -10 and the other running from 0 to 1000? 2) Are you looking for a "break" in the y-axis? If so, plotrix may be able to help, but it's hardly wor

Re: [R] Segment out of the Graph

2011-08-25 Thread David Winsemius
On Aug 25, 2011, at 1:30 PM, Claudio Zanettini wrote: I tried setting ylim=(0, 1000) but the segment that I have draw is at y=-10 so if I set the y origin to 0 I don t have the segment, if a l leave it at -10 I have the segment but the axis start from -10. I would like to have both. So to h

Re: [R] Easier ways to create .Rd files?

2011-08-25 Thread Greg Snow
For the initial .Rd file look at the 'prompt' function, it will create a basic .Rd file based on the current function definition with all the structure in place, then you go in and edit it to have more meaningful information. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain

Re: [R] read.table truncated data?

2011-08-25 Thread jim holtman
When you have an unbalanced quote, it may be hard to determine exactly where it is. It is probably up to the user to determine with there is truncation. In some cases you might have data that goes over several lines that are within quotes and is legal. You might also read up on the 'fill' and 'f

[R] Fwd: Segment out of the Graph

2011-08-25 Thread Claudio Zanettini
-- Forwarded message -- From: Claudio Zanettini Date: 2011/8/25 Subject: Re: [R] Segment out of the Graph To: David Winsemius Thanks David and Michael, In attachment there is one of the graph. the line below the graph is not related to the y label. I draw it at y= -10 just beco

Re: [R] Choropleth in R

2011-08-25 Thread Reza Salimi-Khorshidi
Thanks Jean, but this package appears to be perfect for the US + a couple of other country that unfortunately UK is not one of them. What I want is a solution that takes a list of coordinates/postcodes/cities and a list of values and gives me a colourful UK map. Any thoughts? Cheers On Thu, Aug 25

[R] Adding a normal density curve over the empirical curve

2011-08-25 Thread Newbie
Hi I have created the following plot over the empirical returns.. What I now want to do is to overlay a curve/line with the normal density as a comparison of the two. Does anyone know how to do this? (NB the last two lines are the problem, and are wrong, I know). Thank you in advance! Rikke http

[R] Using write.table i have a table with two columns i would like to save it as an excel file

2011-08-25 Thread kevin123
Using write.table i would like to save data as an excel file to a folder. I am not too sure how to write the file path or what to name the file. I would appreciate any feedback. > write.table(x, file = "", append = FALSE, quote = TRUE, sep = " ", + eol = "\n", na = "NA", dec = ".",

[R] within-groups variance and between-groups variance

2011-08-25 Thread Coghlan, Avril
Hello, I have been looking for functions for calculating the within-groups variance and between-groups variance, for the case where you have several numerical variables describing samples from a number of groups. I didn't find such functions in R, so wrote my own versions myself (see below). I ca

[R] Selections in lists

2011-08-25 Thread mdvaan
Hi, I have produced a list g and I would like to reduce the amount of information contained in each object in g. For each matrix I would like to keep the values where the column name equals g[year][[1]][[x]] and the row names equals g[year][[1]][[-x]]. So in g$`1999`$`8029`, year = 1999 and x = 8

[R] how to match a dataframe and a list

2011-08-25 Thread Serena105
Hi! I am a beginner and I hope you will help me! I have one data frame with numbers as columns names that correspond to codes. In a list I have the codes and the corresponding state. Now I need to match the columns (codes) of the data.frame with the codes in the list, since not all the codes in

Re: [R] Construct a File Path: File Path Unknown

2011-08-25 Thread Jean V Adams
Try the dir() function. ?dir # for example dir("c:/", pattern="foo.pdf", full.names=T, ignore.case=T, recursive=T) Jean Tyler Rinker wrote on 08/25/2011 11:54:28 AM: > > I am not a programmer and am self-taught so I may lack the > language to ask this appropriately (perhaps why an rseek searc

[R] Bivariate normal regression in R

2011-08-25 Thread DHIMAN BHADRA
Hello everyone, I need to fit a bivariate normal regression model to a dataset where the same covariate (say, X) influences two separate but correlated responses (say, Y1 and Y2). So, the bivariate model would look like : Y1 = a1 + b1*X + e1 Y2 = a2 + b2*X + e2 where e1 and e2 are error terms wh

  1   2   >