All functions
|
|
mf_bib()
|
Create a bib file for R packages, including the citations of user-defined packages. |
|
mf_convertc()
|
Convert air concentration between ppm and micromol m-3 |
|
mf_convertRH()
|
Convert the unit of air humidity according to Sonntag 1990. a in g/m-3 , RH in %, q in kg/kg. |
|
mf_daynight()
|
Define daytime or nighttime |
|
mf_dfplot()
|
Plot a dataframe, multiple ys against one x |
|
mf_dfplot2()
|
Plot a dataframe, one y against multiple xs |
|
mf_errorbar()
|
add error bars to a scatterplot. |
|
mf_fillna()
|
Fill a time series with NAs. in: a dataframe and a timestamp vector. out: a dataframe. |
|
mf_findpeaks()
|
Find peaks of a curve
#https://rtricks.wordpress.com/2009/05/03/an-algorithm-to-find-local-extrema-in-a-vector/ |
|
mf_findpeaks2()
|
Find peaks of a curve
http://stackoverflow.com/questions/31220307/calculate-x-value-of-curve-maximum-of-a-smooth-line-in-r-and-ggplot2 |
|
mf_hist()
|
Plot a user-customized hist |
|
mf_hmplot()
|
Plot a Hovmoeller or fingerprint plot |
|
mf_hourrose()
|
Hour rose plot |
|
mf_imagescale()
|
Create a legend |
|
mf_kurtosis()
|
Kurtosis measures whether the data are peaked or flat relative to a normal curve. positive: wide. negative: narrow |
|
mf_list2ascii()
|
Save a list into an ASCII file. in: a list. out: a file. |
|
mf_lm()
|
plot a linear regression figure and return a list of parameters. in: two vectors. out: a figure and a list. |
|
mf_lmdf()
|
calculate linear regression between every two columns in a data frame. in: a dataframes. out: a dataframe showing the linear regression. |
|
mf_names()
|
Enhancement of names() |
|
mf_optimdmodel()
|
optim function for the dmodel in chamber flux calculation. in: initial values for optim function. out: best fitted parameters |
|
mf_outlier()
|
check outliers. in: a vector. out: a plot and a vector of flags. the blank flag means passing all checks. |
|
mf_pairs()
|
plot pair-wise correlations. in: a dataframe. out: a figure. |
|
mf_pairs2()
|
plot pair-wise correlations with p value. in: a dataframe. out: a figure. |
|
mf_pairslm()
|
plot pair-wise correlations with linear regression. in: a dataframe. out: a figure. # not done yet. |
|
mf_pickquantile()
|
A filter to pick out data within a range by quantile of ref. in: a numeric vector. out: a logical vector. |
|
mf_planarfit()
|
calculate the planafit coefficients from 3D wind measurement. in: three vectors or one-column dataframes. out: a list. |
|
mf_plotblank()
|
plot a blank figure |
|
mf_plotcolors()
|
A reminder for colors |
|
mf_plotlty()
|
lty |
|
mf_plotpch()
|
pch numbers |
|
mf_plottype()
|
type |
|
mf_prefix0()
|
fill 0 before a string. in and out: a vector with length 1. if length(x) > 1, use unlist(lapply(x,FUN = mf_fill0)). e.g. 12 --- 012 |
|
mf_rainbow()
|
# give each x a color by groups. |
|
mf_readdir()
|
batch read several files. in: a path. out: a list containing all data frames |
|
mf_satpress()
|
calculation of saturation vapour pressure (Pa) at T(degree C). in and out: a vector. for -45 to 60 degree over water according to Sonntag 1990 |
|
mf_se()
|
standard error |
|
mf_sharedata()
|
A template to create a folder with data files to share |
|
mf_skewness()
|
normality test. if either skewness/se_skew is outside -1.96 -- 1.96, the data are unlikely to be normally distributed. Or Kolmogorov-Smirnov test, Shapiro-Wilks' W test. But a visual examination is the best. Negative values of the skewness indicate data that are skewed to the left(negativelz skewed) |
|
mf_smooth()
|
smooth a series with a giving width. in and out: a vector. |
|
mf_strptime()
|
A simplied version of strptime. only for '%Y-%m-%d %H:%M:%S' |
|
mf_sunriset()
|
calculate sunrise and sunset time in a friendly way. in: a given date and coordinates. out: a dataframe with sunrise and sunset time. |
|
mf_tapply()
|
a friendly version of tapply for dataframes. in and out: same as tapply(). the built-in function tapply returns a matrix with unfriendly row name and colname. mf-tapply returns a friendly dataframe |
|
mf_taylor()
|
mf_taylor: plot a taylor diagram to compare reference (x) and model (y). in: two vectors. out: a figure. |
|
mf_timefiller()
|
Fill time series with NA |
|
mf_timestampfull()
|
convert time stamps into a full format, i.e.
convert 2016-07-14 or 2016-07-14 00:00 into 2016-07-14 00:00:00
convert 14.07.2016 or 14.07.2016 00:00 into 14.07.2016 00:00:00
convert 14.07.2016 or 14.07.2016 00:00 into 14.07.2016 00:00:00
convert 20160714 or 2016071400 or 201607140000 into 20160714000000 |
|
mf_vtapply()
|
a friendly version of tapply |
|
mf_windd()
|
calculate resultant wind direction from u and v. in: a vector of u and v. out: a vector of wind direction. |
|
mf_winddclear()
|
convert wind direction out of the range [0, 360) into the range. |
|
mf_windmean()
|
calculate resultant mean wind spead and resultant mean wind direction.# in: a wind speed vector and a wind direction vector. # out: a mean speed and a mean direction. |
|
mf_windmean2()
|
not recommended! alculate resultant mean wind speed and resultant mean wind direction. especially used in tapply(). this calc is sometimes confusing. better to calc step by step. in: a character vector with wind speed and direction separated with semi colon ';'. out: a chracter vector with resultant mean wind speed and resultant mean direction separated with semi colon ';'. |
|
mf_windrose()
|
Draw my windrose |
|
mf_windu()
|
calculate u component of wind. in: a vector of wind speed and a vector of wind direction. out: a vector of u |
|
mf_windv()
|
calculate v component of wind. in: a vector of wind speed and a vector of wind direction. out: a vector of v |
|
mf_write()
|
save csv file with asking if the file already exists. |