2016-09-03 100 views
1

我使用写入命令将xgBoost的重要性矩阵打印到日志中(写入与文件连接一起使用,并将其指向stderr以及)。下面是我使用的命令:打印日志时的不同间距

importance_matrix <- xgb.importance(names, model=bst) 
write("The top 30 variables are:",stderr()) 
write(paste0("Feature",'\t','\t','Gain','\t','Cover','\t','Frequency'),stderr()) 
write(t(as.matrix(importance_matrix[1:30,])),sep="\t",ncolumns = length(names(importance_matrix)),stderr()) 

输出进来格式:

Feature   Gain Cover Frequency 
pctTillDate 0.560359696 0.1314074664 0.024278250 
colr_per 0.183149483 0.0962457545 0.049618673 
date 0.050528297 0.1143752021 0.066395735 
GREG_D 0.025648433 0.0381476142 0.018070143 
LNGTD_I 0.020346020 0.0485235001 0.101322109 
LATTD_I 0.019241497 0.0421892270 0.093867103 

,这使得它(在日志笨拙得多比出现在这里SO)看起来有点笨拙。因此,为了使它更好看,我想改变t(as.matrix(importance_matrix[1:30,])),sep="\t"的最后一行,这样第一个sep将是2个选项卡('\ t','\ t')并休息单个选项卡('\ t');而不是当前的统一间距。简单但搜索并没有任何想法。有什么建议么?

+0

制表符分隔文件将看起来像锯齿输出由于内容的长度单元格,但放心等量的标签位于列之间。请注意:文本文件仅仅是用于数据迁移或存档的转储输出,因此演示文稿不是本意。 – Parfait

+0

@Parfait我同意它不是为了演示,但在我的情况下,我做了400+的mapper只有hadoop流,后来想要将所有的YARN日志合并成一个日志。那么我想保留重要参数的可读性。否则400多个日志本身很难通过,它们笨拙只会增加分析器的痛苦。 – abhiieor

回答

2

考虑填充的列名,并用空格来使每一个到第一列的最大字符大小矩阵的第一个字符列:

write.table(importance_matrix, sep="\t", row.names = FALSE, quote = FALSE) 
# Feature Gain Cover Frequency 
# pctTillDate 0.56035970 0.13140747 0.02427825 
# colr_per 0.18314948 0.09624575 0.04961867 
# date 0.05052830 0.11437520 0.06639573 
# GREG_D 0.02564843 0.03814761 0.01807014 
# LNGTD_I 0.02034602 0.04852350 0.10132211 
# LATTD_I 0.01924150 0.04218923 0.09386710 

new_matrix <- importance_matrix 

# FIRST COLUMN LARGEST CHAR LENGTH 
charmax <- max(nchar(new_matrix[,1])) 

# PAD COLUMN HEADERS 
colnames(new_matrix) <- lapply(1:ncol(new_matrix), function(i) 
     paste0(colnames(new_matrix)[i], 
       paste(rep(" ", charmax - nchar(colnames(new_matrix)[i])), collapse="")) 
) 

# PAD FIRST COLUMN 
new_matrix[,1] <- sapply(1:nrow(new_matrix), function(i) 
     paste0(new_matrix[i,1], 
       paste(rep(" ", charmax - nchar(new_matrix[i,1])), collapse="")) 
) 

write.table(new_matrix, sep="\t", row.names = FALSE, quote = FALSE) 
# Feature  Gain  Cover  Frequency 
# pctTillDate 0.56035970 0.13140747 0.02427825 
# colr_per  0.18314948 0.09624575 0.04961867 
# date   0.05052830 0.11437520 0.06639573 
# GREG_D  0.02564843 0.03814761 0.01807014 
# LNGTD_I  0.02034602 0.04852350 0.10132211 
# LATTD_I  0.01924150 0.04218923 0.09386710