[R] how to build a heatmap like this

[escaping]

Hi, everyone,

I am a new R user, and was required to build several heatmaps as shown in attached file. To reach it, I need to normalize the RPKM data  of exon-level expression to generate Z-score using R. The data looks like
mpleId

			
				geneName

			
				start

			
				end

			
				rpkm

		


				PC10-6

			
				KIF1B

			
				10270764

			
				10270936

			
				4.00059

		


				PC10-6

			
				KIF1B

			
				10292308

			
				10292492

			
				35.3023

		


				PC10-6

			
				KIF1B

			
				10316305

			
				10316381

			
				85.9611

		


				PC10-6

			
				KIF1B

			
				10318551

			
				10318730

			
				64.5961
Exons and introns are required to shown on the heatmap.

Does anyone know how to do it? 

Thanks a lot!

Yizi