PMMA

PMMA g2 alpha term all terms (t^2 to t^20)

alpha.ngp.plot(filename1="MSD.colmean.matrix.1.txt",filename2="NGP.1.txt", temp=c(320,345,370,395,420,445,470,500,550,600),polymer="PMMA",method="all", unit=0.01)
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
## Warning: package 'ggplot2' was built under R version 3.3.2

alpha.ngp.plot(filename1="MSD.colmean.matrix.1.txt",filename2="NGP.1.txt", temp=c(320,345,370,395,420,445,470,500,550,600),polymer="PMMA",method="all", unit=0.0001)

PMMA g2 alpha term even orders (t^2 to t^20)

alpha.ngp.plot(filename1="MSD.colmean.matrix.1.txt",filename2="NGP.1.txt", temp=c(320,345,370,395,420,445,470,500,550,600),polymer="PMMA",method="even", unit=0.01)

alpha.ngp.plot(filename1="MSD.colmean.matrix.1.txt",filename2="NGP.1.txt", temp=c(320,345,370,395,420,445,470,500,550,600),polymer="PMMA",method="even", unit=0.0001)

### PMMA g2 alpha term even orders solomon method (t^2 to t^20)

alpha.ngp.plot(filename1="MSD.colmean.matrix.1.txt",filename2="NGP.1.txt", temp=c(320,345,370,395,420,445,470,500,550,600),polymer="PMMA",method="solomon",unit=0.01)

alpha.ngp.plot(filename1="MSD.colmean.matrix.1.txt",filename2="NGP.1.txt", temp=c(320,345,370,395,420,445,470,500,550,600),polymer="PMMA",method="solomon",unit=0.0001)

Here I used the time from 1:5001, which could be changed to 0:5000 with no change visually in shape of curve.