similar to: "statistic" term in boot function

Hi Everyone, I have a 2-dim data.matrix(e.g., table1) in which row1 specifies a range of values. row2 - rown specify the number of times I want to replicate each corresponding value in row1. I can do this with the following function: rep(c(table1[1,]),c(table1[X,])) #where X would go from 2 - n. Now, I can do this manually by changing the values of X and save each resulting array/vector in

Buen día, no sé si estoy utilizando bien la lista, es la primera vez. Si lo hago mal me corrigen por favor. Sobre tu comentario Pedro, muchas gracias. Lo qeu entiendo con tu sugerencia de set.seed es qeu de esa forma fijas los resultados, pero no estoy seguro si otra agrupación funcione mejor. Es decir me interesa un método de agrupación que genere la "mejor" agrupación y como los

Hi Folks! I have a question regarding an issue on which I´ve read a few threads but can´t resolve --> Labelling every other tick mark on the x-axis using factors. Here´s an example: db<-data.frame(date=as.factor(c(1:50)),var=rnorm(50)) barplot(db$var,names.arg=levels(db$date),las=2,cex.axis=0.8,cex=0.8) This labels all the tickmarks in my barplot, which is the default behavior. I would

Hi All, I'd like to randomly sample a vector N times, where each successive random sample increases in size. I have realised that the function sample does not take vectors for arguments. For example, x<-rnorm(20,0,1) sample(x,c(1,2,3)) ## will only return one random sample of size 1. The trick seems to be getting past the size argument of the function. I've tried different things

Hi folks, I'm trying to estimate bias-corrected percentile (BCP) confidence intervals on a vector from a simple for loop used for resampling. I am attempting to follow steps in Manly, B. 1998. Randomization, bootstrap and monte carlo methods in biology. 2nd edition., p. 48. PDF of the approach/steps should be available here: https://wyocoopunit.box.com/s/9vm4vgmbx5h7um809bvg6u7wr392v6i9 If

Hi All, I have got matrix 'data' of dimension 22000x600. I want to make 50 independent samples of dimension 22000x300 from the original matrix 'data'. And then want to calculate pearsons CC for each of the obtained 50 matrices. It seems it is possible to do this using 'boot' function from library boot but I am not able to figure out how? I am really stuck. Please help!

Hi all, Are there functions in R that could help me do the following? We have a special type of regression which is called Geometric Mean Regression. We have done some search and found the following: https://stat.ethz.ch/pipermail/r-help/2011-July/285022.html The question is: how to do the statistical inference on GMR results? More specifically, we are looking for the prediction interval:

Hi, I'm currently trying to find/define a relationship between one dependent and several independant variables. The problem is that i cannot use the normal multiple regression/correlation in Spss because the data is not normal distributed. i calculated the spearman roh and Kendalls tau Correlation and also some partial correlations in R. Now i wanna find out the the multiple correlation

Hi everybody, If I've a data frame like this: dataframe a X0 X2 X4 X6 X8 X10 X12 X14 X16 1957 0 0 0 0 0 0 0 0 0 1958 0 0 0 0 0 0 0 0 0 1959 831 0 0 0 0 0 0 0 0 1960 544 282 0 0 0 0 0 0 0 1961 446 365 0 0 0 0 0 0 0 1962 442 473 0 0 0 0 0 0 0 1963 595 468 0 0 0 0 0 0 0 1964

Dear R-Helpers, I am trying to bootstrap the coefficient of variation on a suite of vectors, here I provide an example using one of the vectors in my study. When I ran this script with the vector x <-c(0.625, 0.071428571, 0.133333333, 0.125, 0), it returned CV(boot) [the second one], and stderr(boot) [the second one] without problem. However, when I ran it with the vector in the

Hi Folks, Ran into something I'd really like to do in R simply/elegantly, but my R - coding skills seem surpassed. This is the thing. Imagine the following data: labs<-c("abcdef","abcgg","tgthefdk","tgtijuel","tgtnjmoi","gbnt","dlift") dat<-c(0.5,0.25,1,2,16,0.250,4) dframe<-data.frame(labs,dat) I would like to

Hi, I am trying to run some bootstraps with the boot package. When I run it with 400 replicates it does it ok, but then I need to run the same analysis but with 89, 86, 102 and 106 samples (for four different environments), and then is when I get the error message: > mybootstrap <- boot(Datos, mystat, 2000) Error in t.star[r, ] <- statistic(data, i[r, ], ...) : number of items to

I run the following code: library(rpart) data(kyphosis) fit <- rpart(Kyphosis ~ ., data=kyphosis) plot(fit) text(fit, use.n=TRUE) The text labels represent the count of each class at the leaf node. Unfortunately, the numbers are rounded and in scientific notation rather than the exact number of examples sorted by that node in each class. The plot is supposed to look like

His Friends, I've been using dovecot-1.0.alpha5 for a long time, now i would like to update to v.1.0.rc24 on my Red Hat ES 4 machine, so i need an rpm version that i couldn't find on the web. Also now i can build the rpm package via SPECS, which i don't have. Plis, anybody would help? Thanks in advanced, Alfredo -------------- next part -------------- An HTML attachment was

Hi, I'd like some advice on bootstrapping in R. I have a species x with 20 individuals and a factor containing 0 and 1's (in this case 5 zeros and 15 ones). I want to compare the frequency of the occurrence of 1 with a probability value. This code seems to work to do this in R. attach(test) p <- c(0.5272, (1-0.5272)) sp1_1 <- length(subset(x, x==1)) sp1_0 <- length(subset(x,

Hi all. I wonder if you can help me here. I have centos7 with 1 network interface and on that IFwe have 2 vlan. >From both vlan we'd like to reach the internet independently so basically with 2 different gateways. we tried with all the routes,rules etc but only on one vlan we are able to ping 8.8.8.8 for instance. Any suggestions/ideas? Cheers -- */Alfredo*

I'd to match-merge 2 tables in such a manner that I keep all the rows in table 1, but not the rows that are in both table 1 and 2. Thank you for your help, Alfredo > master <- data.frame(ID=2001:2011) > train <- data.frame(ID=2004:2006) > valid <- ??? in this example table valid should have the following > str(valid) Year: int 2001 2002 2003 2007 2008 2009 2010 2011

I've a matrix like this: 1985 1.38 1.27 1.84 2.10 0.59 3.47 1986 1.05 1.13 1.21 1.54 0.21 2.14 1987 1.33 1.21 1.77 1.44 0.27 2.85 1988 1.86 1.06 2.33 2.14 0.55 1.40 1989 2.10 0.65 2.74 2.43 1.19 1.45 1990 1.55 0.00 1.59 1.94 0.99 2.14 1991 0.92

> -----Mensaje original----- > De: Alfredo Hurtado Pereda > Enviado el: Jueves, 23 de Agosto de 2001 11:38 a.m. > Para: 'wine-users-request@winehq.com'; Alfredo Hurtado Pereda > Asunto: help > > > > Help me Please > > > I have a Problem . I Install Sybase Opclient so i run the > programs isql.exe this program put the Windows this say

Dear All, In a package, I want to use PCA function. The structure I used follow this page: http://www.statmethods.net/advstats/factor.html. fit<-principle(mydata, nfactors=9, rotation=TRUE) or: result<-PCA(mydata) But I don't known why R language in my computer noticed: "not found principle", "not found PCA". I download and installed