(*EnestenotebookhacemoselhistogramadelIPRyleajustamosunadistribuciondeprobabilidad*)
SetDirectory[NotebookDirectory[]];
(***********************************************************************************************)iproc3Er=Flatten[Import["IPR_octa_3Er.csv"]];iproc4Er=Flatten[Import["IPR_octa_4Er.csv"]];iproc5Er=Flatten[Import["IPR_octa_5Er.csv"]];iproc6Er=Flatten[Import["IPR_octa_6Er.csv"]];iproc7Er=Flatten[Import["IPR_octa_7Er.csv"]];iproc8Er=Flatten[Import["IPR_octa_8Er.csv"]];iproc9Er=Flatten[Import["IPR_octa_9Er.csv"]];iproc10Er=Flatten[Import["IPR_octa_10Er.csv"]];iproc11Er=Flatten[Import["IPR_octa_11Er.csv"]];iproc12Er=Flatten[Import["IPR_octa_12Er.csv"]];iproc13Er=Flatten[Import["IPR_octa_13Er.csv"]];iproc14Er=Flatten[Import["IPR_octa_14Er.csv"]];iproc15Er=Flatten[Import["IPR_octa_15Er.csv"]];iproc16Er=Flatten[Import["IPR_octa_16Er.csv"]];iproc17Er=Flatten[Import["IPR_octa_17Er.csv"]];iproc18Er=Flatten[Import["IPR_octa_18Er.csv"]];iproc19Er=Flatten[Import["IPR_octa_19Er.csv"]];iproc20Er=Flatten[Import["IPR_octa_20Er.csv"]];
(*buscamosladistribucionanuestrosdatos*)
res3=FindDistributionParameters[iproc3Er,NormalDistribution[μ,σ],WorkingPrecision20]res4=FindDistributionParameters[iproc4Er,NormalDistribution[μ,σ],WorkingPrecision20]res5=FindDistributionParameters[iproc5Er,NormalDistribution[μ,σ],WorkingPrecision20]res6=FindDistributionParameters[iproc6Er,NormalDistribution[μ,σ],WorkingPrecision20]res7=FindDistributionParameters[iproc7Er,NormalDistribution[μ,σ],WorkingPrecision20]
{μ7.94107×,σ4.5138×}
-8
10
-8
10
{μ7.95925×,σ4.5038×}
-8
10
-8
10
{μ8.09678×,σ4.58923×}
-8
10
-8
10
{μ8.25883×,σ4.74587×}
-8
10
-8
10
{μ8.44308×,σ4.95685×}
-8
10
-8
10
res8=FindDistributionParameters[iproc8Er,NormalDistribution[μ,σ],WorkingPrecision20]res9=FindDistributionParameters[iproc9Er,NormalDistribution[μ,σ],WorkingPrecision20]res10=FindDistributionParameters[iproc10Er,NormalDistribution[μ,σ],WorkingPrecision20]res11=FindDistributionParameters[iproc11Er,NormalDistribution[μ,σ],WorkingPrecision20]res12=FindDistributionParameters[iproc12Er,NormalDistribution[μ,σ],WorkingPrecision20]
{μ8.7062×,σ5.36149×}
-8
10
-8
10
{μ9.01641×,σ5.99658×}
-8
10
-8
10
{μ9.4455×,σ7.09963×}
-8
10
-8
10
{μ1.00732×,σ9.13748×}
-7
10
-8
10
{μ1.11139×,σ1.29625×}
-7
10
-7
10
res13=FindDistributionParameters[iproc13Er,NormalDistribution[μ,σ],WorkingPrecision20]res14=FindDistributionParameters[iproc14Er,NormalDistribution[μ,σ],WorkingPrecision20]res15=FindDistributionParameters[iproc15Er,NormalDistribution[μ,σ],WorkingPrecision20]res16=FindDistributionParameters[iproc16Er,NormalDistribution[μ,σ],WorkingPrecision20]res17=FindDistributionParameters[iproc17Er,NormalDistribution[μ,σ],WorkingPrecision20]
{μ1.31013×,σ2.12383×}
-7
10
-7
10
{μ1.75864×,σ4.41135×}
-7
10
-7
10
{μ3.30857×,σ1.43105×}
-7
10
-6
10
{μ1.22901×,σ8.76404×}
-6
10
-6
10
{μ4.87855×,σ0.000047654}
-6
10
res18=FindDistributionParameters[iproc18Er,NormalDistribution[μ,σ],WorkingPrecision20]res19=FindDistributionParameters[iproc19Er,NormalDistribution[μ,σ],WorkingPrecision20]res20=FindDistributionParameters[iproc20Er,NormalDistribution[μ,σ],WorkingPrecision20]
{μ9.92834×,σ0.000112431}
-6
10
{μ0.0000140522,σ0.000170848}
{μ0.0000174199,σ0.000220682}
(******************************************************************************)
(*Graficamos*)
aj1=Show[Histogram[iproc3Er,40,"ProbabilityDensity"],Plot[1.5PDF[NormalDistribution[μ,σ]/.res3,x],{x,0,Max[iproc3Er]},PlotStyleRed]]
aj2=Show[Histogram[iproc4Er,40,"ProbabilityDensity"],Plot[1.5PDF[NormalDistribution[μ,σ]/.res4,x],{x,0,Max[iproc4Er]},PlotStyleRed]]
aj3=Show[Histogram[iproc5Er,40,"ProbabilityDensity"],Plot[1.5PDF[NormalDistribution[μ,σ]/.res5,x],{x,0,Max[iproc5Er]},PlotStyleRed]]
aj4=Show[Histogram[iproc6Er,40,"ProbabilityDensity"],Plot[1.5PDF[NormalDistribution[μ,σ]/.res6,x],{x,0,Max[iproc6Er]},PlotStyleRed]]
aj5=Show[Histogram[iproc7Er,40,"ProbabilityDensity"],Plot[1.5PDF[NormalDistribution[μ,σ]/.res7,x],{x,0,Max[iproc7Er]},PlotStyleRed]]