[R] Printing vector
Michael Friendly
|r|end|y @end|ng |rom yorku@c@
Tue Jul 23 19:23:40 CEST 2019
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)
> }
>
> Compare
>> print1((-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.com
>
>
> On Mon, Jul 22, 2019 at 12:19 AM Rui Barradas <ruipbarradas using 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). This
>>>> is 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 specified
>>>>> column 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 than
>>>>> the NA). I am too embarrassed to even show the procedure. In short, I
>>>>> like 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.89
>>>>> 0.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 using 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 using 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 using 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]]
>
More information about the R-help
mailing list