similar to: Fast trace of inverse

Dear R-users I would like to integrate something like \int_k^\infty (1 - F(x)) dx, where F(.) is a cumulative distribution function. As mentioned in the "integrate" help-page: integrate(dnorm,0,20000) ## fails on many systems. This does not happen for an adaptive Simpson or Lobatto quadrature (cf. Matlab). Even though I am hardly familiar with numerical integration the implementation

Hi, I am using SAM (from siggenes_1.2.11 package) method to select genes from a microarray data set. After installing the latest R2.4.0 on my computer, to my surprise the results are totally different from that calculated using R2.1.1. Even the example code doesn't work the same way under these two versions of R. Does anybody know what is going on? Thanks for any suggestions.

Dear R-users, I have one question about using SVD to get an inverse matrix of covariance matrix Sometimes I met many singular values d are close to 0: look this example $d [1] 4.178853e+00 2.722005e+00 2.139863e+00 1.867628e+00 1.588967e+00 [6] 1.401554e+00 1.256964e+00 1.185750e+00 1.060692e+00 9.932592e-01 [11] 9.412768e-01 8.530497e-01 8.211395e-01 8.077817e-01 7.706618e-01 [16]

Hi, i'm running CentOS 6.0 on my server and installed mrtg from the rpm-package mrtg-2.16.2 . I also installed the depending packages perl-IO-Socket-INET6 perl-Socket6 .... mrtg works fine with IPV4-Addresses. When i specify a Target by IPV6-Address (or hostname resolving to a V6-address) mrtg fails. Here i have a small sample-config for V4 which is working: LogDir: /tmp ThreshDir: /tmp

I am running bsnmpd with basic snmpd.config (only community and location changed). When there is a problem with HDD and disk disapeared from ATA channel (eg.: disc physically removed) the bsnmpd always dumps core: kernel: pid 1188 (bsnmpd), uid 0: exited on signal 11 (core dumped) I see this for a long rime on all releases of 7.x and 8.x branches (i386 and amd64). I did not tested 9.x. Is

https://bugs.freedesktop.org/show_bug.cgi?id=98630 Bug ID: 98630 Summary: Bad color and glitches on PPC64 Product: Mesa Version: 13.0 Hardware: Other OS: All Status: NEW Severity: normal Priority: medium Component: Drivers/DRI/nouveau Assignee: nouveau at lists.freedesktop.org

Hi all, Does anybody know any function that performs gaussian adapative quadrature integration of univariate functions? Thanks in advance, Regards, Caio __________________________________________________ [[alternative HTML version deleted]]

Guys I was so frustrated about installing samba with ADS and winbind support on solaris 8. After fixing many problems, I am now facing another problem: I created a directory "test" in the samba share and was editing a file under the directory "test" and in the meantime the machine was rebooted. After the machine is started, the whole directory "test" is gone. I

Hello, interesting to hear. What version of NUT, what way of installation on Windows, what sort of UPS (USB or some other)? Asking because this info is lacking and currently there are a few variants to consider: * an MSI package based on NUT v2.6.5 (last released about 9 years ago as "2.6.5-6"), * results of re-integration of that codebase into current NUT master (some months after

Hi, between March 16 and now, bge on a Sun X2200 stopped working, turning off msi (via hw..pci.enable_msi=0) got it working again. I tried first replacing bge with an older version but that did not help. please advice :-) Danny

If I try the breakpoints() function (strucchange package) with a minimum segment size = the number of regressors, there appears the following error message: "minimum segment size must be greater than the number of regressors" According to the documentation: "breakpoints implements the algorithm described in Bai & Perron (2003) for simultaneous estimation of multiple

Hi Everyone, I am trying to figure out the spatstat package for the first time and am having some trouble. Unfortunately, I can't post my data set but I'll hopefully post enough details for some help. I want to model the intensity of a spatial point process using 2 covariates from my data. After reading through the documentation, I have successfully created 2 "ppp" objects. The

Can anyone help me on how to get the nodes and weights of the adaptive quadrature using R. Best wishes Boikanyo. ----- The University of Glasgow, charity number SC004401

https://bugs.freedesktop.org/show_bug.cgi?id=98712 Bug ID: 98712 Summary: [PPC64] GPU Lock Up when Nouveau 2D acceleration is enabled + X failures Product: xorg Version: 7.7 (2012.06) Hardware: Other OS: All Status: NEW Severity: normal Priority: medium Component:

Hi all, We know that Hermite polynomial is for Gaussian, Laguerre polynomial for Exponential distribution, Legendre polynomial for uniform distribution, Jacobi polynomial for Beta distribution. Does anyone know which kind of polynomial deals with the log-normal, Student抯 t, Inverse gamma and Fisher抯 F distribution? Thank you in advance! David [[alternative HTML version deleted]]

Hi Guys I compiled Samba 3.0.26a --with-ads and --with-winbind with krb5-1.6.2 and openldap-2.3.38 using gcc version 3.4.6 on Solaris 8. and added winbind in passwd entry in /etc/nsswitch.conf passwd: files winbind The problem is the entry "passwd: files winbind" in /etc/nsswitch.conf caused "ls -l" and Samba hang. If winbind is removed from nsswitch.conf, "ls

I've written a series of functions that evaluates an integral from -inf to a or b to +inf using equally spaced quadrature points along a normal distribution from -10 to +10 moving in increments of .01. These functions are working and give very good approximations, but I think they are computationally wasteful as I am evaluating the function at *many* points. Instead, I would prefer to use

Hi, I would like to solve a double integral of the form \int_0^1 \int_0^1 x*y dx dy using Gauss Quadrature. I know that I can use R's integrate function to calculate it: integrate(function(y) { sapply(y, function(y) { integrate(function(x) x*y, 0, 1)$value }) }, 0, 1) but I would like to use Gauss Quadrature to do it. I have written the following code (using R's statmod package)

Together with Jan Sjodin (in copy of this email), we begin an implementation of the JIT with MC. The idea, suggested by Jan, is to develop a MCJIT in parallel of the current JIT and to keep the two implementations until (at least) the new MC one is mature enough. Currently code is kept on gitorious (http://gitorious.org/llvm-mc-jit/llvm-mc-jit). Following this, a boolean "bool MCJIT =

Hi All, I am trying to use A Gaussian quadrature over the interval (-infty,infty) with weighting function W(x)=exp(-(x-mu)^2/sigma) to estimate an integral. Is there a way to do it in R? Is there a function already implemented which uses such weighting function. I have been searching in the statmode package and I found the function "gauss.quad(100, kind="hermite")" which uses