similar to: bootstrap generalization error

Dear helpers in this forum, This is a clarified version of my previous questions in this forum. I really need your generous help on this issue. > Suppose I have the following data set: > > id x y > 023 1 2 > 023 2 5 > 023 4 6 > 023 5 7 > 412 2 5 > 412 3 4 > 412 4 6 > 412 7 9 > 220 5 7 > 220 4 8 > 220 9 8 > ...... > Now I want to compute the

Well, I could argue that it's not *completely* OT since my question is motivated by an enquiry that I received in respect of a CRAN package "Iso" that I wrote and maintain. The question is this: Given observations y_1, ..., y_n, what is the solution to the problem: minimise \sum_{i=1}^n (y_i - y_i^*)^2 with respect to y_1^*, ..., y_n^* subject to the "isotonic"

Dear R Help, I am trying to fit a nonlinear model for a mean function $\mu(Data_i, \beta)$ for a fixed covariance matrix where $\beta$ and $\mu$ are low- dimensional. More specifically, for fixed variance-covariance matrices $\Sigma_{z=0}$ and $\Sigma_{z=1}$ (according to a binary covariate $Z $), I am trying to minimize: $\sum_{i=1^n} (Y_i-\mu_(Data_i,\beta))' \Sigma_{z=z_i}^{-1} (Y_i-

Hello, I have a bootstrap implementation loop that I would like to replace by a faster "batch operator". The loop is as follows: for (b in 1:B) { indices <- sample(1:N, size=N, replace=TRUE) # sample n elements with replacement theta_star[b,] = statistic(data,indices) # execute statistic callback function } I would like to replace this by something like (does

Spencer Thanks for your thoughts on this. I did a bit of work and did end up with a method (more a trick), but it did work. I am certain there are better ways to do this, but here is how I resolved the issue. The integral I need to evaluate is \begin{equation} \frac{\int_c^{\infty} p(x|\theta)f(\theta)d\theta} {\int_{-\infty}^{\infty} p(x|\theta)f(\theta)d\theta} \end{equation} Where

Hello, I have the following function that receives a "function pointer" formal parameter name "fnc": loocv <- function(data, fnc) { n <- length(data.x) score <- 0 for (i in 1:n) { x_i <- data.x[-i] y_i <- data.y[-i] yhat <- fnc(x=x_i,y=y_i) score <- score + (y_i - yhat)^2 } score <- score/n

Author: duvander Date: Sun May 31 17:34:27 2009 New Revision: 12 Modified: trunk/source/mxn.cloudmade.core.js Log: Implemented most Cloudmade functions. Still requires image, tile and KML/GeoRSS overlays Modified: trunk/source/mxn.cloudmade.core.js ============================================================================== --- trunk/source/mxn.cloudmade.core.js (original) +++

Revision: 30 Author: freyfogle Date: Sun Sep 20 06:26:44 2009 Log: jslint now passing for google and googlev3 http://code.google.com/p/mapstraction/source/detail?r=30 Modified: /trunk/build.xml /trunk/source/mxn.google.core.js /trunk/source/mxn.googlev3.core.js ======================================= --- /trunk/build.xml Sun Sep 20 06:13:41 2009 +++ /trunk/build.xml Sun Sep 20 06:26:44

...as promised. This describes the mdct used in my d.m.l patch. I think it is the same as the Lee fast-dct. I typed it in a kind of pseudo-TeX, 'cause the ascii art would kill me. Hope you can read TeX source; if not, ask someone who can to make a .ps/.gif/.whatever of the TeX output, and put it on a webpage or something. I'm to lazy to do it (and besides, I don't have access to TeX,

Author: robert.moran Date: Sun May 31 11:13:23 2009 New Revision: 11 Modified: trunk/source/mxn.core.js trunk/source/mxn.js trunk/source/mxn.yahoo.core.js trunk/tests/index.htm Log: Added functionality to automatically load scripts needed. Allowed mapstraction default constructor to have provider omitted (will use first loaded provider). Minor change to yahoo implementation -

Revision: 65 Author: ajturner Date: Sun Dec 13 11:10:04 2009 Log: Use the GeoCommons addOverlay method http://code.google.com/p/mapstraction/source/detail?r=65 Modified: /trunk/source/mxn.geocommons.core.js /trunk/tests/index.htm ======================================= --- /trunk/source/mxn.geocommons.core.js Sun Dec 13 11:10:00 2009 +++ /trunk/source/mxn.geocommons.core.js Sun Dec 13

Hi Everybody I want to print a latex table containing math characters in the column heading These are the formulae I want to use as column headings. It prints OK from TeX $\sum_{i}\sum_{j}C_{P,i,j,y}\times\mathit{FC}_{i}$, $XU_{alt,y}$, $n$, $\bar{C}_{P,y}$ My plan was to create a character vector with these and later rbind the values to them. When I create the vector like:

Dear R users, I would like to optimize a linear approximation of a quadratic function using lpSolve. My code runs without any error or warning message but the constraints that I set don't seem to work properly. Nevertheless, I am certain that my code is somewhere wrong. I would like to solve the following problem: max 2x-x^2+y subject to 2x^2 + 3y^2 <= 6 2>= x,y >= 0 I would

I am trying to use the marker method addlistener(). A simplified version of what I''m trying to do: var map; var marker; map = new mxn.Mapstraction(''myMap'',''google''); map.setCenterAndZoom(new mxn.LatLonPoint(45.559242,-122.636467), 15); marker = new mxn.Marker(new mxn.LatLonPoint(45.559242,

Revision: 53 Author: freyfogle Date: Wed Nov 4 08:34:55 2009 Log: applied patch from John McKerrell http://code.google.com/p/mapstraction/source/detail?r=53 Modified: /trunk/source/mxn.cloudmade.core.js ======================================= --- /trunk/source/mxn.cloudmade.core.js Sun Sep 20 06:30:15 2009 +++ /trunk/source/mxn.cloudmade.core.js Wed Nov 4 08:34:55 2009 @@ -7,6 +7,21 @@

I am trying to test out several mgcv::gam models in a scalar-on-function regression analysis. The following is the 'hierarchy' of models I would like to test: (1) Y_i = a + integral[ X_i(t)*Beta(t) dt ] (2) Y_i = a + integral[ F{X_i(t)}*Beta(t) dt ] (3) Y_i = a + integral[ F{X_i(t),t} dt ] equivalents for discrete data might be: 1) Y_i = a + sum_t[ L_t * X_it * Beta_t ] (2) Y_i

Hi the list, According to what I know, the Canberra distance between X et Y is : sum[ (|x_i - y_i|) / (|x_i|+|y_i|) ] (with | | denoting the function 'absolute value') In the source code of the canberra distance in the file distance.c, we find : sum = fabs(x[i1] + x[i2]); diff = fabs(x[i1] - x[i2]); dev = diff/sum; which correspond to the formula : sum[ (|x_i - y_i|) /

for learning purposes and also to help someone, i used roger peng's document to get the mle's of the gamma where the gamma is defined as f(y_i) = (1/gammafunction(shape)) * (scale^shape) * (y_i^(shape-1)) * exp(-scale*y_i) ( i'm defining the scale as lambda rather than 1/lambda. various books define it differently ). i found the likelihood to be n*shape*log(scale) +

Dear all, Given a LME model (following the notation of Pinheiro and Bates 2000) y_i = X_i*beta + Z_i*b_i + e_i, is it possible to extract the variance-covariance matrix for the estimated beta_i hat and b_i hat from the lme fitted object? The reason for needing this is because I want to have interval prediction on the predicted values (at level = 0:1). The "predict.lme" seems to

Revision: 83 Author: dezfowler Date: Fri May 21 11:57:43 2010 Log: Making the geocoder pass JSLint http://code.google.com/p/mapstraction/source/detail?r=83 Modified: /trunk/source/mxn.google.geocoder.js ======================================= --- /trunk/source/mxn.google.geocoder.js Wed May 19 16:33:35 2010 +++ /trunk/source/mxn.google.geocoder.js Fri May 21 11:57:43 2010 @@ -28,7 +28,7 @@