similar to: optimize function with integral form ?

Hi List, I was not able to make this work. I know it is a simple one, sorry to bother. Give me some hints pls. Thanks! Jen if(length(real.d)>=30 && length(real.b)>=30 && beta1*beta2*theta1*theta2>0 ) { r <- 1; corr <- 1; } real.d and real.b are two vectors, beta1,beta2,theta1,and theta2 are constants. The error occurred like this: Error in if

Hi List, I've got a count matrix like this: [-11.9,-11] (-11,-10.2] (-10.2,-9.3] (-9.3,-8.44] [-11.9,-10.6] 0 0 0 0 (-10.6,-9.3] 0 0 0 0 (-9.3,-8.01] 0 0 0 0 (-8.01,-6.73] 0 0 0 0 I want to

Hi List, Is there anyway to estimate the parameters with integral forms? Here is the function. The idea is to calculate the real mean from the censored mean given normal distribution. I have the function as trun.mean<- function(mu,sigma,t) # t is the threshold { f <- function(x) (1/(sigma*sqrt(2*pi)))*exp(-(x-mu)^2/(2*sigma^2)); pdf.fun <- function(x) x*f(x); mut <-

Hi, I've been running some simulations for a while and the performance of R has been great. However, I've recently changed the code to perform a sort of chi-square goodness-of-fit test. To get the observed values for each cell I've been using table() - specifically I've been using cut2 from Hmisc to divide up the range into a specified number of cells and then using

Hi all, today i tried to enable HA on xcp 1.6 but when i enable it through xencenter I received the following error in xensource.log: Dec 5 17:56:26 xcp1 xapi: [ info|xcp1|490661|Async.pool.enable_ha R:9119ce88ffcb|xapi_ha] Using statefile: /etc/xensource/static-vdis/0/disk Dec 5 17:56:26 xcp1 xapi: [debug|xcp1|490661|Async.pool.enable_ha R:9119ce88ffcb|xapi_ha]

Hi - Just downloaded and installed the memcache-client and cached_model gems and am trying to test it out on a development setup. I added this to my environment/development.rb CACHE = MemCache.new :c_threshold => 10_000, :compression => true, :debug => true, :namespace => ''eztrip'',

Dear R-Help friends, I am unable to get the latest version of R (2.4.1) to compile on my solaris 10 system - has anybody else experienced this problem and are you able to offer me any advice? I appreciate your time, many thanks, Jenny Barnes Here are my CURRENT specifications: platform sparc-sun-solaris2.10 arch sparc os solaris2.10

Dear R helpers, I have following loss data and I need to fit LEFT truncated Log Normal distribution to this data which is Truncated at 1000000. dat = c(1333834,5710254,9987567,7809469,6940935,3473671,1270209,1102523,1124002, 5830159,4302300,3925242,2638409,2324421,7238436,9088709,7439250,4976551,4864319, 8741334,1863770,7098310,4942288,4971829,4986372) library(gamlss.tr) gen.trun(5, LOGNO)

Full_Name: Franz Quehenberger Version: 2.10.1 OS: Windows XP Submission from: (NULL) (145.244.10.3) aggregate is supposed to produce a data.frame that contains a line for each combination of levels of the variables in the by list. The first columns of the result contain these combinations of levels. With two by variables the second by-variable takes always only one value. However, it works fine

>From the manpage: --csum-length=LENGTH By default the primary checksum used in rsync is a very strong 16 byte MD4 checksum. In most cases you will find that a trun- cated version of this checksum is quite efficient, and this will decrease the size of the checksum data sent over the link, mak- ing things faster. You can choose the

Hi I am trying to turn on register_globals, but I am failing. someone trold me that I should change php.ini and I did it. ? - register_globals = Off - register_globals = On I made a php test page html/test.php <?php phpinfo(); ?> and checked it, but I can not make it. output_buffering no value no value output_handler no value no value post_max_size 8M 8M precision 14 14

"Jenny Stadt" <jennystadt at yahoo.ca> asked: > > I am struggling in a bivariate Weibull distribution although I > searched R-Site-Help and found suggestion with Copula. Seems the > maximum likelihood estimate is beyond what I can understand. > > My case is: given two known marginal distribution (both are Weibull), > and the correlation between them. How can I

On Thu, Oct 27, 2016 at 5:39 PM, Maxime Chevalier-Boisvert via llvm-dev <llvm-dev at lists.llvm.org> wrote: > Hello, > > We’d like to keep track of which clang version was used to build our LLVM binaries. We use cmake and ninja with clang to build. What do you people think would be the cleanest way to know which version of clang is used, on a user’s machine, to build those binaries.

Friends I am trying to format a number to a string so 2189.745 goes to "2,189.35" and 309283.929 goes to "309,283.93" I have tried to use formatC(X, big.mark=",",drop0trailing=FALSE, format="f") but it does not get the number of decimals correct. Specifying digits does not work as that is significant digits. I could use a switch statement switching on

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Dan Phillips wrote: > Charles Steinkuehler wrote: >> Dan Phillips wrote: >> >> What we are left with is a requirement in abcde to overcome this and >> >> until then we have the manual method. Have you any thoughts on the best >> >> way to overcome this problem apart from the hacked toc3cue (do you have

Charles Steinkuehler wrote: > Dan Phillips wrote: > >> What we are left with is a requirement in abcde to overcome this and > >> until then we have the manual method. Have you any thoughts on the best > >> way to overcome this problem apart from the hacked toc3cue (do you have > >> a copy of this?) > > I do not have a copy of toc3cue...I made my own

Dear list, Given a list of elements like: aa <- list(one=c("o", "n", "e"), tea=c("t", "e", "a"), thre=c("t", "h", "r", "e")) Is there a function that returns the intersection between all? Both match() and intersect() only deal with two arguments, but sometimes I

> On Oct 27, 2016, at 6:38 PM, Davide Italiano via llvm-dev <llvm-dev at lists.llvm.org> wrote: > > On Thu, Oct 27, 2016 at 5:39 PM, Maxime Chevalier-Boisvert via > llvm-dev <llvm-dev at lists.llvm.org> wrote: >> Hello, >> >> We’d like to keep track of which clang version was used to build our LLVM binaries. We use cmake and ninja with clang to build.

Hey all, I've been seeing this repeatedly and am wondering if anyone has a clue what's causing it.. at least once a day I see two zap fxo channels being bridged, and hanging..now, these two channels should never bridge, but they keep doing it.. any leads on where to look for what's causing it? here's what I see when I 'show channels' Zap/37-1 s@from-pstn:1

Dear all, I need to do some calculation where the code used are below. I get error message when I choose k to be large, say greater than 25. The error message is "Error in integrate(temp, lower = 0, upper = 1, k, x, rho, m) : the integral is probably divergent". Can anyone give some help on resolving this. Thanks. Hannah m <- 100 alpha <- 0.05 rho <- 0.1 F0