search for: nrec

..."datum") # # Not sure if I got the row.names correct, I may need to use as.character (I believe I did use #that in the interactive session). attach(glucose) summary(glucose) ncol <- length(names(glucose)) xrange <- range(1, ncol) yrange <- range(0, max(!is.na(glucose[1:ncol]))) nrecs <- nrow(glucose) colors <- rainbow(nrecs) linetype <- c(1:nrecs) plot(xrange, yrange, type="n", xlab="Measurement moment",  ylab="Glucose (mmol/liter)") for (i in 1: nrecs) {   daily_values <- glucose[i,]   lines(daily_values, type="b", lwd=1.5,...

...ots_1.jpg",height=600,width=600) pdf("Z:\\My Documents\\PROJECTS\\Access Policy\\AccessDAPlots_1.pdf") #png("Z:\\My Documents\\PROJECTS\\Access Policy\\AccessDAPlots_1.png") #tiff("Z:\\My Documents\\PROJECTS\\Access Policy\\AccessDAPlots_1.tiff") while (irec <= nrec) { con_new<-consultants[irec] spec_new<-specialty[irec] spec_cons_new<-spec_cons[irec] if (spec_cons_new!=spec_cons_old || irec==nrec ) { strip_name[con_count]<-paste(spec_old,'\n',con_old) if (spec_new!=spec_old || irec==nrec) { spec_old&...

...ke this string up with the file name I get from file.choose(): filename = file.choose() #get the number of lines in the file # first make a command string with the filename in it cmd = paste('findstr /R /N "^" ', filename, ' | find /C ":"', sep=""') nrec = as.numeric(shell(cmd)) But R puts in escape characters for the ": > cmd [1] "findstr /R /N \"^\" D:\\my_dir\\my_file | find /C \":\"" and shell(cmd) does not work properly. And > nrec=shell('findstr /R /N "^"D:\\my_dir\\my_file | find /C &q...

...he memory leak is in readChar header <-paste("This is a very long and boring text header which appears", "at the beginning of each chunk of an even longer and duller", "file which contains mostly binary data") nchars <- nchar(header) nrecs <-2000 ntries <- 100 ## Write a test file stream <- file("/tmp/gunge") open(stream,open="wb") for(irec in 1:nrecs){ ## This writes null-ended strings unless you use the eos=NULL option writeChar(header,stream,eos=NULL) } close(stream) ## Read in the file ntries t...

...vs <- as.matrix(covs) }else{ n.covs <- 0 } } ordered.t <- t(apply(cbind(onset,time.dep),1,sort,na.last=TRUE)) tot.time.dep <- apply(ordered.t,1,function(x) sum(!is.na(x))) ordered.t <- cbind(rep(0, nrow(ordered.t)), ordered.t) npars <- 4+n.time.dep+n.covs nrecs <- sum(tot.time.dep) new.x <- as.data.frame(matrix(nr=nrecs, nc=npars)) names(new.x) <- c("start", "stop", "event", names(time.dep),names(covs),"episode") this.rec<-0 for(i in 1:length(onset)) { for(j in 1:tot.time.dep[i]) { thi...

Hi List! I?m using regsimq() from the ?lmomRFA?-package to calculate error bounds for diverse distributions. For example: regsimq(gumfit$qfunc, nrec = lmom.data$n, f = lcdfgum, boundprob = c(0.025, 0.975)) Several times I got this error massage: Fehler in quantile.default(ou, probs = boundprob, type = 6) : missing values and NaN's not allowed if 'na.rm' is FALSE So my question is, can I change 'na.rm' = FALSE into '...

....beg -360)} #since sp prefers -180 to 180 encoding lons <- seq(lon.beg, lon.end, length=nx) lons[lons > 180] <- lons[lons > 180]-360 #since sp prefers -180 to 180 encoding #read flat file into R (header first) dims <- scan("/tmp/wgrib.output", nlines=1, quiet=T) nrec <- dims[1]*dims[2] indata <- paste("indata.", i, j, mon2, sep="") indata <- scan("/tmp/wgrib.output", 'numeric', nlines=nrec, skip=1, quiet=T) indata <- matrix(as.numeric(indata), nrow=dims[1], ncol=dims[2], byrow=F) #to get precip into mm/m...

...>to 180 encoding > lons <- seq(lon.beg, lon.end, length=nx) > lons[lons > 180] <- lons[lons > 180]-360 #since sp >prefers -180 >to 180 encoding > > #read flat file into R (header first) > dims <- scan("/tmp/wgrib.output", nlines=1, quiet=T) > nrec <- dims[1]*dims[2] > indata <- paste("indata.", i, j, mon2, sep="") > indata <- scan("/tmp/wgrib.output", 'numeric', nlines=nrec, >skip=1, >quiet=T) > indata <- matrix(as.numeric(indata), nrow=dims[1], ncol=dims[2], > >byrow=F)...

Hello, I am using the lmomco package (lmom.ub and pargev) to compute the GEV parameters (location, scale, and shape), which are used to estimate return values. I was wondering how/if I can calculate upper and lower confidence (CI_u, CI_l) intervals for each return frequency using the GEV parameters to fill-in the table below? Xi (location) = 35.396 Alpha (scale) = 1.726 Kappa (shape) =

...nths ago, but only recently discovered this service problem. The message we get when this happens is the predictable "trust relationship failed" error while trying to access a non-associated service. jonathan -- Jonathan Woytek w: 412-681-3463 woytek+@cmu.edu NREC Computing Manager c: 412-401-1627 KB3HOZ PGP Key available upon request

Hi, Searched for similar errors and nothing came up. Anybody know what this is? When using aggregate stddev, and then trying to print it at END, using printa, I get this error: dtrace: processing aborted: Invalid return value from callback avg works fine. Tried on two recent builds of solaris (build 89 and build 94), on two difference sparc systems. Here''s a sample script:

...number %lu", ino); + + free(buf); + free(leaf); + + return ip; + +failed: + free(ip); + free(buf); + free(leaf); + + return ip; +} + +static xfs_fsblock_t +select_child(xfs_dfiloff_t off, + xfs_bmbt_key_t *kp, + xfs_bmbt_ptr_t *pp, + int nrecs) +{ + int i; + + for (i = 0; i < nrecs; i++) { + if (be64_to_cpu(kp[i].br_startoff) == off) + return be64_to_cpu(pp[i]); + if (be64_to_cpu(kp[i].br_startoff) > off) { + if (i == 0) + return be64_to_cpu(pp[i]); + else +...