search for: nrecs

..."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, l...

Hi.....I am relatively new to R but was quite pleased with myself at having generated a series of lattice plots as PDFs. I was very surprised when plotting these out as jpegs (or png or tiff) that the strip title information above each lattice plot vanished. The pdf was fine. Has anybody any ideas? I can't add an image as the information is sensitive. Many Thanks Steve Creamer Here is the

On a Windows platform I am trying to count the number of lines in a file. In a DOS window, the following works: C:\Users\jar>findstr /R /N "^" D:\my_dir\my_file | find /C ":" 5317 (it works with double \\ also) But in R, I need to make 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

...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 ti...

...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]) { this...

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

Dear R-help community, My main goal of this message is to find a way of skipping a file of a month/year in a loop that does not exist (and making it's output into an data.out array would be NA) and moving onto the next year/month in the loop to carry on filling data.out with real precipitation data. The situation so far: I downloaded 50 years worth of GRIB data files from the NCEP data

Hi Jenny If you want a general solution I understand. However I just downloaded the file fine (as far as I can tell) so you are welcome to have a copy. I can email it to you if you want. I do not think your test for NA is valid. i.e if(test != "NA"){ } I think you should use if(is.na(test)){ } Or something similar. J --- John Seers Institute of Food Research Norwich

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) =

File this under the "this used to work" section. Particulars: Version: 3.0.13 Build source: Source RPM OS/Distribution: RedHat Enterprise Linux AS 3 Hardware: dual-Xeon, 4GB RAM, SCSI OS disks Security configuration: NT4 domain, NT4 PDC and BDC This system is being run under RedHat Cluster Services. We have three samba instances running on one physical machine, providing windows

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 +...