Re: [R] Need to append vector to all levels of nested list WITHOUT a loop

2012-06-18 Thread Rui Barradas
Hello, Try (b <- list(list(m=4, v=5), list(m=6, v=7))) (m <- matrix(1:6, 2)) b2 <- lapply(seq_along(b), function(i){ b[[i]]$line <- m[i, ]; b[[i]]}) b2 Also, next time, post a data example, using dput(). Hope this helps, Rui Barradas Em 18-06-2012 14:54, Johannes Reichl escreveu:

Re: [R] Inconsistency using seq

2012-06-18 Thread Petr Savicky
On Mon, Jun 18, 2012 at 12:58:13PM -0700, hamoreno wrote: > Hi all, > > Is there any problem of precision when using seq?. For example: > > x<- seq(0,4,0.1) > x[4]=0.3 > > BUT: > > x[4]-0.3=5.551115e-17 > > It means when I use this condition within an if clause, it does not find > values with

Re: [R] (1-1e-100)==1 true?

2012-06-18 Thread Petr Savicky
On Mon, Jun 18, 2012 at 02:26:41PM -0700, whf1984911 wrote: > Hi, > > This problems has bothered me for the lase couple of hours. > > > 1e-100==0 > [1] FALSE > > (1-1e-100)==1 > [1] TRUE > > How can I tell R that 1-1e-100 does not equal to 1, actually, I found out > that > > (1-1e-16)==1 >

Re: [R] how to get help for a package. Thanks!

2012-06-18 Thread Jorge I Velez
James, Try library(osmar) library(help = "osmar") HTH, Jorge.- On Mon, Jun 18, 2012 at 9:39 PM, james pruett <> wrote: > Hi, > > I installed "osmar". > Q: How do I learn what methods it contains. > > I tried > > Help, things like that > > Thanks > cellurl > >[[alternative HTML vers

Re: [R] how to get help for a package. Thanks!

2012-06-18 Thread Pascal Oettli
Hello, Search by yourself. http://cran.r-project.org/web/packages/osmar/index.html Regards. Le 12/06/19 10:39, james pruett a écrit : Hi, I installed "osmar". Q: How do I learn what methods it contains. I tried Help, things like that Thanks cellurl [[alternative HTML version

Re: [R] how to get help for a package. Thanks!

2012-06-18 Thread R. Michael Weylandt
The help function in R is just that -- a function -- so you have to invoke it, e.g., as help(mean) rather than just typing Help. You might also try help.start() which opens a fuller help system in your browser. To see the available functionality of the package, the easiest way is to check its CRAN

[R] Partial correlation analysis for censored data in R?

2012-06-18 Thread Marco Chiaberge
Hi does anyone know if there is any R function that preforms partial correlations analysis for censored data? I found a fortran program on the PSU website at http://www2.astro.psu.edu/statcodes/cens_tau.f but I was wondering if there is anything like that in any R package. Thanks m. _

[R] how to get help for a package. Thanks!

2012-06-18 Thread james pruett
Hi, I installed "osmar". Q: How do I learn what methods it contains. I tried > Help, things like that Thanks cellurl [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PL

[R] help with xy.coords(x,y)

2012-06-18 Thread adamshun
i am working on the project to analyze hedge fund performance, i would appreciate that if you guys could spare some time helping me out with the R code. Thanks. The senario is: i applied BOXPLOT() to plot the performance of all hedge funds with 7 strategies. And right now in this boxplot I nee

Re: [R] Changing strip text

2012-06-18 Thread Duncan Mackay
Hi Julie Try strip.custom(factor.levels = paste("Treat", c("0","2","4","8","16"), sep = " ") i.e xyplot(log(pfuml)~Transfer|Treat,groups=Pop,ylab="Log10 Pfu/ml",as.table=T, type="a", auto.key=list(space="right",points=F,lines=T) ), strip=strip

[R] how to use by function

2012-06-18 Thread cowboy
hi all, Assume I have data like data<-rbind(c(1,2),c(1,3),c(2,1),c(3,2),c(3,4)) I want to get some matrix like 1,2,3 2,NA,NA 3,2,4 I'm using by mat<-matrix(NA,3,3) by(data,data[,1],mat[data[,1],]<-c(data[,2])) but it doesn't work. Any ideas? thanks, cowboy _

Re: [R] Help with readBin

2012-06-18 Thread kapo coulibaly
Yes, the numbers look right. Thank you very much. I think the "junk" header was throwing my code off. For some reason I was missing it. Now everything is clear. Once again thanks a bunch. On Mon, Jun 18, 2012 at 9:43 PM, William Dunlap wrote: > You didn't give much of a description of what sort

Re: [R] Error while installing R on RHEL 6

2012-06-18 Thread John
On Mon, 18 Jun 2012 02:26:32 -0700 (PDT) Manish Gupta wrote: > Hi, > > I am trying to install R from source and using configure command but > gettting one error message. > > *configure: error: --with-readline=yes (default) and headers/libs are > not available* > > I am installing on server wh

Re: [R] Help with readBin

2012-06-18 Thread William Dunlap
I suspect that "junk" (the four bytes after the 2nd record marker), read as a 4-byte integer instead of four 1-byte integers, is the number of bytes of data following it. Its value in your example is 18920 = 8 * 43 * 55, where 43 and 55 are two integers in the header, probably the dimensions o

Re: [R] Help with readBin

2012-06-18 Thread William Dunlap
You didn't give much of a description of what sort of numbers you expected in the header so this is pretty much a guess. However, by reading the tail of the file with offsets 0 through 7 bytes we get numbers in the 30-40 range for an offset of 4 bytes. I called that field "junk" below and placed

Re: [R] Help with readBin

2012-06-18 Thread kapo coulibaly
No. But here it is: c(52L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, -16L, 63L, 0L, 0L, 0L, 0L, 0L, 0L, -16L, 63L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 72L, 69L, 65L, 68L, 43L, 0L, 0L, 0L, 55L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 52L, 0L, 0L, 0L, -24L, 73L,

Re: [R] Help with readBin

2012-06-18 Thread kapo coulibaly
I still haven't found a working solution. Is it allowed to attach a file so that somebodyelse can reproduce the problem? On Thu, May 3, 2012 at 5:00 PM, William Dunlap wrote: > You can do the following to allow others to recreate your problem. > > yourFileBytes <- readBin("yourFile", what="inte

Re: [R] Efficient distance calculation on big matrix

2012-06-18 Thread Stefan Evert
> It so happens I have been looking at very similar changes, as well as > adding multi-threading support for dist(); these should make it into > R-devel later this summer. That's good to hear! I was thinking about giving OpenMP a try in my package, but am not sure whether it's worth the overhea

Re: [R] Trying to speed up an if/else statement in simulations

2012-06-18 Thread Rui Barradas
Hello, Try creating an index vector. In this case, I've called it 'one'. Then, 'n.one' is the number of ones in flag1. (Since data is an R function, I've renamed your example dat). n <- 10 # number of observations # First make up some data set.seed(123) dat <- data.frame(flag1=rbinom(n, 1,

Re: [R] Inconsistency using seq

2012-06-18 Thread Henrik Bengtsson
FYI, > isZero <- function(x, neps=1, eps=.Machine$double.eps, ...) { + (abs(x) < neps * eps) + } > x <- seq(from=0, to=4, by=0.1) > isZero(x[4]-0.3) [1] TRUE > isZero(x[4]-0.3, neps=1) [1] TRUE > isZero(x[4]-0.3, neps=0.1) [1] FALSE You could also have called isZero() isValueSmallEnoughForWhatI

Re: [R] (1-1e-100)==1 true?

2012-06-18 Thread Ted Harding
On 18-Jun-2012 21:26:41 whf1984911 wrote: > Hi, > > This problems has bothered me for the lase couple of hours. > >> 1e-100==0 > [1] FALSE >> (1-1e-100)==1 > [1] TRUE > > How can I tell R that 1-1e-100 does not equal to 1, actually, > I found out that > > (1-1e-16)==1 > [1] FALSE >> (1-1e-1

Re: [R] Listing all binary trees of an ordinal set

2012-06-18 Thread David Winsemius
On Jun 18, 2012, at 5:25 PM, elmo wrote: Hi all, I'm new to R, have been reading books and trying to get started coding too. The first thing of substance I've been trying to do is to create a function to return a list of all binary trees of a list of ordinals. So, for example, an input o

Re: [R] S4 into a data frame

2012-06-18 Thread Martin Morgan
On 06/18/2012 01:49 PM, spf385 wrote: Greetings, I'm still very new to R; however I have been tasked with converting an S4 object into a data.frame. We are just trying to learn the basics of S4. Could anyone offer a 'simple' example of coercing an S4 object (such as might be stored in a package

Re: [R] Changing many csv files using apply?

2012-06-18 Thread R. Michael Weylandt
sapply() won't be notably faster. What you might try (no guarantee) is wrapping things up as if you were going to use sapply() but instead use mclapply() from the parallel package -- that will parallelize the results and should be faster by roughly as many cores as you use. The "no guarantee" disc

Re: [R] (1-1e-100)==1 true?

2012-06-18 Thread David Smith
This is standard behaviour for a floating-point computational system like R. You might like to take a look at this blog post for a backgrounder on floating-point arithmetic in R. http://blog.revolutionanalytics.com/2009/03/when-is-a-zero-not-a-zero.html # David Smith On Mon, Jun 18, 2012 at 2:26

Re: [R] Error w/ Start up Script

2012-06-18 Thread Lopez, Dan
Hi Michael, I wasn't getting how to actually use it from ?suppressPackageStartupMessages example online. Thank you very much. Dan -Original Message- From: R. Michael Weylandt [mailto:michael.weyla...@gmail.com] Sent: Monday, June 18, 2012 2:53 PM To: Lopez, Dan Cc: R help (r-help@r-p

Re: [R] Changing strip text

2012-06-18 Thread Peter Ehlers
On 2012-06-18 10:02, Julie Truman wrote: I'm trying to create a 5-panel xyplot with each strip labelled by treatment. The code is as follows: xyplot(log(pfuml)~Transfer|Treat,groups=Pop,ylab="Log10 Pfu/ml",as.table=T, type="a", auto.key= list(spac

[R] Changing many csv files using apply?

2012-06-18 Thread Chang, Emily@OEHHA
Dear all, I have many csv files whose contents I want to change a bit en masse. So far, I've written code that can change them in a for loop, like so: # Subset of files in the folder I want to change subset = "somestring" # Retrieve list of files to change filelist=list.files() filelist = fileli

Re: [R] Error w/ Start up Script

2012-06-18 Thread R. Michael Weylandt
You are perhaps not using suppressPackageStartupMessages() correctly if you're getting those warnings: my .Rprofile has these lines: suppressPackageStartupMessages(library("quantmod")) suppressPackageStartupMessages(library("ggplot2")) suppressPackageStartupMessages(library("lattice")) but produc

Re: [R] Trying to speed up an if/else statement in simulations

2012-06-18 Thread David Winsemius
On Jun 18, 2012, at 1:29 PM, nqf wrote: Dear R-help, I am trying to write a function to simulate datasets of size n which contain two time-to-event outcome variables with associated 'Event'/'Censored' indicator variables (flag1 and flag2 respectively). One of these indicator variables nee

Re: [R] Trying to speed up an if/else statement in simulations

2012-06-18 Thread R. Michael Weylandt
I might try something like: data[data$flag1 == 1, "flag2"] <- runif(sum(data$flag1 == 1)) < 0.95 and similarly for the other case. Hope this helps, Michael On Mon, Jun 18, 2012 at 12:29 PM, nqf wrote: > Dear R-help, > > I am trying to write a function to simulate datasets of size n which conta

Re: [R] Inconsistency using seq

2012-06-18 Thread Jeff Newmiller
actually, you should build your sequences with integers and scale those to get floating point sequences. --- Jeff NewmillerThe . . Go Live... DCN:Basics: ##.#. ##.#. Li

Re: [R] Trying to speed up an if/else statement in simulations

2012-06-18 Thread Peter Alspach
Tena koe Natalie Does this do what you want? set.seed(123) natalie <- matrix(NA, nrow=10^5, ncol=2, dimnames=list(NULL, paste0('flag', 1:2))) natalie[,'flag1'] <- rbinom(nrow(natalie), 1, 0.5) natalie[natalie[,'flag1']==1, 'flag2'] <- rbinom(sum(natalie[,'flag1']), 1, 0.95) natalie[natalie[,'fl

Re: [R] Inconsistency using seq

2012-06-18 Thread Bert Gunter
Nope. round() is not failsafe either. Basically, there is theoretically no way to guarantee immunity to floating point error in conditional comparisons. You need to either switch to integers -- which **are** exactly represented -- or figure out another way. Of course, practically speaking, suffici

Re: [R] How can I "declutter"/make a biplot less messy in R

2012-06-18 Thread Sarah Goslee
Hi Derek, It's good form to provide an example with code so that we can see what functions you're using and what you've tried. I'm going to assume that you're using princomp(). Here's one way to reduce the clutter: pc.cr <- princomp(USArrests, cor = TRUE) biplot(pc.cr, xlabs=rep("x", nrow(USArre

[R] Listing all binary trees of an ordinal set

2012-06-18 Thread elmo
Hi all, I'm new to R, have been reading books and trying to get started coding too. The first thing of substance I've been trying to do is to create a function to return a list of all binary trees of a list of ordinals. So, for example, an input of list(3,1,2,4) would return: list(list(1, list(

Re: [R] How can I "declutter"/make a biplot less messy in R

2012-06-18 Thread Bryan Hanson
Don't do the biplot. Darn hard to make sense of anyway. Plot the scores and the loadings separately. You can see how to do that in this thread: http://r.789695.n4.nabble.com/How-to-plot-PCA-output-td4614732.html Good Luck. Bryan *** Bryan Hanson Professor of Chemistry & Biochemistry

[R] (1-1e-100)==1 true?

2012-06-18 Thread whf1984911
Hi, This problems has bothered me for the lase couple of hours. > 1e-100==0 [1] FALSE > (1-1e-100)==1 [1] TRUE How can I tell R that 1-1e-100 does not equal to 1, actually, I found out that > (1-1e-16)==1 [1] FALSE > (1-1e-17)==1 [1] TRUE The reason I care about this is that I was try to u

[R] S4 into a data frame

2012-06-18 Thread spf385
Greetings, I'm still very new to R; however I have been tasked with converting an S4 object into a data.frame. We are just trying to learn the basics of S4. Could anyone offer a 'simple' example of coercing an S4 object (such as might be stored in a package attribute) into a data.frame? Ultimately

[R] Biplot - how can I make it more readable?

2012-06-18 Thread Derek . Schanze
Hello, I am doing a principle component analysis on a dataset with a lot of different variables and have constructed a biplot of the data. Unfortunately, as can be seen on the attached image, the biplot is very messy, cluttered, and hard to read. I have performed a few modifications including

Re: [R] Trying to speed up an if/else statement in simulations

2012-06-18 Thread William Dunlap
inline below Bill Dunlap Spotfire, TIBCO Software wdunlap tibco.com > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On > Behalf > Of nqf > Sent: Monday, June 18, 2012 10:30 AM > To: r-help@r-project.org > Subject: [R] Trying to speed up an

Re: [R] Using Sub

2012-06-18 Thread arun
Hi Dan, I am using linux.  So, not sure whether this will work or you tried it before.  After you paste the path at (1:), then press enter two times. One more point.  In your path, there were spaces within slashes.  So, if you do the scan, it will return as 5 items, which you have to join again

[R] Installing xlsx package on Mac OS X

2012-06-18 Thread Victoria Xiao
Hi all, I'm using a Mac port of R on Mac 10.6. I need to install xlsx package but it does not install. Do any of you know a way to install it? Thanks very much. Sincerely, Victoria Xiao [[alternative HTML version deleted]] __ R-help@r-project.

Re: [R] Using Sub

2012-06-18 Thread arun
Hi, I guess you need something like this: > str1<-"G:\\Compensation Audits_Reports\\Step Audit\\Steps off Step.accdb" > gsub('','/',str1) [1] "G:/Compensation Audits_Reports/Step Audit/Steps off Step.accdb" A.K. - Original Message - From: "Lopez, Dan" To: Duncan Murdoch Cc: R

[R] How can I "declutter"/make a biplot less messy in R

2012-06-18 Thread dss
I am doing a principle component analysis on a dataset with a lot of different variables and have constructed a biplot of the data. Unfortunately, as can be seen on the attached image, the biplot is very messy, cluttered, and hard to read. I have performed a few modifications including outlier re

Re: [R] Inconsistency using seq

2012-06-18 Thread hamoreno
Thanks everyone... Seems that I will have to use round before seq to make sure everything has the correct precision. -- View this message in context: http://r.789695.n4.nabble.com/Inconsistency-using-seq-tp4633739p4633750.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] Error w/ Start up Script

2012-06-18 Thread Lopez, Dan
Ouch Bert. Not patience for the newbies, I see. I know what ?anything means which is why I didn't ask about ?Startup. But being that I just recently discovered what Rprofile.site is I thought it might do something different within the start-up file and tried it both ways. Besides I think what I

Re: [R] Inconsistency using seq

2012-06-18 Thread R. Michael Weylandt
No, this is rather the nature of floating point calculations. You may perhaps be looking for ?all.equal or R FAQ 7.31 (I think that's the one) which is google-able. It's a complicated subject, but those should get you started. Best, Michael On Jun 18, 2012, at 2:58 PM, hamoreno wrote: > Hi

Re: [R] Inconsistency using seq

2012-06-18 Thread Sarah Goslee
Please read R FAQ 7.31. The problem is with your computer (and everyone else's), and has nothing to do with seq(). On Mon, Jun 18, 2012 at 3:58 PM, hamoreno wrote: > Hi all, > > Is there any problem of precision when using seq?. For example: > > x<- seq(0,4,0.1) > x[4]=0.3 > > BUT: > > x[4]-0.3=5

Re: [R] Inconsistency using seq

2012-06-18 Thread Milan Bouchet-Valat
Le lundi 18 juin 2012 à 12:58 -0700, hamoreno a écrit : > Hi all, > > Is there any problem of precision when using seq?. For example: > > x<- seq(0,4,0.1) > x[4]=0.3 > > BUT: > > x[4]-0.3=5.551115e-17 > > It means when I use this condition within an if clause, it does not find > values with 0.

Re: [R] Error w/ Start up Script

2012-06-18 Thread Bert Gunter
I gave you the references you need. It's up to you to read and understand how to use them. You apparently need to work harder to do so. BTW ?anything is an alternative way to type help("anything") at the command line for R's help pages. -- Bert On Mon, Jun 18, 2012 at 1:34 PM, Lopez, Dan wro

[R] Trying to speed up an if/else statement in simulations

2012-06-18 Thread nqf
Dear R-help, I am trying to write a function to simulate datasets of size n which contain two time-to-event outcome variables with associated 'Event'/'Censored' indicator variables (flag1 and flag2 respectively). One of these indicator variables needs to be dependent on the other, so I am creating

[R] Inconsistency using seq

2012-06-18 Thread hamoreno
Hi all, Is there any problem of precision when using seq?. For example: x<- seq(0,4,0.1) x[4]=0.3 BUT: x[4]-0.3=5.551115e-17 It means when I use this condition within an if clause, it does not find values with 0.3 for x[4] as it is not precisely 0.3. Is there any bug in seq() ? -- View this

Re: [R] noob requesting help

2012-06-18 Thread capital_P
Rui Barradas wrote > > sorry about the misleading tip. > This should do it. > > do.call(rbind, dep) > No problem at all. That did the trick. Thank you both so much for your help. I am forever in your debt. -- View this message in context: http://r.789695.n4.nabble.com/noob-requesting-help-

[R] Changing strip text

2012-06-18 Thread Julie Truman
I'm trying to create a 5-panel xyplot with each strip labelled by treatment. The code is as follows: xyplot(log(pfuml)~Transfer|Treat,groups=Pop,ylab="Log10 Pfu/ml",as.table=T, type="a", auto.key= list(space="right",points=F,lines=T)) I've tried:

[R] multiplicative error model

2012-06-18 Thread saraberta
HI, i'm trying to estimate a multiplicative error model in R, i've seen the dynamo package but it seems that i'm not able to use it. does someone know how can i do it? thanks -- View this message in context: http://r.789695.n4.nabble.com/multiplicative-error-model-tp4633730.html Sent from the

Re: [R] Loop Help

2012-06-18 Thread @ngel
Thank you Jeff, that sounds like the solution to my problem but how would I do that code-wise? could you please help? -- View this message in context: http://r.789695.n4.nabble.com/Loop-Help-tp4633558p4633701.html Sent from the R help mailing list archive at Nabble.com. _

[R] How to average time series data within certain time periods

2012-06-18 Thread Dominik van Pinxteren
Hi, I am new to R and wondering how I can subset and average time series data based on given irregular time periods: Say I have a time series of a measured variable "x" (hourly data) in "values.csv" with columns "date" and "x" which I would like to average within certain irregular time pe

Re: [R] Error w/ Start up Script

2012-06-18 Thread Lopez, Dan
Hi Bert, I tried entering: "?suppressPackageStartupMessages" and then "suppressPackageStartupMessages" into my Rprofile.site file and the messages still came up. I did however go to the Misc menu and select "List search path" and got the below results. Not exactly sure what it means but I do n

Re: [R] Error w/ Start up Script

2012-06-18 Thread Bert Gunter
?suppressPackageStartupMessages Put this into your appropriate startup file. ?Startup ## for details on starting R and what files are executed at startup. -- Bert On Mon, Jun 18, 2012 at 12:46 PM, Lopez, Dan wrote: > David, > > Sorry I should have been more specific. They are not as much erro

Re: [R] Interactive Graphics

2012-06-18 Thread MacQueen, Don
As I recall, there is a package that can be used to create interactive svg files. I don't remember its name, but I seem to recall it was fairly easy to find on the CRAN packages page. -Don -- Don MacQueen Lawrence Livermore National Laboratory 7000 East Ave., L-627 Livermore, CA 94550 925-423-1

Re: [R] Error w/ Start up Script

2012-06-18 Thread Lopez, Dan
David, Sorry I should have been more specific. They are not as much errors as they are extra script that is running at startup after "Type q() to quit R". It wasn't there before. And am thinking eventually it could become a problem. Thanks. Dan Here it is again...what it looks like when I firs

Re: [R] Error w/ Start up Script

2012-06-18 Thread David Winsemius
On Jun 18, 2012, at 11:25 AM, Lopez, Dan wrote: I noticed about a couple weeks ago that my R start up script starting showing the below errors, following the line "Type 'q()' to quit R" I didn't see any errors reported. I didn't see any warnings, either. -- David. Not sure what I did d

Re: [R] How do anova() and Anova(type="III") handle incomplete designs?

2012-06-18 Thread John Fox
Dear Justin, On Mon, 18 Jun 2012 11:24:33 -0400 Justin Montemarano wrote: > Thanks for your response, John. That was helpful. > > I was using Type III from Anova() as a comparison to some results I had > obtained JMP, which I've lost access to and have moved on to R, and I was > confused by th

Re: [R] multiplicative error model

2012-06-18 Thread Özgür Asar
What do you mean by using the package? Aren't the examples in the package manual at http://cran.r-project.org/web/packages/dynamo/dynamo.pdf enough for you? Ozgur -- View this message in context: http://r.789695.n4.nabble.com/multiplicative-error-model-tp4633730p4633732.html Sent from the R

Re: [R] Using Sub

2012-06-18 Thread Lopez, Dan
Jeff - Thanks. I think I'll just stick to using file.choose(). Dan -Original Message- From: Jeff Newmiller [mailto:jdnew...@dcn.davis.ca.us] Sent: Monday, June 18, 2012 10:29 AM To: Lopez, Dan; arun Cc: R help Subject: Re: [R] Using Sub Please study how escapes work. The representation

Re: [R] Process XML files

2012-06-18 Thread Santosh
Please ignore my previous request for clarification. I was able to resolve my problem using "getNodeSet". Thanks, Santosh On Thu, Jun 14, 2012 at 10:54 AM, Santosh wrote: > Dear Rxperts, > I am back to favoRite! I would need your favoR, please! > > Is there a way to use read "format" arguments

Re: [R] Using Sub

2012-06-18 Thread William Dunlap
The parser (the part of R that interprets text as R code) interprets things like "\t" as the tab character and there is no way to avoid that. However, when the functions readLines(), scan(), and readline() handle input they do not assign any special meanings to backslashes. Hence, if you want to

Re: [R] Using Sub

2012-06-18 Thread Jeff Newmiller
Please study how escapes work. The representation you see in the source code is not the same as what is actually stored in memory. As to your problems with the clipboard, I recommend writing scripts that do not interact with the clipboard directly. Violating this rule puts an error-prone human

Re: [R] Package of EM and MI for IRT in R

2012-06-18 Thread Bert Gunter
SEARCH! Go to rseek.org Type "missing data imputation" . You will bring up several alternatives. -- Bert On Mon, Jun 18, 2012 at 8:04 AM, R. Michael Weylandt < michael.weyla...@gmail.com> wrote: > Possibly Amelia: http://cran.r-project.org/web/packages/Amelia/index.html > > Best, > Michael >

Re: [R] Loop Help

2012-06-18 Thread Özgür Asar
@angel, You can try append option in write.table for e.g. for(i in 1:20) write.table(i,"out.txt",append=T,sep="\t",row.names=F,col.names=F) This will help you save your all results by writing each results on subsequent rows. Hope this helps Ozgur >Hello Jim >thank you so much for your respons

Re: [R] Cholesky decomposition error

2012-06-18 Thread Kjetil Halvorsen
see inline! On Mon, Jun 18, 2012 at 12:38 AM, wrote: > Thanks Kjetil for your detailed code for log-Cholesky but my situation is > like optimization of the variables inside the matrix. But you must of course call my function makemat inside the call to optim() or whatever function you use for o

Re: [R] SAP ABAP Consultant available

2012-06-18 Thread Q
Please DO NOT put ads for your services on this site. This is a violation of the R forum terms. I will be forwarding this issue to the appropriate moderator to take further action. On Tue, Jun 12, 2012 at 1:37 PM, Kanna wrote: > > [1]Click here to unsubscribe if you no longer wish to receive

Re: [R] Efficient distance calculation on big matrix

2012-06-18 Thread luke-tierney
It so happens I have been looking at very similar changes, as well as adding multi-threading support for dist(); these should make it into R-devel later this summer. Best, luke On Sun, 17 Jun 2012, Stefan Evert wrote: I'm working on analyzing a large data set, lets asume that dim(Data)=c(10

Re: [R] Using Sub

2012-06-18 Thread Lopez, Dan
Hi A.K., That works but when I copy and paste from Window explorer it will have a single backslashes and not double. It works with double as what you have in your example below but not with single as what I would end up passing. Do you have any other suggestions or methods to accomplish this?

Re: [R] noob requesting help

2012-06-18 Thread Rui Barradas
Hello, Ok, sorry about the misleading tip. This should do it. do.call(rbind, dep) Rui Barradas Em 18-06-2012 13:32, capital_P escreveu: Rui Barradas wrote Try names(departures) <- NULL dep <- data.frame(departures[ !sapply(departures, is.null) ]) You now have data.frame 'dep'. dep <

Re: [R] Loop Help

2012-06-18 Thread David Winsemius
On Jun 18, 2012, at 5:47 AM, @ngel wrote: I think that for now I'll keep it There is no "it" in this message. Users of Nabble often fail to read the POsitng Guide where it is requested that you include context. as is and try to clean the data in excel. Could someone please tell me how c

[R] Obtaining r-squared values from phylogenetic autoregression in ape

2012-06-18 Thread Louise Mair
Hello, I am trying to carry out a phylogenetic autoregression to test whether my data show a phylogenetic signal, but I keep calculating bizzare R-squared values. My script is: > library(ape) > x <- "R:1,GK:1)d:1,(MW:1,G:1)n:1)a:1,SPW:1)a:1,WB:1)b:1,(((SPBF:1,PBF:1)n:1,(HBF:1,SWF:1)c:1)

Re: [R] How do anova() and Anova(type="III") handle incomplete designs?

2012-06-18 Thread Justin Montemarano
Thanks for your response, John. That was helpful. I was using Type III from Anova() as a comparison to some results I had obtained JMP, which I've lost access to and have moved on to R, and I was confused by the error. Given that I do have a continuous covariate, the analyses are not likely comp

[R] Error w/ Start up Script

2012-06-18 Thread Lopez, Dan
I noticed about a couple weeks ago that my R start up script starting showing the below errors, following the line "Type 'q()' to quit R" Not sure what I did different and so far it hasn't caused major problems but kind of fear that it will. Can you please help? I also posted my Rprofile below.

[R] Need to append vector to all levels of nested list WITHOUT a loop

2012-06-18 Thread Johannes Reichl
Hi everyone, I have a list betaMoments with 2 levels, e.g. > beta1 = list( +m = 4, +v = 5 +) > beta2 = list( +m = 6, +v = 7 +) > > betaMoments = list() > betaMoments[[1]] = beta1 > betaMoments[[2]]

Re: [R] Package of EM and MI for IRT in R

2012-06-18 Thread R. Michael Weylandt
Possibly Amelia: http://cran.r-project.org/web/packages/Amelia/index.html Best, Michael On Mon, Jun 18, 2012 at 3:19 AM, Kamontip Srihaset wrote: > > > Dear all members, > > I am Phd. candidate student at Chulalongkorn U., Thailand. I am interested in > expectation maximization algorithm (EM) a

[R] Installation problems in Pander from Github

2012-06-18 Thread avinash barnwal
Hello all, I was trying to create nice tables in HTML on windows and then got to know about Pander. I have installed Pandoc, devtools, Rtools (which has same path as of R ) . I have also gone through the depends of Pander and fixed it all. Still i am not able to install and getting a conflicting e

Re: [R] checking the validity of the expressions

2012-06-18 Thread R. Michael Weylandt
What does any of this mean? Can you give a concrete example of the transformation you hope to accomplish? Michael On Mon, Jun 18, 2012 at 3:10 AM, arunkumar wrote: > hi > > I have a data and have to convert it into the given expression > > like > > data=1000:2000 and expression is exp() > >

[R] making svm work harder

2012-06-18 Thread luiceur
Hi all, I am working on an optimization of the svm classification function based on the package e1071. So far we have not been able to find any dataset that really slows down the application. However, from R users surveys seem that svm would be a good function to speed up. Is there any one working

Re: [R] out put of loop in list

2012-06-18 Thread Rui Barradas
Hello, Please read the posting guide: to draw samples from what? A vector? Matrix? Etc... Anyway, here it goes. In all cases the result is a list with 4 elements. # Draw 4 samples of integers from 1 to 10 x <- rep(10, 4) # without replacement set.seed(1) lapply(x, sample, size=6) # with rep

Re: [R] Error while installing R on RHEL 6

2012-06-18 Thread manish gupta
I used it but still same error. Regards On Mon, Jun 18, 2012 at 7:14 PM, Pascal Oettli wrote: > Hello, > > The error message gives you the solution. > > Regards > > > > Le 18/06/2012 18:26, Manish Gupta a écrit : > > Hi, >> >> I am trying to install R from source and using configure command bu

Re: [R] aligning axis labels in a colorkey from levelplot

2012-06-18 Thread Deepayan Sarkar
On Sat, Jun 16, 2012 at 1:41 PM, Stephen Eglen wrote: > R does a great job with the fine details regarding plots.  e.g in the > following: > > library(lattice) > y <- -4:4/10 > xyplot(y~1, las=1) > > the y axis is labelled with numbers -0.4, -0.2, 0.0, 0.2, 0.4 with the > numbers aligned on the de

Re: [R] out put of loop in list

2012-06-18 Thread jim holtman
I would be good if you showed us what you are doing now. There are ways of drawing samples and storing in a list, but we have to see what you are trying to accomplish. You could use 'lapply' or do it in a loop, but the solution would depend on how you are doing things now and what is the problem

Re: [R] Using Sub

2012-06-18 Thread Lopez, Dan
Hi Duncan, That's right but I want to be able to use this to replace back slashes with front slashes when I copy and paste file paths. How can I make this work so I don't get this error? Thanks. Dan -Original Message- From: Duncan Murdoch [mailto:murdoch.dun...@gmail.com] Sent: Friday

Re: [R] Efficient distance calculation on big matrix

2012-06-18 Thread Boel Brynedal
Stefan, that looks wonderful! I am most certainly going to try to download and use your 'wordspace', even though I am unsure on how to even download it at this point. Many thanks! But yes, I need the Canberra distance unfortunately. But decreasing the time by more than 60% will save me days! If I c

Re: [R] Loop Help

2012-06-18 Thread Jeff Newmiller
Then don't do that. write.csv is designed to write a complete matrix or data frame in one call. Combine your list of matrices into one (e.g. using sapply) and write that in one call afterward rather than piecewise as you do the calculations. -

Re: [R] Loop Help

2012-06-18 Thread @ngel
Hello Jim thank you so much for your response. When I use write.csv, on the csv file is saved only the last loop So when I write results.matrix <- apply(degree_w(net.static[[i]]), 2, mean) write.csv(results.matrix, file = "results.csv") instead of getting a 550 line document with all loop result

Re: [R] Loop Help

2012-06-18 Thread jim holtman
What do you mean by "only the last j"? write.csv writes to a file you specify and I assume that you are changing the file name if you are making successive writes. You need to provide more detail. On Mon, Jun 18, 2012 at 5:47 AM, @ngel wrote: > I think that for now I'll keep it as is and try to

Re: [R] How to specify "newdata" in a Cox-Modell with a time dependent interaction term?

2012-06-18 Thread Terry Therneau
I've been out for a week, with one more to go. I'll look at this in earnest when I return. Terry T On 06/17/2012 04:07 AM, Jürgen Biedermann wrote: > Dear John, > > Thank you very much for your help! It was very important for getting > along further. > > I found out some additional things whi

Re: [R] some help to improve "hist to plot relative frequencies"

2012-06-18 Thread Deepayan Sarkar
On Fri, Jun 15, 2012 at 9:22 PM, gianni lavaredo wrote: > Dear Researches, > > sorry for disturb. I wish to improve my figure in R plotting the relative > frequencies of my data set. > > library(lattice) > a <- c(0,0,0,1,1,2,4,5,6,7,7,7,7,7,8,8,8,8,9,9,9,9,10,10,11) > histogram(a, xlab="myData") >

Re: [R] noob requesting help

2012-06-18 Thread capital_P
Rui Barradas wrote > > Try > > names(departures) <- NULL > dep <- data.frame(departures[ !sapply(departures, is.null) ]) > > You now have data.frame 'dep'. > dep <- data.frame(departures[ !sapply(departures, is.null) ]) Error in data.frame(`121.1` = list(device_info_serial = integer(0), hour

Re: [R] noob requesting help

2012-06-18 Thread Rui Barradas
Hello, Have you tried what I wrote in my last post? Revised: # Make a copy first, and operate on the copy. dep <- departures names(dep) <- NULL dep <- data.frame(dep) Rui Barradas Em 18-06-2012 11:09, capital_P escreveu: David Winsemius wrote That's going to remove a lot of rows. Since y

Re: [R] triangular matrix

2012-06-18 Thread Miguel Manese
Hi Lucia, On Mon, Jun 18, 2012 at 6:11 PM, lucinka wrote: > Hello, > > I got this matrix of gentic distances between my samples. it is 85x85 but > only lower half (without diagonal) contains my distances. How can I make a > mean and standard deviation on these distances, please ? You can try so

Re: [R] Error while installing R on RHEL 6

2012-06-18 Thread Prof Brian Ripley
On 18/06/2012 10:26, Manish Gupta wrote: Hi, I am trying to install R from source and using configure command but gettting one error message. *configure: error: --with-readline=yes (default) and headers/libs are not available* I am installing on server which is not connected by net so can inst

  1   2   >