similar to: Puzzled by what Rprof is telling me

Hi I am trying to print out means, STDs and histograms under two sets of factors. I can manage it for one set - see below but not for two sets. That is, I want ot print out the mean STD and Histogram for each ITEM code within each DELIVERABLE code. In addition I can only get to view the histogram for the last item. How do you get R to stop overriding the histogram for eg level 1 for factor 1

Dear members, I would like a help for extracting the values from a step function (stepfun). >From help(stepfun) we have the following example: Y0<-c(1.,2.,4.,3.) y0<-c(1.,2.,3.,4.) sfun<-stepfun(1:3,y0,f=0) plot(sfun) Now, suppose instead I was given the object (*sfun*, say) from which I wanted to extract the values generated by the function *stepfun*. More precisely, I want to

I have a dataframe from an On-Board Survey with weights ("expwgt") and variables for up to 8 used lines: VEH1 through VEH8. The lines are labeled "MT-..1" through "MT-902". I want to know how many transfers there are between MT-802 and MT-901. That is, when one of them is VEHx and the other is VEHx+1 or VEHx-1 E.g. {VEH1 = MT-802 AND VEH2 = MT-901 } plus {VEH2 =

I have several data input functions that include a construct like the following: indat <- function(stuff) { [...] awkp <- paste("awk ...yada yada...", stuff, etc.) scan(pipe(awkp),quiet=TRUE) } I use these functions a lot, and as I do, a lengthening list of closed, but undestroyed, connections accumulates. Ultimately this appears to stuff things up and cause

Hello everyone, I have been trying to figure out away to integrate under a spline produced by the package tps(fields). As the package does not output functions I am trying to do something similar to the trapezium rule. My data are 3D (x, y & z). I have extracted from the surface output by Tps the values of z at regular intervals so that I have a grid of figures, for example: 1 4 6 6 8 8

I've been experimenting with R during the last week and I thought I'd pass on some comments on issues on the documentation that left me a bit confused. To start with (literally), the ReadMe.rw1021 file says that the easiest way to install R is to download the rwinst.exe installer program. That gave me the impression that the program would connect to an R download site, download the

I have 2 dataframes. "mydata" contains numerical data. "mybys" contains information on the "group" each row of the data is in. I wish to aggregate each column in mydata using the corresponding column in mybys. Please see the example below. What is a more elegant or "better" way to accomplish this task? Thanks! mydata =

Hi All I need some help with plotting a step function, currently I'm using sfun <- stepfun(c(1, 2, 5,10, 20), c(0, 11, 22, 33, 44, 0), f=0) plot(sfun, pch=NA, main="", xlim=c(1,20)) which I working fine, but my data is in the following format: Min Max Value 1 2 11 2 5 22 510 33 10 20 44 1. To save time, is there a way to directly plot the above data as it is (without any

Hi all, I have a dataset. Each time I want to sample N(i) elements from it and I want to repeated sampling M times. N(i) is varied from time to time. For example, dataset = 1:50; a = list(); M = 1000; I want to do something like for(i in 1:M) { a[[i]] = sample(dataset,

Hi, I have been examining large data and need to do simple linear regression with the data which is grouped based on the values of a particular attribute. For instance, consider three columns : ID, x, y, and I need to regress x on y for each distinct value of ID. Specifically, for the set of data corresponding to each of the 4 values of ID (76,111,121,168) in the below data, I should invoke

First, I love R and am grateful to be using this free and extremely high quality software. Recently I have been comparing two algorithms and naturally I programmed in R first. It is so slow that I can almost feel its pain. So I decided to do a comparison with Java. To draw 500,0000 truncated normal, Java program takes 2 second and R takes 72 seconds. Not a computer science major, I find it hard

Please note is with replacement From: taby gathoni <tabieg@yahoo.com> To: R help <r-help@r-project.org> Sent: Thursday, March 10, 2011 11:53 AM Subject: [R] random sampling steps in R Dear all, Could someone assist me in random sampling steps/code in R? I have a main sample of 42 males and 165 females and I want to come up with about 1000 samples of 20 males and 20 females from

This commit breaks loading nouveau on my optimus laptop. I get the following spam in dmesg: [ 24.102587] MXM: GUID detected in BIOS [ 24.103167] ACPI Exception: AE_AML_PACKAGE_LIMIT, Index (0x000000003) is beyond end of object (length 0x0) (20131218/exoparg2-420) [ 24.103242] ACPI Error: Method parse/execution failed [\_SB_.PCI0.VGA_._DSM] (Node ffff88007c82e1e0), AE_AML_PACKAGE_LIMIT

Hello, Given a vector I would like to rapidly identify duplicated non contiguous elements. Given for example c(1, 1, 2, 3, 2, 4, 5, 6, 4) I would like to get: FALSE FALSE TRUE FALSE TRUE TRUE FALSE FALSE TRUE In fact I need to check this on the columns of a matrix! I can do that of couse with loops but is there any function already available? Thanks

Dear list, I observed two problems that I suppose are generic. First, using 'Rprof' to profile a parallel (based on the package 'parallel') code caused .... Error in unserialize(node$con) : error reading from connection.... Second, on a multicore desktop, I concurrently opened two terminals and ran two separate R sessions, both running (actually identical) parallel code (which

Full_Name: Kevin C. Bartz Version: 1.5.1 OS: Solaris 2.6 Submission from: (NULL) (192.223.226.5) A while ago, I noticed this oddity about R profiling and setMethod. First, I "test out" Rprof. > require(methods) Loading required package: methods [1] TRUE > > Rprof("test.out") > data.frame("a") X.a. 1 a > Rprof(NULL) So far, so good. Next, I

I'm having difficulty with Rprof. I have a documentation example test.qss that runs fine without profiling, but under Rprof, > Rprof() > source("test.qss") Error in standardGeneric("model.matrix") : UseMethod used in an inappropriate fashion Luke wrote about a similar circumstance last summer: # From: Luke Tierney (luke@stat.umn.edu) # Date: Fri Jul

I find it difficult to find Rprof using the usual search tools... could I suggest that it be cross-referenced to profile Rprofile Profile to make this easier? url: www.econ.uiuc.edu/~roger/my.html Roger Koenker email rkoenker@uiuc.edu Department of Economics vox: 217-333-4558 University of Illinois fax: 217-244-6678 Champaign, IL 61820

Full_Name: John Garvin Version: 1.9.0 OS: Linux Submission from: (NULL) (128.42.129.78) This may or may not technically be a bug, but it's certainly an annoyance. Rprof only takes into account computation that occurs inside functions. If a time-consuming operation occurs outside a function, it doesn't record the time it takes. Consider this program 'array.r': Rprof() foo <-

Hi, >From reading ?Rprof, I checked R CMD Rprof --help and learned that there are options to specify the min % to print. This is currently (R-devel r52975) displayed with the --help option as --min%total minimum % to print for 'by total' --min%self minimum % to print for 'by self' So I tried R CMD Rprof --min%total 5 and got an error. After looking at