similar to: integral depending on a parameter

Hi all, I have the data from the posterior distribution for some parameter. I want to find the 95% credible interval. I think "t.test(data)" is only for the confidence interval. I did not fine function for the Bayesian credible interval. Could some one suggest me? Thanks [[alternative HTML version deleted]]

Would like to know how much data can R process - number of rows and columns? [[alternative HTML version deleted]]

Dear R useRs, i have the function f1(x, y, z) which i want to integrate for x and y. On the one hand i do this by first integrating for x and then for y, on the other hand i do this the other way round and i wondering why i doesn't get the same result each way? z <- c(80, 20, 40, 30) "f1" <- function(x, y, z) {dgamma(cumsum(z)[-length(z)], shape=x, rate=y)}

I have been told by the CRAN administrators that the following code generated an error on 64-bit Fedora Linux (gcc, clang) and on Solaris machines (sparc, x86), but runs well on all other systems): > fn <- function(x, y) ifelse(x^2 + y^2 <= 1, 1 - x^2 - y^2, 0) > tol <- 1.5e-8 > fy <- function(x) integrate(function(y) fn(x, y), 0, 1,

Dear List I apologize for the multiple postings. After being in the weeds on this problem for a while I think my original post may have been a little cryptic. I think I can be clearer. Essentially, I need the following a <- c(2,3) b <- c(4,5,6) (2*4) + (2*5) + (2*6) + (3*4) + (3*5) +(3*6) But I do not know of a built in function that would do this. Any suggestions? -----Original

The integral of any probability density from -Inf to Inf should equal 1, correct? I don't understand last result below. > integrate(function(x) dnorm(x, 0,1), -Inf, Inf) 1 with absolute error < 9.4e-05 > integrate(function(x) dnorm(x, 100,10), -Inf, Inf) 1 with absolute error < 0.00012 > integrate(function(x) dnorm(x, 500,50), -Inf, Inf) 8.410947e-11 with absolute error <

I can confirm git head builds and passes /make fullcheck/ (skipping the noise part of test_streams.sh) on - armv6-linux (Raspbian), GCC 4.6 & GCC 4.8 - i686-pc-mingw32, GCC 4.8.1 - MSVC 2005 Express 32-bit (make fullcheck with MinGW) - MSVC 2013 Express 32-bit (make fullcheck with MinGW) Furthermore, I've also ran the test_streams.sh of FLAC 1.2.1 with FLAC 1.2.1 as encoder and git as

Has anyone successfully setup a shorewall Ucarp solution?

Dear r-users, Is there any command in R allowing to evaluate a double integral? for instance let say I want to evaluate the following integral: integrate[lo=(0,1),up=(2,3)] f(x,y)=x^2+y^2 where lo is the vector of lower bounds and up that of upper bounds. I thaught the function "adapt" would work but it did not. Many thanks, Dominique K.

Good day, all, I've gleaned that dovecot has support for verifying SSL client certs on connect. What IMAP clients have people found that support sending a specific SSL cert? I can't find this feature in Thunderbird, at least. Thanks, Peter

Hello everyone, I'm trying to determine if there is an equivalent to the WinXP/Vista/ Server command SC.exe (service control). From what I have gleaned, the svcctl API in Samba tries to replicate the other direction (controlling a Samba server's services/inits from Windows/MMC) but I'm looking to control a Windows server's services from (in my case) Linux. Essentially, I'm

Hello, I hope it's the correct ML to get support for "advanced" ssh use (sorry if it's not the case) And I would be very grateful if someone could help me on this issue. Here is my challenge : - I have X devices (around 30) and one SSH server - Each of them have a unique public key and create one dynamic reverse port forwarding on the server - All of them connect with the

Hi, > I was thinking of checking whether pdev is a VGA class device and whether > it decodes vga access, and in that case automatically calling How can I figure that? Ok, class is easy, but decode? pci.h offers functions to set vga decode but not to get that info ... thanks, Gerd

Hi, > I was thinking of checking whether pdev is a VGA class device and whether > it decodes vga access, and in that case automatically calling How can I figure that? Ok, class is easy, but decode? pci.h offers functions to set vga decode but not to get that info ... thanks, Gerd

On Tue, Sep 1, 2015 at 10:57 AM, peter dalgaard <pdalgd at gmail.com> wrote: > > On 01 Sep 2015, at 16:21 , Avraham Adler <avraham.adler at gmail.com> wrote: > >> Hello. >> >> There used to be changelog of sorts for R-devel [1] and R-release [2]. >> Neither have been updated since 2015-07-24. Have these moved >> elsewhere, or are they no longer

To enhance my understanding, and that of my students, I have a question about cox.zph in the survival package. If I have correctly gleaned the high-level point from the 1994 Biometrika paper of Grambsch and Therneau, it looks to me like cox.zph provides a mechanism to test for a simple trend in plots of a function of time, g(t) versus the scaled schoenfeld residuals and it also provides some

On 11/21/2015 12:57 PM, Marcelo Roccasalva wrote: > Boot from DVD/USB. Select rescue installed (on centos7 or vmlinuz > rescue on older), select configure network and mount installed, then > "chroot /mnt/sysimage" and yum install kernel... This should do... > > On Sat, Nov 21, 2015 at 2:05 PM, Siva Prasad Nath > <shivaprasadnath21 at gmail.com> wrote: >> My

Hi everyone, Ok, so I admit I'm asking this to avoid having to do my own research. but I'm also asking in the hopes of avoiding doing my own scripting. ACB Radio is looking at moving across to Icecast. There's really not a lot to stop us at this point. But we are keen to retain access to total listener hours, average listen time and amount of data served per mountpoint. I

> > I have checked in all the latest into CVS and am going to start the > > test suite again. if all goes well I will probably release this as > > beta 10. > > > > anyway, try it out and let me know if anything bad happens! it > > should be a short jump from beta 10 to 1.0. > > I've just checked out the latest from scratch. There is no configure

Dear list members, I am trying to estimate parameters of the AR(1)-GARCH(1,1) model. I have one additional dummy variable for the AR(1) part. First I wanted to do it using garchFit function (everything would be then estimated in one step) however in the fGarch library I didn't find a way to include an additional variable. That would be the formula but, as said, I think it is impossible to add