similar to: integration of a discrete function

Hello list! I would like to fit a distribution to each of the peaks in a histogram, such as this: http://photos1.blogger.com/blogger/7029/2724/1600/DU145-Bax3-Bcl-xL.png . The peaks are identified using Petr Pikal peaks function ( http://finzi.psych.upenn.edu/R/Rhelp02a/archive/33097.html), but after that I am quite stuck. Any idea as to how I can: Fit a distribution to each peak Integrate the

Dear R users, I have spectral data (say, wavelength vs. extinction coefficient) for which I??d like to calculate an integral (i.e. the area underneath the curve). Suppose the (artificial) dataset is lambda E 1 2 2 4 3 5 4 8 5 1 6 5 7 4 8 9 9 8 10 2 How can I calculate an integral for these values using R? Many thanks for any help! Regards Christoph

Dear all I would like to have some data in the form of a 2-dimensional array (matrix) of lists, so that I can easily find the desired list object by indexing the structure by rows and columns. In matlab there exists a data type called "cell array": a matrix of "cells", which are composite objects very similar to R lists. I know that in R you can create 1-dimensional

Dear R users I have a really simple question (hoping for a really simple answer :-): Having estimated the spectral density of a time series "x" (heart rate data) with: x.pgram <- spectrum(x,method="pgram") I would like to compute the power in a specific energy band. Assuming that frequency(x)=4 (Hz), and that I am interested in the band between f1 and f2, is the

Hi all I am running a 2-way repeated measure anova with 1 between-subjects factor (Group=treatment, control), and 1 within-subject factor (Time of measurement: time1, time2). I extract the results of the anova with: summary(aov(effect ~ Group*Time + Error=Subj/Time, data=mydata)) Now, this must be clearly a dumb question, but how can I quickly extract in R all the post-hoc t-tests for the

Dear all, I am unsure about how to specify a model in R and I thought of asking some advice to the list. I have two groups ("Group"= A, B) of subjects, with each subject undertaking a test before and after a certain treatment ("Time"= pre, post). Additionally, I want to enter the age of the subject as a covariate (the performance on the test is affected by age),

Hi all I am running a 2-way repeated measure anova with 1 between-subjects factor (Group=treatment, control), and 1 within-subject factor (Time of measurement: time1, time2). I extract the results of the anova with: summary(aov(effect ~ Group*Time + Error=Subj/Time, data=mydata)) Now, this must be clearly a dumb question, but how can I quickly extract in R all the post-hoc t-tests for the

Hello, I am trying to plot a few correlograms on the same figure, with the function corrgram() from the package corrgram. However, the function does not seem to use the base graphic system, as setting out the multiple figure layout with, e.g., par(mfrow=c(2, 2,)) does not work. Does anybody know a workaround for this? Many thanks in advance for any advice best giuseppe -- Giuseppe Pagnoni,

Dear all, I am estimating a mixed-model in Ubuntu Raring (13.04¸ amd64), with the code: fm0 <- lme(rt ~ run + group * stim * cond, random=list( subj=pdSymm(~ 1 + run), subj=pdSymm(~ 0 + stim)), data=mydat1) When I check the approximate variance-covariance matrix, I get: > fm0$apVar [1] "Non-positive definite

Dear all, I am having quite a hard time in trying to figure out how to correctly spell out a model in R (a repeated-measures anova with a within-subject covariate, I guess). Even though I have read in the posting guide that statistical advice may or may not get an answer on this list, I decided to try it anyway, hoping not to incur in somebody's ire for misusing the tool. For the sake of

An embedded and charset-unspecified text was scrubbed... Name: inte tillg?nglig URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20110610/b259a660/attachment.pl>

Full_Name: Morten Welinder Version: 2 OS: Solaris/space/gcc2.95.2 Submission from: (NULL) (65.213.85.217) I changed src/nmath/standalone/test.c to read: --------------------------------------------------------------------------------- #define MATHLIB_STANDALONE 1 #include <Rmath.h> #include <stdio.h> int main() { double x; for (x = 99990; x <= 100009; x++) printf

Hi, Please I have a question: If my server is out of space (disk full) and the postfix is LDA, the new messages a deffered or refused,so postfix send 450 to sender ... But When I use dovecot (1.2.16) like LDA, postfix accept message (even there is no space), forward to dovecot, dovecot REFUSE the message and message is lost ! How can I solve this ? I don't want tu use quota and the disk

Is there some sort of replot function, i.e., reproducing the current plot like gnuplot's command replot? Of course it's easy to issue a plot command twice, but after several lines statements it gets a little bit annoying ... and the plot disappears each time one resizes the window, changes from single plot to multiple plots etc. Just dreaming ... Fritz --

_ platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 1 minor 6.2 year 2003 month 01 day 10 language R > Until this version, I've not had to explicitly start the x11 device. Now,

I thought this was going to be easy ... Can the label size of identify() be controlled by setting par(cex.*) because I'm having no luck? My only recourse is to save the index and position of the labels from identify() and use text() to replot them. Regards Alex Alex Hanke Department of Fisheries and Oceans St. Andrews Biological Station 531 Brandy Cove Road St. Andrews, NB Canada E5B 2L9

Greetings, I am trying to get samba3 to work well with W95. The windows 95 box can see the samba server, but cannot authenticate (claims wrong password) to get to any shares. I know the smbpasswd is good as I can authenticate and get to the shares using the same account and password on a W2000 machine. I suspect I am doing something wrong in the [global] section of the samba server. Would

Hello, i am working in origin pro, i want to plot a graph as like a pdf attached but with black and white lines. here radial axis varies from 0 to 1. and angular axis from 0 degree to 60 degree.and third axis which is depend on both radial and axial gives non intersecting lines. how can i read the data from plot for replot. vikas -------------- next part -------------- A non-text attachment was

Hi I want to create a vector of colors that are as different from one another as possible. ?rainbow states "Conceptually, all of these functions actually use (parts of) a line cut out of the 3-dimensional color space...". This suggests to me that the resulting colors are all placed on this "line" and are equi-distant along it. The resulting color palette is a range of

Sorry, Should have used plain text to begin with and better wrapping: Hi People, Trying to be as informative as possible: I'm trying to add printers drivers to cups-owned, samba shared printer. Called "PDF-Creator" But when I issue the command "cupsaddsmb -U root -v PDF-Creator" I get a "NT_STATUS_LOGON_FAILURE" Weird thing is: I was able to do this a few