Very nice indeed. Thank you gentlemen. Steven
Michael Friendly 於 2019/7/24 上午 01:23 寫道:
Nice to see William Dunlap take the trouble to mimic the classic Fortran behavior of printing **** for numbers that don't fit in the given width :)-Michael On 7/22/19 6:33 p.m., William Dunlap via R-help wrote:The following mimics Fortran printing with format <perLine>F<fWidth>.<fPrecision>. print1 <- function (x, perLine = 10, fWidth = 8, fPrecision = 2, fortranStars = TRUE) { format <- paste0("%", fWidth, ".", fPrecision, "f") oldWidth <- getOption("width") on.exit(options(width = oldWidth)) options(width = perLine * fWidth) fx <- sprintf(format, x) if (fortranStars) { fx[nchar(fx) > fWidth] <- strrep("*", fWidth) } cat(fx, sep = "", fill = TRUE) invisible(x) } Compareprint1((-1.7)^(1:24))-1.70 2.89 -4.91 8.35 -14.20 24.14 -41.03 69.76 -118.59 201.60 -342.72 582.62 -990.46 1683.78-2862.42 4866.12-8272.4014063.08********40642.31 ******************************** with the output from the Fortran % cat a.f double precision x(24); integer i do 10 i=1,24 x(i) = (-1.7d0) ** i 10 continue write(6, "(10f8.2)") x stop end % gfortran a.f % ./a.out -1.70 2.89 -4.91 8.35 -14.20 24.14 -41.03 69.76 -118.59 201.60 -342.72 582.62 -990.46 1683.78-2862.42 4866.12-8272.4014063.08********40642.31 ******************************** Compare Bill Dunlap TIBCO Software wdunlap tibco.comOn Mon, Jul 22, 2019 at 12:19 AM Rui Barradas <ruipbarra...@sapo.pt> wrote:Simpler, no loops: print0b <- function(x, len = 10, digits = 2, fill = ""){ n <- length(x) x <- round(x, digits = digits) m <- n %/% len remainder <- n %% len A <- matrix(x[seq_len(len*m)], ncol = len) if(remainder > 0){ A <- rbind(A, c(x[(len*m + 1):n], rep(fill, len*(m + 1) - n))) } A } Hope this helps, Rui Barradas Às 07:47 de 22/07/19, Rui Barradas escreveu:Hello, Maybe something like the following is what you want. I have added an extra argument 'fill' to allow to choose what to print in the end. It's default value is "" making the entire matrix elements characters but it can be NA or 0. print0 <- function(x, len = 10, digits = 2, fill = ""){ n <- length(x) x <- round(x, digits = digits) passes <- n %/% len remainder <- n %% len A <- matrix(fill, nrow = passes + (remainder > 0), ncol = len) for(i in seq_len(passes)){ A[i, ] <- x[(len*(i - 1) + 1):(len*i)] } A[nrow(A), 1:remainder] <- x[(len*passes + 1):n] A } print0(rnorm(23), 10) print0(rnorm(23), 10, fill = 0) Hope this helps, Rui Barradas Às 21:34 de 20/07/19, Steven escreveu:Dear All:Below is what I meant. Procedure print0 allows me to print a vector of length 53 in four rows of 10 plus 1 row of 3 (Ido not like the NA). Thisis silly. I am hoping that there is a candid way to print the matrix. Thank you. Steven Yen === n<-53; x<-runif(n); # x<-round(x,2) print0<-function(x,c=10,digits=2){ # ****************************************** # Print vector in rows of a specified length # ****************************************** n<-length(x) r<-n/c; if(n%%c>0) r<-as.integer(r)+1 y<-rep(NA,r*c) y[1:n]<-x y<-matrix(y,r,c,byrow=T) y<-round(y,digits=digits) y } print0(x,c=10,digits=3) # result # [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] # [1,] 0.576 0.291 0.600 0.515 0.135 0.335 0.296 0.911 0.454 0.696 # [2,] 0.699 0.728 0.442 0.469 0.996 0.539 0.772 0.768 0.652 0.882 # [3,] 0.614 0.228 0.748 0.071 0.788 0.428 0.885 0.722 0.432 0.881 # [4,] 0.422 0.148 0.459 0.870 0.044 0.421 0.282 0.337 0.751 0.579 # [5,] 0.468 0.659 0.446 0.199 0.388 0.576 0.829 0.186 0.823 0.960 # [6,] 0.880 0.944 0.709 NA NA NA NA NA NA NA Steven 於 2019/7/20 下午 02:00 寫道:Is there a convenient way to print a vector into rows of a specifiedcolumn length? What I need is to print in the old FORTRAN format, viz.,format(10F8.2) which would print, for instance, a vector of 25 into two rows of 10 plus an incomplete row of 5. I managed to write a procedure for that task, as shown below (except that I prefer simply blanks rather thanthe NA). I am too embarrassed to even show the procedure. In short, Ilike to print in the above FORTRAN format. Thank you. ----[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 0.66 0.26 0.82 0.73 0.13 0.05 0.56 0.67 0.74 0.87 [2,] 0.91 0.25 0.40 0.39 0.50 0.890.07 0.84 0.14 0.75 [3,] 0.38 0.08 0.86 0.97 0.56 NA NA NA NA NA[[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.[[alternative HTML version deleted]]
______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
