similar to: Jackknife and rpart

Dear R users, I?m student of Master in Statistic and Data analysis, in New University of Lisbon. And now i?m writting my dissertation in variance estimation.So i?m using Survey Package to compute the principal estimators and theirs variances. My data is from Incoming and Expendire Survey. This is stratified Multi-stage Survey care out by National Statistic Institute of Mozambique. My domain of

I have a matrix with NAs and want to extract the longest column. > is.matrix(foo) [1] TRUE > dim(foo) [1] 2000 75 > GetLength <- function(x) {length(na.omit(x))} > junk <- apply(foo, 2, GetLength) > junk [1] 1004 512 432 523 691 396 607 838 [9] 730 389 388 445 609 333 637 1024 [17] 1163 823 718 466 799 459 701 833 [25] 456 549 376 728 539 384

Greetings. I'm trying to determine whether to use rpart or randomForest for a classification tree. Has anybody tested efficacy formally? I've run both and the confusion matrix for rf beats rpart. I've looking at the rf help page and am unable to figure out how to extract the tree. But more than that I'm looking for a more comprehensive user's guide for randomForest including

Hi, I have the following jackknife code which is much slower than my colleagues C code. Yet I like R very much and wonder how R experts would optimize this. I think that the for (i in 1:N_B) part is bad because Rprof() said sum() is called very often but I have no idea how to optimize it. #O <- read.table("foo.dat")$V1 O <- runif(100000); k=100 # size of block to delete

Hi R-experts, I want to use the jackknife function from the bootstrap package onto a linear model. I can't figure out how to do that. The manual says the following: # To jackknife functions of more complex data structures, # write theta so that its argument x # is the set of observation numbers # and simply pass as data to jackknife the vector 1,2,..n. # For example, to jackknife #

Hi dear all, Can aynone help me about delete-d jackknife usually normal jackknife code for my data is: n <- nrow(data) y <- data$y z <- data$z theta.hat <- mean(y) / mean(z) print (theta.hat) theta.jack <- numeric(n) for (i in 1:n) theta.jack[i] <- mean(y[-i]) / mean(z[-i]) bias <- (n - 1) * (mean(theta.jack) - theta.hat) print(bias) but how i can apply delete-d jackknife

Dear all I have used the lda and predict functions to classify a set of objects of unknown origin. I would like to use a jackknife reclassification to assess the degree to which the outcomes deviate from that expected by chance. However, I can't find any function that allows me to do this. Any suggestions of how to generate the jackknife reclassification to assess classification accuracy?

Hi dear all, Can someone help me about detection of outliers using jackknife after bootstrap algorithm? -- View this message in context: http://r.789695.n4.nabble.com/jackknife-after-bootstrap-tp3041634p3041634.html Sent from the R help mailing list archive at Nabble.com. [[alternative HTML version deleted]]

Folks, Thanks for the help with the hier.part analysis. All the problems stemmed from an import problem which was solved with file.chose(). Now that I have the variables that I'd like to use I need to run some GLM models. I think I have that part under control but I'd like to use a jackknife approach to model validation (I was using a hold out sample but this seems to have fallen out

Is there a routine in the boot package to get the jackknife-after- bootstrap influence values? That is, the influence values of a jackknife of the bootstrap estimates? I can see how one would go about it from the jack.after.boot code, but that routine only makes pretty pictures. It wouldn't be hard to write, but I find it hard to believe this isn't part of the package already. Thanks

Hi All, I'm not able to figure out how to perform a Jackknife test for a 2-sample dispersion test in R. Is there a built-in function to perform this or do we have to take a step by step approach to calculate the test statistic? Any help would be awesome. Thanks! -- View this message in context: http://r.789695.n4.nabble.com/Jackknife-for-a-2-sample-dispersion-test-tp4444274p4444274.html

Greetings, I?ve fit the following mixed effects model using the NLME package: hd.impute.lme <- lme(I(log(HEIGHT_M - 1.37)) ~ SPECIES + SPECIES:I(1/(DBH_CM + 2.54)), random = ~ I(1/(DBH_CM + 2.54)) | PLOTID, data = trees, na.action = na.exclude) I would now like to extract a jackknife estimate of model error. I tried the following code, however, the estimate produced seems too

Dear R friends I´m interested into apply a Jackknife analysis to in order to quantify the uncertainty of my coefficients estimated by the logistic regression. I´m using a glm(family=’binomial’) because my independent variable is in 0 - 1 format. My dataset has 76000 obs, and I´m using 7 independent variables plus an offset. The idea involves to split the data in let’s say 5 random subsets and

I'm using the lme function in nmle to estimate the variance components of a fully nested two-level model: Y_ijk = mu + a_i + b_j(i) + e_k(j(i)) lme computes estimates of the variances for a, b, and e, call them v_a, v_b, and v_e, and I can use the intervals function to get confidence intervals. My understanding is that these intervals are probably not that robust plus I need intervals on the

# I need to generate some data. I'm modeling some time series that follow a # negative exponential decay (mostly). I have 20 samples that can easily be fit with cubic splines. # What I want to do is generate many thousands of similar samples using the parameters from the splines # For instance one data sample looks not unlike this: foo.curve <- 1 * exp(-0.01 * 1:500) + 0.5

Hello R Community, I've been run the following codes. However, I've been getting an unusual segfault that I'm unable to trace its origin. Please give me a light to decipher the "caught segfault" Thanks for you attention. Bernardo. > options(STERM='iESS', editor='emacsclient') > rm(list = ls()) > > source("fgenIGLD.R") #RNG

Dear List I am using R with mcgv package to model spatial variation in density estimates of dorcas gazelle in Sinai. I have 59 points of data and 4 explanatory variables(distance from mountain edge, camel presence, Latitude & Longitude). I want to test the model fir via bootstraping. I have used the jacknife bootstraping but it have the limitation of allowing only 58 trials. I tried to use the

Afternoon R help, I want to run Rasch/IRT analyses using the ltm package, however, I am using large scale survey data which requires weighting for accurate results. I attempted to create a weighted object to insert into the formulae of the ltm packages, however, the survey data only includes 30 replicate weights and a sampling weight. The svrepdesign requires additional information such as

Hello All! I''m novice at ROR. Would you like to explain what mean such constuction in Ruby. 1: class DemoController < ActionController::Base 2: scaffold :item 3: end Is 2th line calling of method ''scaffold''? If it is so, when it called? For all instances of class, or only for class? I watch sources of ActionController, but don''t understand completely.

Hello, please excuse the simplicity of this question as I am not very good with stats. I am taking a class, using R which I am learning at the same time, and the questions asks us to "Estimate the average abundance using Poisson regression with a log link". I can estimate the abundance from "x", but I can seem to figure out how to get the average abundance in this method. Any