Upload New File

FractionModel.Rmd

+---
+title: "Fraction model"
+output: html_notebook
+---
+
+#Import data 
+```{r}
+library(readxl)
+FractionData <- read_excel("Fractiondata.xlsx", 
+    sheet = "Alldata_average")
+attach(FractionData) 
+    
+#rename data
+viscosity = log(`Final viscosity (mPa.s)`)
+dehulledcoarse = `Coarse fraction (%)` #dehulled from Nienke
+dehulledfine = `Fine fraction (%)` #dehulled from Nienke
+flour = `Flour (%)` #also dehulled 
+hull = `Hull (%)`
+protisol = `Protein isolate (%)`
+starchisol = `Starch isolate (%)`
+fibreisol = `Fibre isolate (%)`
+```
+
+#Load packages 
+```{r}
+library(MASS)
+library(tidyverse)
+library(ggplot2)
+library(heatmap3)
+library(corrplot)
+library(cowplot)
+library(car)
+library(R2HTML)
+library(readxl)
+```
+
+#Model with only parameters that make sense
+
+```{r}
+model= lm(viscosity~ dehulledcoarse + dehulledfine + flour + hull+ protisol + starchisol + fibreisol, data=FractionData)
+
+
+# model= lm(viscosity~coarse* fine* flour*hull* protisol* starchisol* fibreisol, data=FractionData)
+summary(model)
+```
+
+```{r}
+plot(model)
+
+```
+```{r}
+HTML2clip(summary(model))
+```