https://github.com/CNG-LAB/cngopen
Tip revision: b8031194ea8737e7aae34249d9376dbf1259f4e3 authored by Bin Wan on 09 August 2022, 14:49:26 UTC
update_issue_license_openaccess
update_issue_license_openaccess
Tip revision: b803119
heritability_solar_collect_data.txt
#! /bin/sh
# this is a universal ICBM processing script
echo This script will assemble the data from VOXEL_X_Y_Z into X Y Z H2r P format
if [ $# -ne 1 ]
then
echo $1
echo Give the name for the output file e.g. brain.txt
exit
fi
# don't modify after this line
pheno_list=`find . -type d -name "node*"`
# Lets make solar_run file
echo node","name","H2r","rp >> $1
for cur_dir in $pheno_list
do
echo $cur_dir
XX=`echo $cur_dir | awk -F "_" '{print $2}'`
echo $XX
X=`echo $cur_dir`
echo $X
line=`more $cur_dir/polygenic.out | grep "H2r is"`
H2r=`echo $line | awk '{print $3}'`
echo $H2r
rp=`echo $line | awk '{print $6}'`
echo $rp
echo se $se $cur_dir/polygenic.out
echo $XX","$X","$H2r","$rp >> $1
done