Hi, I'm have a problem about the script below, please help.

What is the problem you see?
Where does the panel() function come from? It is not a function in base R
What is the object shape_data? Can you post that or a subset of it? See

FAQ: How to do a minimal reproducible example ( reprex ) for beginners - meta / Guides & FAQs - RStudio Community