similar to: highlight overlapping region of two densities

Hi all, this is OT, but maybe someone can give me a clue. I've got data from eye tracker experiments (750 data points). These figures show how the data lock like |** * ***** |*** ****** | ** **** | * * | | * * | | * |-------------- |--------------- Y axes display velocity, x axes display time, * are the data points.

Are there any functions in R for fitting distributions? In particular, I would like to fit the weibull and extreme-value distributions when: i) the data are given and the maximum likelihood estimates are required; ii) the sample moments are available but the complete data record is not (so moment estimates are required). Any suggestions will be gratefully received.

Dear list, I have two things I am struggling... # First set.seed(123) myD <- data.frame( Light = sample(LETTERS[1:2], 10, replace=T), Feed = sample(letters[1:5], 20, replace=T), value=rnorm(20) ) # Mean for Light myD$meanLight <- unlist( lapply( myD$Light, function(x) mean( myD$value[myD$Light == x]) ) ) # Mean for Feed

Using: # rsync --version rsync version 3.0.7 protocol version 30 Copyright (C) 1996-2009 by Andrew Tridgell, Wayne Davison, and others. Web site: http://rsync.samba.org/ Capabilities: 64-bit files, 64-bit inums, 32-bit timestamps, 64-bit long ints, socketpairs, hardlinks, symlinks, IPv6, batchfiles, inplace, append, ACLs, xattrs, iconv, symtimes rsync comes with ABSOLUTELY NO

Hello everybody, Data is myd <- data.frame(id1=rep(c("a","b","c"),each=3),id2=rep(1:3,3),val=rnorm(9)) I want to get a cumulative sum over each of id1. trying aggregate does not work myd$pcum <- aggregate(myd[,c("val")],list(orig=myd$id1),cumsum) Please suggest a solution. In real the dataframe is huge so looping with for and subsetting is not a

Dear All, I have a database of 200 observations named myD. In the dataframe there are a column named code (with codes varying from 1 to 77), a column named "prevalence" with some quantitative measurements are given and an column named Pr_mean, with no values. I would like to set a cycle to compute the average of prevalence values for each different code and store the averages under the

An object of class `density' has `x' and `y' components, which are what you need to do the plot. Just bind them into a matrix and use write() to write to a file. HTH, Andy -----Original Message----- From: Hinnerk Boriss [mailto:boriss at izbi.uni-leipzig.de] Sent: Thursday, November 14, 2002 7:12 AM To: R-help at stat.math.ethz.ch Subject: [R] evaluating density objects Hi! Is

Hi all, ok, i know i can cut a dendrogram, which i did. all i get is three objects that a dendrograms itself. for example: myd$upper, myd$lower[[1]], myd$lower[[2]] and so on. of course i can plot them seperately now. but the lower parts still have hundreds of branches. i?ll need a 30 " widescreen to watch the whole picture. what i?d like to is group the lower branches , so that i get a

Hi, I must be missing something here...Essentially, a short piece of code works if it's standalone, but doesn't work if it's divided into two functions. The code that works is: ################### WORKS ############### library(pamr) set.seed(120) x <- matrix(rnorm(1000*20),ncol=20) y <- sample(c(1:4),size=20,replace=TRUE) mydata <- list(x=x,y=y)

Hi Everyone, This example code results in R 'crashing'; that is the R application closes with no warnings or error messages. #----------------------- myD <- read.table(stdin(), header=TRUE, nrows=20) Broth Salt pH Temp N Y Growth 1 310 9.0 2.92 10 90.0 NA 0 2 615 6.0 7.82 30 1.0 2 1 3 217 2.0 7.34 10 7.0 8

At office, using the internal LAN at my disposal, I'm having a go at parallel computing - to begin with - with pvm, rpvm & snow. The two boxes are as follows Remote machine uffbsd: CPU: Intel(R) Pentium(R) 4 CPU 2.00GHz (1994.13-MHz 686-class CPU) Origin = "GenuineIntel" Id = 0xf24 Stepping = 4 real memory = 260046848 (248 MB) This machine NbBSD: CPU: Mobile Intel(R)

Hi, > tmpdata<-subset(myd, TYPE=="A") > levels(tmpdata$TYPE) > [1] "A" "B" "C" I'd like to get only "A" as output... Thanks for your help Marc

Hi, R is very new for me, so excuse if my questions are too basic... BTW, are there any forum where new R users could get help without annoying this huge mailing list ? In following code, I'd like to choose the color for each of the curve diplayed. png(filename = filename, width = 950, height = 600, pointsize = 10, bg = "white") xyplot(HITS+MS1*3+FREQ~TIME

Dear friends, Do you know how to calculate the CONCAVE hull of a set of points (2- dimensional or n-dimensional)? is that possible in R? (With a "smoothing" parameter of course). Best, -- Corrado Topi Global Climate Change & Biodiversity Indicators Area 18,Department of Biology University of York, York, YO10 5YW, UK Phone: + 44 (0) 1904 328645, E-mail: ct529 at york.ac.uk

sorry for cross-posting Dear all, I have tow (several) bivariate distributions with a known mean and variance-covariance structure (hence a known density function) that I would like to compare in order to get an intersect that tells me something about "how different" these distributions are (as t-statistics for univariate distributions). In order to visualize what I mean hear a little

Dear all.... any one can help me??? I am aware of monitoring disk space using nagios (disk partitions with check_disk). Now, I have a situation that I would like to monitor size of an individual file through nagios. e.g /var/lib/mysql/somedatabase/vlog.myd, this file will grow very fast and i would like to monitor if the size of this file crosses 1.5 G. I'm already used 'check_file'

https://bugzilla.samba.org/show_bug.cgi?id=8712 Summary: --link-dest doesn't work if target file exists (but needs updating) Product: rsync Version: 3.0.7 Platform: All OS/Version: Linux Status: NEW Severity: normal Priority: P5 Component: core AssignedTo: wayned at

Hello group, I'm having a very hard time rsyncing efficiently a MySQL database which contains very large binary blobs. (Actually, it's the database of Mantis bug tracker [http://www.mantisbt.org/], with file attachments stored directly in the table rows. I know it's a bad idea from many other reasons, but let's say it was given to me as such.) First, I was dumping the

I''ve "inherited" the administration of a puppet-dashboard (version 1.1.0, installed on RHEL 5.6 from puppetlabs RPM), and have hit a problem I''m hoping for some help with. In short, one table, "resource_statuses" appears to be growing at a rate far higher then the other tables: mysql> select count(*) from nodes; +----------+ | count(*) | +----------+ |

Hi, We're encountering a rather weird, sporadic hang during a large-file transfer. A bit of background on the setup: - Host A rsyncs a MySQL database from Host B. Both sides use rsync 3.1.1 on Debian Jessie. - Host B serves the files using rsync --daemon from an ext4 filesystem residing on a thin LVM snapshot. - Host A uses the following command to rsync the files to a btrfs