similar to: distribution overlap - how to quantify?

Hi All, I have a problem to create a variable that is a function of an integral of another function. The problem is the following: I have a variable called cip. I have to create another variable called bip that is a function of the former variable cip and also the cumulative distribution function of the variable cip. Consider the equation below:

I think there are some missing words in "R Internals"'s description of the P_ macro. It currently has "A macro as a wrapper for ngettext", which I think ought to be something like "The macro P_ may be used as a wrapper for ngettext". The following patch also makes the 2 alternate definitions of P_ have the same argument names, StringS and StringP. Expanding the

Hi all, I am wondering if anyone has implemented (or at least tried to) an automatic reparametrization in order to satisfy "trivial" constraints (in the sense of Dennis & Schnabel, 1983) in optimization problems. To be perhaps clearer let us consider a simple bi-exponential model for some recorded signal (sorry for the LaTex notations I hope they aren't too confusing): $s(t) = A

I need to refresh my memory on Probability Theory, especially on conditional probability. In particular, I want to solve the following two problems. Can somebody point me some good books on Probability Theory? Thank you! 1. Z=X+Y, where X and Y are independent random variables and their distributions are known. Now, I want to compute E(X | Z = z). 2.Suppose that I have $I \times J$ random number

Hi all, I wanto to run a plot about the levels of a variable parting on an ols regression. The regression in done on the rate of return of the variable. Imagine R_{t}=a+b*R_{t-1} So If P, the "estimated" price would be P_{t}=P_{t-1}*R_{t} Imagine that I obtain lm fitted values and the original R_{t}, R_{t-1} are [1000,1] dimmension... How it is posible to obtain also 1000,1 fitted

The Helmholtz Alliance of Earth System Dynamics aims at the development of innovative remote sensing products based on the newest satellite missions. The focus here is on quantifying dynamic processes within the bio-, hydro-, geo- and cryosphere to monitor and understand environmental change on a global scale. The Alliance will be the nucleus within an innovative network of more than 60

Dear All, I am documenting a R package. That means writing the *.Rd files inside the \man folder of the package structure I was wondering how to write the symbol for an integral function in a formula. Similar to this one in LaTeX: \int_{0}^{10} \Omega(t)dt I already tried \deqn{\int_{0}^{10} \Omega(t)dt} but it does not work. Any idea? Which math symbols does R-help recognise? Regards,

Hi, Can Someone inform me how to generate RV's using the below CDF, by inverse technique. Thanks for your help and time. My CDF is as follows \[ F(x)=0 \ \text{if} \ x < 0\]\[ F(x)=\{\frac{x-x_i}{x_{i+1}-x_{i}}*(p_{i+1}-p_{i})\}+p_{i}\ \forall \ x_{i}\leq x < x_{i+1} \] \[ F(x)=1 \ \text{if} \ x > x_{i+1} \] Regards Murthy

Ok, I think my first mail didn't get through so I'll try this again. Please disregard. I'm trying to use subscripts on plot text. The help files say this can be done and give several example. I'm trying to do: P_{2} (for example) and so I tried: text(5,2,expression(P_{2})) which doesn't work. There are no subscript examples. I recall someone telling me that latex-like

How can I get the log odds associated with the levels in strata() within a clogit() model? I'm running R-1.9.0 on a Linux platform. I am using clogit() to run a Rasch model in Item Response Theory in psychometrics. Symbolically, the model is: logit(p_{j,k}) = \log({\Pr(p_{j,k}) \over \Pr(1-p_{j,k})}) = \theta_j - \alpha_k, That is, the log odds of answering an test item correctly is

Dear R-Helpers, this is my first post ever to a mailing list, so please feel free to point out any missunderstandings on my side regarding the conventions of this mailing list. My problem: Assuming the following character vector is given: names <- c("filia Maria", "vidua Joh Dirck Kleve (oo 02.02.1732)", "Bernardus Engelb Franciscus Linde j.u.Doktor referendarius

Hello all, I'm hoping to find a way to evaluate the following sort of integral in R. \int_a^b \int_{g(y)}^Inf f(x,y) dx dy. The integral has no closed form and so must be evaluated numerically. The "adapt" package provides for multidimensional integration but does not appear to allow the limits of integration to be a function. I need to evaluate a number of integrals of this

I tried to manually assign each of 3 arrays a unique TBAA node. But it does not seem to help: alias analysis still considers arrays as may-alias, which most likely prevents the desired optimization. Below is the sample code with TBAA metadata inserted. Could you please suggest what might be wrong with it? Many thanks, - D. marcusmae at M17xR4:~/forge/llvm$ opt -time-passes -enable-tbaa -tbaa

Hi, I am compiling R-0.49 on my SGI Indigo 2 R8000 machine and it stopped in the middle. Could anyone help me to fix it? Many thanks! -- Xiaoming Cai School of Geography .&______~*@*~______&. University of Birmingham "w/%%%%%%%%%%%%%%%%%%%\w" Edgbaston, Birmingham

I know this has been dealt with before on this list, but the previous messages lacked detail, and I haven't figured it out yet. The model is: \x_{ij} = \mu + \alpha_i + \beta_j \alpha is a random effect (subjects), and \beta is a fixed effect (condition). I have a link function: p_{ij} = .5 + .5( 1 / (1 + exp{ -x_{ij} } ) ) Which is simply a logistic transformed to be between .5 and 1.

Hello, I'm trying to fit a special kind of proportional odds model from: Whitehead et al. (2001). Meta-analysis of ordinal outcome using individual patient data. Statistics in medicine 20: 2243-2260. (model 2) The data are as follows: library(nlme) library(geepack) library(Design) library(MASS) options(contrasts=c("contr.SAS","contr.poly")) counts <-

Hi, In the code below, MyNamespace::MyClass_ptr will be named "struct.boost::MyClass_ptr" by llvm-g++, and not "struct.MyNamespace::MyClass_ptr" as I expected. I observed this with the real boost::shared_ptr but used the code below (also attached as structnametest1.cc) on the demo page (http://llvm.org/demo/index.cgi) to reproduce the behavior. When I extended my test code

On Sun, Aug 23, 2009 at 8:28 PM, Chris Lattner<clattner at apple.com> wrote: > > On Aug 23, 2009, at 5:50 PM, OvermindDL1 wrote: > >> On Sun, Aug 23, 2009 at 6:32 PM, Daniel Dunbar<daniel at zuster.org> wrote: >>> >>> This is too heavy, and we don't need the extra features, and regexec >>> is well tested and much more standard. Unless there

Hi, I want to calculate expected likelihood over a 2D discrete distribution, something like \sum_k \sum_l p_{kl} L(v_k, v_l) It is trivial to write a cycle like for(k in 1:K) for(l in 1:L) sum <- sum + p[k,l]*L(v[k], v[l]) But is there a more clever way for R? Best wishes, Ott

Hi I'm starting to use R in windows. I need unzip some files. Each of them are in different folders and subfolders. Names of folders have the same beginning of the name. [example:( ~ / M_ /P_ (...)/ (...)R(...)/*.zip)] I do not know the expression (character) necessary to call this direction and to unzip files are in each subfolder. Thank you. -- View this message in context: