A Percolation Experiment.nb

(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 11.3' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     26789,        651]
NotebookOptionsPosition[     23793,        589]
NotebookOutlinePosition[     24383,        609]
CellTagsIndexPosition[     24340,        606]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["A Percolation Experiment", "Title",
 CellChangeTimes->{{3.735401190130281*^9, 
  3.735401207854809*^9}},ExpressionUUID->"54110861-f639-45a8-8036-\
360d6e82038a"],

Cell["Alan Wang, Eric Shao, Ian Klein, Ziyang Gao", "Chapter",
 CellChangeTimes->{{3.7354012490351543`*^9, 
  3.735401302173149*^9}},ExpressionUUID->"5ef4eb4a-8dfe-4b68-95c3-\
44175b0db53d"],

Cell[CellGroupData[{

Cell["15/5/2018", "Chapter",
 CellChangeTimes->{{3.735401214227606*^9, 3.735401220277458*^9}, 
   3.735401316302207*^9},ExpressionUUID->"65893379-dc3c-4c4c-84a9-\
6433669c7719"],

Cell[CellGroupData[{

Cell["Objective", "Section",
 CellChangeTimes->{{3.7354014807834578`*^9, 3.735401531025219*^9}, {
  3.735401607806436*^9, 3.735401611395979*^9}, {3.735401837663166*^9, 
  3.735401844075989*^9}, {3.7354024076090937`*^9, 
  3.73540241643454*^9}},ExpressionUUID->"8109fa25-7a21-4f3d-b444-\
f05df3cc3222"],

Cell[TextData[{
 StyleBox["Our object is to study the phenomenon of percolation process \
through experiment with a two-dimensional lattice made of identical small \
squares (\[OpenCurlyDoubleQuote]sites\[CloseCurlyDoubleQuote]). Through \
measuring the resistance of the lattice conductive paper with different \
conductance, we try to derive the relationship between ",
  FontSize->18],
 StyleBox["p",
  FontSize->18,
  FontSlant->"Italic"],
 StyleBox[", the concentration of the sites of the lattice and \
\[CapitalOmega], the finite effective resistance.",
  FontSize->18]
}], "Text",
 CellChangeTimes->{{3.7354018490854197`*^9, 3.7354019002828703`*^9}, {
   3.735401978159585*^9, 3.735401992229109*^9}, {3.735402035061934*^9, 
   3.735402043835126*^9}, {3.735402094604944*^9, 3.735402153830616*^9}, {
   3.7354022018940353`*^9, 3.7354022306802692`*^9}, {3.7354023000524483`*^9, 
   3.735402378303405*^9}, {3.7354024430048523`*^9, 3.735402517329321*^9}, {
   3.735402554340952*^9, 3.7354026789778748`*^9}, 3.73540283846567*^9, 
   3.7354029503457108`*^9, {3.735403054194621*^9, 3.735403056592243*^9}, {
   3.735580949122334*^9, 
   3.735580957945977*^9}},ExpressionUUID->"cad76a21-b6d4-4a07-ae5f-\
f4e5e0cdfbe5"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Procedure", "Section",
 CellChangeTimes->{{3.7354015329870872`*^9, 3.735401546040732*^9}, {
  3.735401615684142*^9, 
  3.735401626147417*^9}},ExpressionUUID->"ca632133-c69b-4241-85af-\
26a7bbad3f0b"],

Cell[TextData[{
 StyleBox["The description below are the steps we followed through our \
experiment: \n1. Measure the resistance between the two ends of the \
conductive paper using a digital ohmmeter.\n2. Use Mathematica to generate a \
large set of x and y coordinates for sites that will be cut out.\n3. Cut out \
small squares of conductive paper one by one. Since the paper is 20 cm by 28 \
cm, we measure the resistance for every 6 squares cut off (about 1 % of the \
total area).\n4. Repeat step 4 until the resistance goes to infinity \
(conductance goes to 0)\n5. Because we expect an exponential dependence of \
the conductance to ",
  FontSize->18],
 StyleBox["(P - Pc)",
  FontSize->18,
  FontSlant->"Italic"],
 StyleBox[", we plot ",
  FontSize->18],
 StyleBox["log",
  FontSize->18,
  FontSlant->"Italic"],
 StyleBox["(Conductance) (conductance is the inverse of the measured \
resistance) against ",
  FontSize->18],
 StyleBox["log(P - Pc)",
  FontSize->18,
  FontSlant->"Italic"],
 StyleBox[".The slope of this plot should give us a good approximation of t \
(critical exponent).",
  FontSize->18]
}], "Text",
 CellChangeTimes->{{3.73540276161692*^9, 3.735402863858609*^9}, {
  3.735402896320311*^9, 3.735403002468384*^9}, {3.7354030455278597`*^9, 
  3.7354030472331047`*^9}, {3.73540308322407*^9, 3.735403123823505*^9}, {
  3.735403209534327*^9, 3.7354032484484587`*^9}, {3.735574627731543*^9, 
  3.7355746306614227`*^9}, {3.735576019229609*^9, 3.735576029964684*^9}, {
  3.735576073165058*^9, 3.735576088489911*^9}, {3.735580987226368*^9, 
  3.7355810040538273`*^9}},ExpressionUUID->"40588d5b-62ab-4a67-92ba-\
39b47998541b"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Data and Plot", "Section",
 CellChangeTimes->{{3.73540166426609*^9, 
  3.73540167701036*^9}},ExpressionUUID->"359c02e6-bee8-4efe-9a44-\
1aeff4cd2473"],

Cell[TextData[StyleBox["Our first plot is the the relationship between the \
resistance(\[CapitalOmega]) and P. ",
 FontSize->18]], "Text",
 CellChangeTimes->{{3.735575716782013*^9, 3.7355757522591333`*^9}, {
  3.73557578538304*^9, 3.7355758563468447`*^9}, {3.735575963375819*^9, 
  3.735575983516292*^9}, {3.735578064692472*^9, 3.7355780720382967`*^9}, {
  3.735581016320281*^9, 
  3.7355810274025717`*^9}},ExpressionUUID->"f0c95c33-c5f9-4de6-a090-\
312e26086b7d"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"init", " ", "=", " ", 
   RowBox[{"{", 
    RowBox[{
    "108", ",", "116", ",", "127", ",", "135", ",", "149", ",", "165", ",", 
     "174", ",", "187", ",", "195", ",", "204", ",", "216", ",", "224", ",", 
     "236", ",", "246", ",", "254", ",", "266", ",", "278", ",", "290", ",", 
     "303", ",", "312", ",", "323", ",", "338", ",", "355", ",", "392", ",", 
     "437", ",", "480", ",", "560"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"dataSet", " ", "=", " ", 
   RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "0"}], ",", 
   RowBox[{"i", "<", 
    RowBox[{"Length", "[", "init", "]"}]}], ",", 
   RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"AppendTo", "[", 
     RowBox[{"dataSet", ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{"560", "-", 
           RowBox[{"6", "*", "i"}]}], ")"}], "/", "560"}], ",", 
        RowBox[{"init", "[", 
         RowBox[{"[", "i", "]"}], "]"}]}], "}"}]}], "]"}], ";"}]}], 
  "\[IndentingNewLine]", "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"P", " ", "=", 
   RowBox[{"dataSet", "[", 
    RowBox[{"[", 
     RowBox[{"All", ",", "1"}], "]"}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.7355725495529413`*^9, 3.735572694522958*^9}, {
  3.735572733096553*^9, 3.735572821304522*^9}, {3.7355731301186037`*^9, 
  3.73557314323393*^9}, {3.7355751657784357`*^9, 3.7355751718107224`*^9}, {
  3.735575656965296*^9, 3.7355756673244963`*^9}, {3.735577459454273*^9, 
  3.735577464394644*^9}, {3.735577597900056*^9, 3.735577673225155*^9}, {
  3.735577914976804*^9, 3.7355779318576202`*^9}, {3.7355793801763372`*^9, 
  3.735579441696149*^9}},
 CellLabel->"In[1]:=",ExpressionUUID->"0736cf5d-b83f-4a3b-bf0f-2a5e6b816b78"],

Cell[CellGroupData[{

Cell[BoxData[
 StyleBox[
  RowBox[{"plot1", "=", 
   RowBox[{"ListPlot", "[", 
    RowBox[{"dataSet", ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"RGBColor", "[", 
         RowBox[{"1", ",", "0", ",", "0"}], "]"}], ",", 
        RowBox[{"PointSize", "[", "0.01", "]"}]}], "}"}]}], ",", 
     RowBox[{"Frame", "\[Rule]", "True"}], ",", 
     RowBox[{"FrameLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"\"\<P\>\"", ",", "\"\<Resistance(Ohm)\>\""}], "}"}]}], ",", 
     RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
     RowBox[{"PlotLabel", "\[Rule]", 
      RowBox[{"Style", "[", 
       RowBox[{"\"\<Resistance versus P\>\"", ",", 
        RowBox[{"FontFamily", "\[Rule]", "\"\<Helvetica\>\""}], ",", 
        RowBox[{"FontSize", "\[Rule]", "12"}], ",", 
        RowBox[{"FontWeight", "\[Rule]", "\"\<Bold\>\""}]}], "]"}]}], ",", 
     RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", 
     RowBox[{"ImageSize", "\[Rule]", "500"}]}], "]"}]}],
  FontSize->14]], "Input",
 CellChangeTimes->{
  3.7355731619147577`*^9, {3.735580855179528*^9, 3.73558087165746*^9}, {
   3.735581032810854*^9, 3.735581096200323*^9}},
 CellLabel->"In[5]:=",ExpressionUUID->"f2f2b58c-6a41-4d6e-b0e1-0af7a101a5fe"],

Cell[BoxData[
 FormBox[
  GraphicsBox[{{}, {{}, 
     {RGBColor[1, 0, 0], PointSize[0.01], AbsoluteThickness[1.6], 
      StyleBox[PointBox[CompressedData["
1:eJxTTMoPSmViYGCQAmIQ3WK1otlqxXt7BjCIdmAvDQAiGD/WoXv9j671P95B
+AfiHfhfLeB7tQDKf5DgMEnFE4ig/AVJDmLxH0XjP76F8lMcZs0EASj/QKqD
7BVHIILyE9IdFvKBDHwD5Wc4qIKMmwTlN2Q6gFzXbAXlM2Q76Ox7pL3v0Wso
P8dhA8h5XVB+Q66DMRhA+QfyHLbn3t6We/sVlJ/vYA02EMpfUOAANk4byk8o
dHACO/AlhK9Q5HA0vAaIoPwPRQ5g56lA+Q3FDmfPgMALCN+gxCEAHIBQvkKp
A9i7sjD5MocIkHFHn0P1VziAXLctF8oPqHYAB58olM9Q5wAAbJWI2Q==
        "]],
       FontFamily->"Arial"]}, {}}, {}, {}, {}, {}},
   AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
   Axes->{True, True},
   AxesLabel->{None, None},
   AxesOrigin->{0.7172433035714286, 0},
   BaseStyle->{FontFamily -> "Arial"},
   DisplayFunction->Identity,
   Frame->{{True, True}, {True, True}},
   FrameLabel->{{
      FormBox["\"Resistance(Ohm)\"", TraditionalForm], None}, {
      FormBox["\"P\"", TraditionalForm], None}},
   FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
   GridLines->{Automatic, Automatic},
   GridLinesStyle->Directive[
     GrayLevel[0.5, 0.4]],
   ImagePadding->All,
   ImageSize->500,
   Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
         (Identity[#]& )[
          Part[#, 1]], 
         (Identity[#]& )[
          Part[#, 2]]}& ), "CopiedValueFunction" -> ({
         (Identity[#]& )[
          Part[#, 1]], 
         (Identity[#]& )[
          Part[#, 2]]}& )}},
   PlotLabel->FormBox[
     StyleBox[
     "\"Resistance versus P\"", FontFamily -> "Helvetica", FontSize -> 12, 
      FontWeight -> "Bold", StripOnInput -> False], TraditionalForm],
   PlotRange->{{0.7214285714285714, 0.9892857142857143}, {0, 480.}},
   PlotRangeClipping->True,
   PlotRangePadding->{{
      Scaled[0.02], 
      Scaled[0.02]}, {
      Scaled[0.02], 
      Scaled[0.05]}},
   Ticks->{Automatic, Automatic}], TraditionalForm]], "Output",
 CellChangeTimes->{
  3.735577471259715*^9, 3.73557762365224*^9, 3.735577676908205*^9, 
   3.735579453444016*^9, 3.7355808740696363`*^9, 3.7355810452310667`*^9, {
   3.7355810887582083`*^9, 3.7355811014495783`*^9}, 3.7355813543229017`*^9},
 CellLabel->"Out[5]=",ExpressionUUID->"7c26f04e-49ad-42c5-a8dc-d8738f2aa844"]
}, Open  ]],

Cell[TextData[StyleBox["\nOur second plot is the relationship between  log \
(Conductance) (conductance is the inverse of the measured resistance) against \
log (P - Pc).",
 FontSize->18]], "Text",
 CellChangeTimes->{{3.735578077321547*^9, 3.73557808412878*^9}, 
   3.735578119697363*^9},ExpressionUUID->"2771cb0f-1b66-46f1-93f7-\
ad546527184e"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"Pc", "=", 
  RowBox[{
   RowBox[{"1", "-", 
    RowBox[{
     RowBox[{"Length", "[", "init", "]"}], "*", 
     RowBox[{"6", "/", 
      RowBox[{"(", 
       RowBox[{"20", "*", "28"}], ")"}]}]}]}], "//", 
   "N"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"conductanceSet", "=", 
   RowBox[{"Log", "[", 
    RowBox[{"1", "/", "init"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"PPcSet", "=", 
   RowBox[{"Log", "[", 
    RowBox[{"P", "-", "Pc"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"dataSet2", " ", "=", " ", 
   RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"For", "[", 
   RowBox[{
    RowBox[{"i", "=", "1"}], ",", 
    RowBox[{"i", "\[LessEqual]", 
     RowBox[{"Length", "[", "conductanceSet", "]"}]}], ",", 
    RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"AppendTo", "[", 
      RowBox[{"dataSet2", ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"PPcSet", "[", 
          RowBox[{"[", "i", "]"}], "]"}], ",", 
         RowBox[{"conductanceSet", "[", 
          RowBox[{"[", "i", "]"}], "]"}]}], "}"}]}], "]"}], ";"}]}], 
   "\[IndentingNewLine]", "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.735574684089658*^9, 3.735574692124015*^9}, {
   3.735574885118393*^9, 3.735574986426917*^9}, {3.735575117995585*^9, 
   3.735575134109962*^9}, {3.735575240784678*^9, 3.735575292029511*^9}, {
   3.735575330574944*^9, 3.7355753391242847`*^9}, {3.7355754909239902`*^9, 
   3.7355754934666233`*^9}, {3.735575674053743*^9, 3.7355756752147217`*^9}, {
   3.735577502740513*^9, 3.735577532432534*^9}, 3.735577576958418*^9, {
   3.7355776893073683`*^9, 3.7355776963060007`*^9}, 3.735578043038404*^9, {
   3.735579511462364*^9, 3.7355796575651608`*^9}, {3.735580161661293*^9, 
   3.7355801804131203`*^9}, {3.735580210773408*^9, 3.735580246592834*^9}, {
   3.735580829499632*^9, 3.735580830707961*^9}},
 CellLabel->"In[6]:=",ExpressionUUID->"30db69ad-7572-43fc-ae8b-201d526b0eb9"],

Cell[BoxData[
 FormBox["0.7107142857142857`", TraditionalForm]], "Output",
 CellChangeTimes->{
  3.735577476050116*^9, {3.735577511006721*^9, 3.735577534447476*^9}, 
   3.73557757771798*^9, 3.7355776266717873`*^9, {3.735577679328392*^9, 
   3.735577697295014*^9}, 3.7355780448482027`*^9, 3.735579478960401*^9, 
   3.735579514299122*^9, {3.735579581450322*^9, 3.735579606739667*^9}, 
   3.7355796586628523`*^9, {3.735580164374475*^9, 3.7355801864096193`*^9}, {
   3.735580216794422*^9, 3.735580248662573*^9}, 3.735580832137114*^9, 
   3.735580880215074*^9, 3.735581354416181*^9},
 CellLabel->"Out[6]=",ExpressionUUID->"6c38414d-7557-477b-9cbe-dde142a3837e"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 StyleBox[
  RowBox[{
   RowBox[{"plot2", "=", 
    RowBox[{"ListPlot", "[", 
     RowBox[{"dataSet2", ",", 
      RowBox[{"PlotStyle", "\[Rule]", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"RGBColor", "[", 
          RowBox[{"1", ",", "0", ",", "0"}], "]"}], ",", 
         RowBox[{"PointSize", "[", "0.01", "]"}]}], "}"}]}], ",", 
      RowBox[{"Frame", "\[Rule]", "True"}], ",", 
      RowBox[{"FrameLabel", "\[Rule]", 
       RowBox[{"{", 
        RowBox[{"\"\<Log(P-Pc)\>\"", ",", "\"\<Log Conductance\>\""}], 
        "}"}]}], ",", 
      RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
      RowBox[{"PlotLabel", "\[Rule]", 
       RowBox[{"Style", "[", 
        RowBox[{"\"\<Log Conductance versus Log(P-Pc)\>\"", ",", 
         RowBox[{"FontFamily", "\[Rule]", "\"\<Helvetica\>\""}], ",", 
         RowBox[{"FontSize", "\[Rule]", "12"}], ",", 
         RowBox[{"FontWeight", "\[Rule]", "\"\<Bold\>\""}]}], "]"}]}], ",", 
      RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", 
      RowBox[{"ImageSize", "\[Rule]", "500"}]}], "]"}]}], 
   "\[IndentingNewLine]"}],
  FontSize->14]], "Input",
 CellChangeTimes->{{3.735575413288041*^9, 3.735575413404303*^9}, 
   3.73557569947018*^9, {3.7355798708531647`*^9, 3.735580040148498*^9}, {
   3.735580083984898*^9, 3.735580086586328*^9}, {3.735580197182508*^9, 
   3.735580197788735*^9}, {3.7355806571313353`*^9, 3.735580822468248*^9}},
 CellLabel->"In[11]:=",ExpressionUUID->"ab6d307e-bd00-4f46-b66b-82f4fd3bddd1"],

Cell[BoxData[
 FormBox[
  GraphicsBox[{{}, {{}, 
     {RGBColor[1, 0, 0], PointSize[0.01], AbsoluteThickness[1.6], 
      StyleBox[PointBox[CompressedData["
1:eJwBwQE+/iFib1JlAgAAABsAAAACAAAASDLjOW/Y87/EUMCbgLoSwM0FCtQE
c/S/i42lJa0DE8ABIQ+6qhP1v0IuYJRyYBPAduHFsN+69b+VJf1ZAJ8TwIfT
Rq4yafa/LecbgAoEFMDZr267RR/3v7U/Zvh8bBTABd4KgdHd978Qh3Fk36IU
wER0ArOppfi/Y9u8vKfsFMAP1eCdwnf5vzhqPzKNFxXAWxgTNThV+r/7Dmc/
wUUVwHM89SVXP/u/Apif+0iAFcCKx/WrqDf8v8v66IqGpRXAKatPNwJA/b/1
t4W59toVwBfEOnyaWv6/J95RoXUFFsDHIodTJor/v4B1P/Q6JhbAN38hmABp
AMBqN24OgFUWwGfmdZ0zGwHAUFL2G6+CFsCfyD+ZZd4BwJqpwUP1rRbAqRpI
2iy2AsDz+sbF3NoWwEJyuRVlpwPApdzh09X4FsCJ8Nv93bgEwFOZ/+NQHBfA
sw3gV5H0BcCP0xCLzEoXwB1X/lj2aQfAFjogcgx9F8C+AHjV9TIJwOUCt3aS
4hfALpyeFyKAC8COM+sB2lEYwECPNpWGvg7ARQgD/PSxGMDejnqqCyUSwNoW
BanOTxnAUnnbOQ==
        "]],
       FontFamily->"Arial"]}, {}}, {}, {}, {}, {}},
   AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
   Axes->{True, True},
   AxesLabel->{None, None},
   AxesOrigin->{-1.1888429974654933`, -4.599840949293971},
   BaseStyle->{FontFamily -> "Arial"},
   DisplayFunction->Identity,
   Frame->{{True, True}, {True, True}},
   FrameLabel->{{
      FormBox["\"Log Conductance\"", TraditionalForm], None}, {
      FormBox["\"Log(P-Pc)\"", TraditionalForm], None}},
   FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
   GridLines->{Automatic, Automatic},
   GridLinesStyle->Directive[
     GrayLevel[0.5, 0.4]],
   ImagePadding->All,
   ImageSize->500,
   Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
         (Identity[#]& )[
          Part[#, 1]], 
         (Identity[#]& )[
          Part[#, 2]]}& ), "CopiedValueFunction" -> ({
         (Identity[#]& )[
          Part[#, 1]], 
         (Identity[#]& )[
          Part[#, 2]]}& )}},
   PlotLabel->FormBox[
     StyleBox[
     "\"Log Conductance versus Log(P-Pc)\"", FontFamily -> "Helvetica", 
      FontSize -> 12, FontWeight -> "Bold", StripOnInput -> False], 
     TraditionalForm],
   PlotRange->{{-4.536177314501144, -1.240340448496811}, {-6.327936783729195, \
-4.68213122712422}},
   PlotRangeClipping->True,
   PlotRangePadding->{{
      Scaled[0.02], 
      Scaled[0.02]}, {
      Scaled[0.05], 
      Scaled[0.05]}},
   Ticks->{Automatic, Automatic}], TraditionalForm]], "Output",
 CellChangeTimes->{{3.7355807515245533`*^9, 3.7355808239359293`*^9}, 
   3.735580881748811*^9, 3.735581354504572*^9},
 CellLabel->"Out[11]=",ExpressionUUID->"ef5b41cf-6d06-4429-ad8e-0976c062aa4b"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 StyleBox[
  RowBox[{"model1", "=", 
   RowBox[{"LinearModelFit", "[", 
    RowBox[{"dataSet2", ",", 
     RowBox[{"{", "x", "}"}], ",", "x"}], "]"}]}],
  FontSize->14]], "Input",
 CellLabel->"In[12]:=",ExpressionUUID->"d55b034d-8d8a-4c86-b5b8-d155483e6341"],

Cell[BoxData[
 FormBox[
  TagBox[
   RowBox[{"FittedModel", "[", 
    TagBox[
     PanelBox[
      TagBox[
       RowBox[{
        RowBox[{"0.479111432099796`", " ", "x"}], "-", "4.45123979995591`"}],
       Short[#, 2]& ],
      FrameMargins->5],
     Editable -> False], "]"}],
   InterpretTemplate[
   FittedModel[{
     "Linear", {-4.45123979995591, 
      0.479111432099796}, {{$CellContext`x}, {1, $CellContext`x}}, {0, 0}}, {{
     1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 
     1., 1., 1., 1., 1., 1., 1., 1., 
     1.}}, {{-1.240340448496811, -Log[108]}, {-1.2780807764796578`, -
       Log[116]}, {-1.3173014896329394`, -Log[127]}, {-1.3581234841531944`, -
       Log[135]}, {-1.4006830985719902`, -Log[149]}, {-1.4451348611428243`, -
       Log[165]}, {-1.491654876777717, -Log[174]}, {-1.5404450409471488`, -
       Log[187]}, {-1.5917383353346997`, -Log[195]}, {-1.6458055566049754`, -
       Log[204]}, {-1.7029639704449238`, -Log[216]}, {-1.763588592261359, -
       Log[224]}, {-1.8281271133989299`, -Log[236]}, {-1.897119984885881, -
       Log[246]}, {-1.9712279570396036`, -Log[254]}, {-2.05127066471314, -
       Log[266]}, {-2.138282041702769, -Log[278]}, {-2.2335922215070947`, -
       Log[290]}, {-2.3389527371649206`, -Log[303]}, {-2.456735772821305, -
       Log[312]}, {-2.590267165445827, -Log[323]}, {-2.744417845273085, -
       Log[338]}, {-2.926739402067041, -Log[355]}, {-3.1498829533812502`, -
       Log[392]}, {-3.43756502583303, -Log[437]}, {-3.8430301339411983`, -
       Log[480]}, {-4.536177314501144, -Log[560]}}, CompressedData["
1:eJxTTMoPSmViYGCQBmIQDQEf7D2MHlvm3/i8H8Y/y8p1haX4C5zPqMi/a5Xw
Vzi/7OHRDfd3Ifjtl93WGWV+g/Nvrs/b7Sr/Hc5nvcfVePEugu9SwrR55dIf
cD7/1QdzD5X/hPOjJYRNLUJ/wfnFNl9Vw+1/w/ldx7+uXmH+B87XXO1vzuTw
F84XP2JVMyvqH5x/XKk9WK3rP5xvXq84gyGT4QCMn/6sdK6xNCOcP/+E/czU
ewj+SimPWzrbmOB8p6KdoqnLmeH8zg+3/97dwQLnb+Z9ED7xCyucLxv+L+Jb
Jjucv4+h4upXI044X2/OPHGlBm4436HfbGrbPj44/15f1SpuVaEDAKXVkS4=

     "], 
     Function[Null, 
      Internal`LocalizedBlock[{$CellContext`x}, #], {HoldAll}]]& ],
   Editable->False,
   SelectWithContents->True,
   Selectable->True], TraditionalForm]], "Output",
 CellChangeTimes->{3.735580112573188*^9, 3.735580259762125*^9, 
  3.735581354651593*^9},
 CellLabel->"Out[12]=",ExpressionUUID->"0296acb8-6b69-47bb-8ae2-977af5ab82a8"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"model1", "[", "\"\<ParameterTable\>\"", "]"}]], "Input",
 CellChangeTimes->{{3.7355803410259047`*^9, 3.735580371177719*^9}},
 CellLabel->"In[13]:=",ExpressionUUID->"432511a7-6fbd-4573-b856-5d8856c9f5d4"],

Cell[BoxData[
 FormBox[
  StyleBox[
   TagBox[GridBox[{
      {"\<\"\"\>", "\<\"Estimate\"\>", "\<\"Standard Error\"\>", "\<\"t\
\[Hyphen]Statistic\"\>", "\<\"P\[Hyphen]Value\"\>"},
      {"1", 
       RowBox[{"-", "4.45123979995591`"}], "0.08636355317921623`", 
       RowBox[{"-", "51.54072101131568`"}], "6.609188172841369`*^-27"},
      {"x", "0.479111432099796`", "0.03757355369052019`", 
       "12.75129406300144`", "1.9329182816090172`*^-12"}
     },
     AutoDelete->False,
     GridBoxAlignment->{"Columns" -> {{Left}}, "Rows" -> {{Automatic}}},
     GridBoxDividers->{
      "ColumnsIndexed" -> {2 -> GrayLevel[0.7]}, 
       "RowsIndexed" -> {2 -> GrayLevel[0.7]}},
     GridBoxItemSize->{"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}},
     GridBoxSpacings->{
      "ColumnsIndexed" -> {2 -> 1}, "RowsIndexed" -> {2 -> 0.75}}],
    "Grid"], "DialogStyle",
   StripOnInput->False], TraditionalForm]], "Output",
 CellChangeTimes->{3.7355803717599792`*^9, 3.735581354731284*^9},
 CellLabel->"Out[13]=",ExpressionUUID->"5aa13805-063a-48e6-9f69-59f53851d9e5"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Analysis and Conclusion ", "Section",
 CellChangeTimes->{{3.735401697042487*^9, 3.7354017180726624`*^9}, {
  3.735573646265518*^9, 
  3.735573652301127*^9}},ExpressionUUID->"784aba9a-885f-405e-bbe5-\
662180d4a31b"],

Cell[TextData[{
 StyleBox["From the first plot we could see an overall inverse relationship \
between the resistance, \[CapitalOmega] and and concentration P(number of \
square remains divided by the total number of sites of the conductivity \
paper). From the second plot we can conclude that ",
  FontSize->18],
 Cell[BoxData[
  FormBox[
   StyleBox[
    SubscriptBox["P", "c"],
    FontSize->16], TraditionalForm]],
  FormatType->"TraditionalForm",ExpressionUUID->
  "73b8e642-7c6f-4f13-941e-b726d2c846d8"],
 StyleBox[" is about 0.7, and from the linear model fit we got t value equals \
to 0.497. We derived  ",
  FontSize->18],
 Cell[BoxData[
  FormBox[
   StyleBox[
    SubscriptBox["P", "c"],
    FontSize->16], TraditionalForm]],
  FormatType->"TraditionalForm",ExpressionUUID->
  "24c95fdf-fedb-41ba-908e-d7165c0b1bd8"],
 StyleBox[" from the argument of the Peculation article that there is no weak \
links exist in the paper and cutting off a site belongs to dead ends doesn\
\[CloseCurlyQuote]t affect on our resistance much any more.",
  FontSize->18]
}], "Text",
 CellChangeTimes->{{3.7355718395316553`*^9, 3.7355718941503363`*^9}, {
  3.73557198955083*^9, 3.73557203387987*^9}, {3.735572118214024*^9, 
  3.735572120151482*^9}, {3.735572357829474*^9, 3.7355724796787367`*^9}, {
  3.735572947517291*^9, 3.735572957344202*^9}, {3.735573289919222*^9, 
  3.735573437551996*^9}, {3.735573506510138*^9, 3.735573596035473*^9}, {
  3.735575878020979*^9, 3.735575900265914*^9}, {3.7355797022196007`*^9, 
  3.7355797126044197`*^9}, {3.735580279990769*^9, 3.73558033673519*^9}, {
  3.73558038245286*^9, 3.735580401936748*^9}, {3.73558044058223*^9, 
  3.735580629785145*^9}, {3.735581134203927*^9, 
  3.7355811739938087`*^9}},ExpressionUUID->"5d260513-6e50-4bc3-8c06-\
631a333fb078"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
},
WindowSize->{820, 811},
WindowMargins->{{Automatic, 11}, {Automatic, 0}},
SpellingDictionaries->{"CorrectWords"->{
 "nn", "kab", "len", "Shao", "Cn", "Wildavsky", "Hoshen", "Kopelman", "Pc", 
  "ln", "Pc is about 0.7. And from the beca"}},
FrontEndVersion->"11.3 for Mac OS X x86 (32-bit, 64-bit Kernel) (March 5, \
2018)",
StyleDefinitions->FrontEnd`FileName[{"Report"}, "StandardReport.nb", 
  CharacterEncoding -> "UTF-8"]
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 167, 3, 97, "Title",ExpressionUUID->"54110861-f639-45a8-8036-360d6e82038a"],
Cell[750, 27, 190, 3, 69, "Chapter",ExpressionUUID->"5ef4eb4a-8dfe-4b68-95c3-44175b0db53d"],
Cell[CellGroupData[{
Cell[965, 34, 177, 3, 58, "Chapter",ExpressionUUID->"65893379-dc3c-4c4c-84a9-6433669c7719"],
Cell[CellGroupData[{
Cell[1167, 41, 301, 5, 69, "Section",ExpressionUUID->"8109fa25-7a21-4f3d-b444-f05df3cc3222"],
Cell[1471, 48, 1216, 23, 141, "Text",ExpressionUUID->"cad76a21-b6d4-4a07-ae5f-f4e5e0cdfbe5"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2724, 76, 205, 4, 69, "Section",ExpressionUUID->"ca632133-c69b-4241-85af-26a7bbad3f0b"],
Cell[2932, 82, 1642, 36, 271, "Text",ExpressionUUID->"40588d5b-62ab-4a67-92ba-39b47998541b"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4611, 123, 156, 3, 69, "Section",ExpressionUUID->"359c02e6-bee8-4efe-9a44-1aeff4cd2473"],
Cell[4770, 128, 465, 8, 37, "Text",ExpressionUUID->"f0c95c33-c5f9-4de6-a090-312e26086b7d"],
Cell[5238, 138, 1862, 43, 168, "Input",ExpressionUUID->"0736cf5d-b83f-4a3b-bf0f-2a5e6b816b78"],
Cell[CellGroupData[{
Cell[7125, 185, 1255, 28, 110, "Input",ExpressionUUID->"f2f2b58c-6a41-4d6e-b0e1-0af7a101a5fe"],
Cell[8383, 215, 2254, 54, 360, "Output",ExpressionUUID->"7c26f04e-49ad-42c5-a8dc-d8738f2aa844"]
}, Open  ]],
Cell[10652, 272, 345, 6, 89, "Text",ExpressionUUID->"2771cb0f-1b66-46f1-93f7-ad546527184e"],
Cell[CellGroupData[{
Cell[11022, 282, 2016, 48, 168, "Input",ExpressionUUID->"30db69ad-7572-43fc-ae8b-201d526b0eb9"],
Cell[13041, 332, 656, 10, 37, "Output",ExpressionUUID->"6c38414d-7557-477b-9cbe-dde142a3837e"]
}, Open  ]],
Cell[CellGroupData[{
Cell[13734, 347, 1494, 32, 132, "Input",ExpressionUUID->"ab6d307e-bd00-4f46-b66b-82f4fd3bddd1"],
Cell[15231, 381, 2476, 59, 357, "Output",ExpressionUUID->"ef5b41cf-6d06-4429-ad8e-0976c062aa4b"]
}, Open  ]],
Cell[CellGroupData[{
Cell[17744, 445, 272, 7, 43, "Input",ExpressionUUID->"d55b034d-8d8a-4c86-b5b8-d155483e6341"],
Cell[18019, 454, 2338, 47, 51, "Output",ExpressionUUID->"0296acb8-6b69-47bb-8ae2-977af5ab82a8"]
}, Open  ]],
Cell[CellGroupData[{
Cell[20394, 506, 227, 3, 41, "Input",ExpressionUUID->"432511a7-6fbd-4573-b856-5d8856c9f5d4"],
Cell[20624, 511, 1073, 23, 76, "Output",ExpressionUUID->"5aa13805-063a-48e6-9f69-59f53851d9e5"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[21746, 540, 220, 4, 69, "Section",ExpressionUUID->"784aba9a-885f-405e-bbe5-662180d4a31b"],
Cell[21969, 546, 1784, 38, 169, "Text",ExpressionUUID->"5d260513-6e50-4bc3-8c06-631a333fb078"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}
]
*)