r - Issue adding inputs/calculations in FOR Loop -



r - Issue adding inputs/calculations in FOR Loop -

as title states, have problem adding calculated values empty vector, , in case of loop 1000 iterations, insert same value 1000 times, instead of different values each of 100 iterations. sifted through code, line line, setting , j values forcefulness stop iterations. see below code:

# population size n <- 100000 # iterations per sample size n.sim <- 1000 # various sample sizes samples.sim <- seq(from=100, to=5000, by=50) cap.recap.partdeux <- function(n, n.sim, samples.sim){ n <- na # sample size bias.of.est <- na # bias of estimator sd.of.est <- na # standard deviation of estimators sample l.sim.chap <- na # temporary list of estimators utilize later # loop: various sample sizes for(i in 1:length(samples.sim)){ # loop: iterations per sample size (j in 1:n.sim){ <- 1:1 j <- 1:2 # grab 1 sim.one <- sample(n, samples.sim, replace=t, prob=null) sim.one <- as.numeric(sim.one) # convert numeric # grab 2 sim.two <- sample(n, samples.sim, replace=t, prob=null) sim.two <- as.numeric(sim.two) # convert numeric # find mutual elements sim.m.two <- intersect(sim.one, sim.two) # amount of mutual elements sim.l.m.two <- length(sim.m.two) # calculate chapman estimator sim.chap <- ((samples.sim[i]+1)*(samples.sim[i]+1)/(sim.l.m.two+1))-1 l.sim.chap[j] <- list(sim.chap) } # end loop: iterations per sample size # calculate bias of estimator each sample sum.b.est <- sum(unlist(l.sim.chap), na.rm=t) bias.est <- (sum.b.est/n.sim)-n bias.of.est[i] <- bias.est # calculate standard deviation of estimator each sample sd.est <- sd(unlist(l.sim.chap), na.rm = true) sd.of.est[i] <- sd.est # sample size n[i] <- samples.sim[i] } # end loop: various sample sizes # homecoming 3 columns , create info frame 3 <- (data.frame(n, bias.of.est, sd.of.est)) # list of info frame true population output <- (list(three, "pop"=n)) return(output) } # end function output <- cap.recap.partdeux(n=100000, n.sim=1000, samples.sim) test <- data.frame(output)

basically, on line l.sim.chap[i] <- list(sim.chap), vector l.sim.chap loaded repeated values first sim.chap iteration, not sim.chap values each , subsequent iteration.

removed <-1 , j <- 1. code works now. strange.

# population size n <- 100000 # iterations per sample size n.sim <- 1000 # various sample sizes samples.sim <- seq(from=100, to=5000, by=50) cap.recap.partdeux <- function(n, n.sim, samples.sim){ n <- na # sample size bias.of.est <- na # bias of estimator sd.of.est <- na # standard deviation of estimators sample l.sim.chap <- na # temporary list of estimators utilize later # loop: various sample sizes for(i in 1:length(samples.sim)){ # loop: iterations per sample size (j in 1:n.sim){ # grab 1 sim.one <- sample(n, samples.sim, replace=t, prob=null) sim.one <- as.numeric(sim.one) # convert numeric # grab 2 sim.two <- sample(n, samples.sim, replace=t, prob=null) sim.two <- as.numeric(sim.two) # convert numeric # find mutual elements sim.m.two <- intersect(sim.one, sim.two) # amount of mutual elements sim.l.m.two <- length(sim.m.two) # calculate chapman estimator sim.chap <- ((samples.sim[i]+1)*(samples.sim[i]+1)/(sim.l.m.two+1))-1 l.sim.chap[j] <- list(sim.chap) } # end loop: iterations per sample size # calculate bias of estimator each sample sum.b.est <- sum(unlist(l.sim.chap), na.rm=t) bias.est <- (sum.b.est/n.sim)-n bias.of.est[i] <- bias.est # calculate standard deviation of estimator each sample sd.est <- sd(unlist(l.sim.chap), na.rm = true) sd.of.est[i] <- sd.est # sample size n[i] <- samples.sim[i] } # end loop: various sample sizes # homecoming 3 columns , create info frame 3 <- (data.frame(n, bias.of.est, sd.of.est)) # list of info frame true population output <- (list(three, "pop"=n)) return(output) } # end function output <- cap.recap.partdeux(n=100000, n.sim=1000, samples.sim) test <- data.frame(output)

r for-loop vector

Comments

Post a Comment

Popular posts from this blog

java Multi query from Mysql using netbeans -

java Multi query from Mysql using netbeans - i have java code retrieve info mysql appropriate code in netbeans.the code below. want multi query such illustration finding total number of products , on.could help me please. in advance import java.sql.*; public class javamysqlselectexample { public static void main(string[] args) { seek { string mydriver = "com.mysql.jdbc.driver"; string myurl = "jdbc:mysql://localhost/sample"; class.forname(mydriver); connection conn = drivermanager.getconnection(myurl, "root", "mypasscode1"); string query = "select * products"; statement st = conn.createstatement(); resultset rs = st.executequery(query); while (rs.next()==true) { int code = rs.getint("code"); string productname = rs.getstring("product_name"); string producttype = rs.getstring("product_type"); system.out.format(" %d, %s, %s\n", code, pr
Read more

c# - DotNetZip fails with "stream does not support seek operations" -

c# - DotNetZip fails with "stream does not support seek operations" - i using dotnetzip in c# unzip stream follows: public static void unzipfromstream(stream stream, string outdir) { //omit seek grab block using (zipfile zip = zipfile.read(stream)){ foreach (zipentry e in zip){ e.extract(outdir, extractexistingfileaction.overwritesilently); } } } stream obtained using webclient client = new webclient(); stream fs = client.openread(url); however, got next exception exception during extracting zip stream system.notsupportedexception: stream not back upwards seek operations. @ system.net.connectstream.get_position() @ ionic.zip.zipfile.read(stream zipstream, textwriter statusmessagewriter, encoding encoding, eventhandler`1 readprogress) on server side(asp.net mvc 4), returning filepathresult or filestreamresult both caused exception. should obtain stream differently on client side? or how create server homec
Read more

c++ - StartServiceCtrlDispatcher don't can access 1063 error -

c++ - StartServiceCtrlDispatcher don't can access 1063 error - i write on c++ service , have error, , don't can fixed it. my service main function int main (int argc, tchar *argv[]) {dword f; for(int i=0;i<113;i++) f = getlasterror(); servicepath = lptstr(argv[0]); outputdebugstring(_t("my sample service: main: entry")); //installservice(); service_table_entry servicetable[] = { {service_name, (lpservice_main_function) servicemain}, {null, null} }; // startservice(); if(!startservicectrldispatcher(servicetable)) { f = getlasterror(); //addlogmessage("error: startservicectrldispatcher"); } else if( memcmp(argv[argc-1],"install",7)) { installservice(); } else if( memcmp(argv[argc-1],"remove",6)) { removeservice(); } else if( memcmp(argv[argc-1],"start",5)) { startservice(); } else if( memcmp(argv[argc-1],"stop",4))
Read more