Please write a function that takes input of an institution of higher learning and a vector of variables and outputs a multipanel plot showing histograms of those variables with the position of the specified institution marked in the histogram. The plots should be based on the $2015-16$ College Scorecard data.

Choose a list of four quantitative variables that you think would be interesting to plot. It should show plots for all four variables by default, but allow a user to select a smaller subset of the variables. The inputs of your function should also include two variables that specify whether the histograms should include all institutions or only include institutions that match the target institution in the type of control $($public$,$ private$-$not$-$for$-$profit, or for$-$profit$)$ and level $(4-$year, $2-$year, less$-$than$-2-$year$).$

Note that some variables are not defined for all types of institutions. For example, completion rate has separate variables for four$-$year and less$-$than$-$four$-$year institutions. Thus, the histogram for this variable always includes only the same type of institutions as the target institution. Your plots should have descriptive labels and titles. The function should give a descriptive error message if the inputs are not correct $($e$.$g$.$ the institution name or the variables are not known$).$

I have attempted the code in R below but I believe it has some issues still. The output of histograms does not show the frequency well and the institution selected at the end 'University of Alabama' is not visible as specified in the question. Please help me resolve these issues and ensure I am answering the question to completion in R$.$ After this, please replicate the exact same function in Python and provide code.

#### Function

$```\{$r$\}$

plot$\_$histogram function$($INSTNM$,$ variables $=$ c$("$TUITIONFEE$\_$IN$",$ "TUITIONFEE$\_$OUT", "AVGFACSAL", "FAMINC"$),$

include$\_$all $=$ TRUE, control$\_$type $=$ NULL, level $=$ NULL$)\{$

# Read the dataset

collegedata read.csv$("$MERGED$2015\_16\_$PP$.$csv$",$na$.$strings$=$c$("","$NA$",$ "NULL","PrivacySuppressed"$))$

# Check if institution$\_$name exists in the dataset

if $(!$INSTNM $\%$in$\%$ collegedata$INSTNM$)\{$

stop$("$Institution name not found in the dataset."$)$

$\}$

# Check if variables are valid

valid$\_$variables c$("$TUITIONFEE$\_$IN$",$ "TUITIONFEE$\_$OUT", "AVGFACSAL", "FAMINC"$)$

if $(!$all$($variables $\%$in$\%$ valid$\_$variables$))\{$

stop$("$Invalid variable$($s$)$ specified."$)$

$\}$

# Filter data based on control type and level if include$\_$all is FALSE

if $(!$include$\_$all$)\{$

if $($is$.$null$($control$\_$type$)||$ is$.$null$($level$))\{$

stop$("$Both control$\_$type and level must be specified when include$\_$all is FALSE."$)$

$\}$

collegedata subset$($collegedata$,$ CONTROL $==$ control$\_$type & LEVEL $==$ level$)$

$\}$

# Subset data for the specified institution

institution$\_$data subset$($collegedata$,$ INSTNM $==$ INSTNM$)$

# Check if institution exists

if $($nrow$($institution$\_$data$)==0)\{$

stop$("$Institution not found in the specified control type and level."$)$

$\}$

# Adjust plot margins

par$($mar $=$ c$(3,4,2,2))$

# Create multipanel plot

par$($mfrow $=$ c$($length$($variables$),1))$

for $($variable in variables$)\{$

# Check if variable exists in the dataset

if $(!($variable $\%$in$\%$ names$($collegedata$)))\{$

warning$($paste$("$Variable$",$ variable, "not found in the dataset. Skipping..."$))$

next

$\}$

hist$\_$data collegedata$[[$variable$]]$

institution$\_$value institution$\_$data$[[$variable$]]$

# Plot histogram

hist$($hist$\_$data, main $=$ paste$("$Histogram of$",$ variable$),$ xlab $=$ variable,

ylab $=$ "Frequency", col $=$ "lightblue", border $=$ "black", ylim $=$ c$(0,$ max$($hist$($hist$\_$data, plot $=$ FALSE$)$$counts$)*2\mathrm{.}1))$

# Add vertical line for institution value

abline$($v $=$ institution$\_$value, col $=$ "blue", lwd $=2)$

$\}$

# Reset plotting parameters

par$($mfrow $=$ c$(1,1))$

par$($mar $=$ c$(5,4,4,2)+0\mathrm{.}1)$ # Reset plot margins to default

$\}$

$```$

#### Testing the plot$\_$histogram function

$```\{$r$\}$

hist plot$\_$histogram$("$University of Alabama", include$\_$all $=$ TRUE$)$

Please do all of this in R and Python as simply as possible, without the use of any complex code and provide annotation which helps explain