similar to: Estimation of discrete unimodal density

Dear R-users Do you know of an easy way in R of performing "isotonic (unimodal)" regression ala BBBB, Barlow, Bartolomew, Brenner and Brunk. best regards Helgi -- Helgi Tomasson FAX: 354-552-6806 University of Iceland PHONE:354-525-4571 Faculty of Economics and Business Administration email:helgito at

Dear All, Does anybody know if there is a functionality in R to break histograms that show a clear bi-modal (or multi-modal) distribution into a series of unimodal histograms that added up result in the original histogram? I was thinking of using QQ-plots (for which tools are available in R), and then observing the number of times the observed quantiles cross the 1:1 line, but this only gives an

Hi. I hope you can help me... I have fitted the following ARIMA model: arima1<-arima(bigspring$log.volume, order=c(0,1,2)) I need to predict 30 days ahead. I used following code predict(arima1,n.ahead=30,se=T) However I get 30 predictions, but from predictions 2:30 I get the same predictions. Why is this? What am I doing wrong Thanks Catherine KSS Ltd Seventh Floor St

Hola a tod en s, estoy intentando averiguar el número de modas en una distribución. Para ello utilizo diptest. Mi duda es que no acabo de entender cuando la información suministrada por los test suponen la existencia o no de unimodalidad/multimodalidad. Una parte de la salidad de diptest es la que pego a continuación (el resto esta en el fichero adjunto con las distribuciones kernels y las

Dear all i could not estimate the optima value or range value in unimodal plot in glm please help me out thanking you regard madan _________________________________________________________________ Hotmail® has a new way to see what's up with your friends. orial_WhatsNew1_052009 [[alternative HTML version deleted]]

Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: 7bit X-Virus-Scanned: by amavisd-new X-Spam-Checker-Version: SpamAssassin 2.63 (2004-01-11) on hypatia.math.ethz.ch X-Spam-Level: **** X-Spam-Status: No, hits=4.2 required=5.0 tests=MSGID_FROM_MTA_HEADER,RCVD_IN_BL_SPAMCOP_NET autolearn=no version=2.63 Dear r-help list, my question is about ordination technics:

Hi, I'm trying to find an implementation of monotonic regression in R and I haven't been able to find anything that's really related to this. isoMDS in the MASS package uses monotonic regression, however, I was wondering if there is any standalone function for monotonic regression? Basically what I'm trying to do is implement monotonic regression where I can see not just the

Hi List, I'm currently setting up an openldap server and included the following lines in my slapd.conf : modulepath /usr/lib/ldap moduleload back_hdb after finishing up my config and i run slaptest on it i get an error saying that the modulepath doesn't exist. I checked and it indeed isn't there , in fact i can find it anywhere on my system (centos 5.7). the packages i've

I am interested in developing a backend for the PIC microcontrollers. Specifically, I plan to write for the 8-bit MCUs, but would hopefully also support the 16 and 32 bit models (eventually). It is my understanding that a backend for PIC16 was previously a part of the project but was dropped in version 2.9. Is there any development on this front currently? Is there any interest in developing this

I am attempting to fit a full zero-inflated Poisson model then use backward elimination to arrive at the best-fitting model. When I try to fit the model with zeroinfl() I get this error: Error in while (abs((ll_old - ll_new)/ll_old) > control$reltol) { : missing value where TRUE/FALSE needed Any suggestions? Thanks for your help! Paige Barlow MS Student Virginia Tech Dept Fish

Hi! I'm trying to transfer a incomming call from a Queue to another extension. I'm receiving a call from a queue with the AgentCallbackLogin. The queu is as following: Queue(sales|t) Which should allow transfers. So as soon as the call is answered I would like to be able to transfer it When the agent presses the # I get the dialtone but as soon as I press any digit Asterisk tells me

Hi All Currently i'm setting up a 5.4 server and try to create a 3rd raid device, when i run: $mdadm --create /dev/md2 -v --raid-devices=15 --chunk=32 --level=raid6 /dev/sdc /dev/sdd /dev/sde /dev/sdf /dev/sdg /dev/sdh /dev/sdi /dev/sdj /dev/sdk /dev/sdl /dev/sdm /dev/sdn /dev/sdo /dev/sdp /dev/sdq the device file "md2" is created and the raid is being configured. but somehow

Hi All, I've been asked to setup a 3d renderfarm at our office , at the start it will contain about 8 nodes but it should be build at growth. now the setup i had in mind is as following: All the data is already stored on a StorNext SAN filesystem (quantum ) this should be mounted on a centos server trough fiber optics , which in its turn shares the FS over NFS to all the rendernodes

Hi all, I?ve got a HFS+(not journaled) volume connected to my centos6.2 test server, i installed the kmod-hfs(plus) packages and read/write works all fine. but since i?m going to use this for serving mac home folders via netatalk i would like to mount it with support for Extended Attributes and acl?s. so i add user_xattr and acl to my fstab options but then it fails to mount. checking the

Hi all, I've setup a ethernet bond on my centos 5.6 server , when i do a reboot the bond does come up but cleared all the slaves and i've to manually re-add them with ifenslave. does anyone know a solution to this? am i missing something? offcourse i can add it to my rc.local but there must be a more elegant way. please see my configs below Thanks, Wessel ifcfg-bond0: DEVICE=bond0

Hi All, I'm currently trying to install centos 5.4 x86-64bit on a raid 1, so if one the 2 disks fails the server will still be available. i installed grub on /dev/sda using the advanced grub configuration option during the install. after the install is done i boot in linux rescue mode , chroot the filesystem and copy grub to both drives using: grub>root (hd0,0) grub>setup (hd0)

Dear List, I have a set of distributions recorded at an equal interval of time and I would like to plot them as series of horizontal histograms (with the x-axis representing time, and y-axis representing the bins) since the distribution shifts from unimodal to multimodal in several occasions. What I would like to see is something close to a violinplot, but I do not want a kernel density

Hi All, I am in the middle of doing an analysis of a Case-Cohort design. I had three questions about the analysis: a) Does any one know of some public code for developing the patient risk sets (indexed by failure time) or is there a better way to organize the data? b) I was planning to use the Barlow weighting method. Has this or any other weighting method (Prentice, Self-Prentice) been

I am just beginning to learn R. Please help. I have two matrices. One 20 columns by 16 rows represent student responses on a 20 question test. The second is 20 by 16 representing the correct responses 16 times. Subtracting one from the other should yield 0's if a response is correct. Yet Response minus Bigkey doesn't work. Documents indicate vectors should add and subtract, and

I'm not sure where to begin with this, but I was wondering if someone could refer me to an R package that would test to see if a distribution fits a bimodal distribution better than a unimodal distribution. Thanks, Andrew [[alternative HTML version deleted]]