Hi,
I am beginning on R and I am struggleling a bit ^^.
It is my first post, so I may have forgot/don't undertstand some rules. I am sorry if it's the case.
If it's not in the right place, let me know it.
I try to build a little tool for data analysis and I am stuck. It is on a data set including (Name of the proteins/Normalized expression/What organism has been used/with what it has been infected/Organismpathogen***/number of replica/if it has been treat by placebo or real pathogen)
library(dplyr)
#>
#> Attachement du package : 'dplyr'
#> The following objects are masked from 'package:stats':
#>
#> filter, lag
#> The following objects are masked from 'package:base':
#>
#> intersect, setdiff, setequal, union
library(tidyverse)
#Building an AOV
data%>%aov(formula=formula("Norm~Strain*Cultivar*rep"))%>%tidy()%>%mutate(contribVariance=paste(round(100*meansq/sum(meansq),0),"%"))->aov
#> Error in tidy(.): impossible de trouver la fonction "tidy"
summary(aov)
#> Error in object[[i]]: objet de type 'closure' non indiçable
Created on 2021-03-03 by the reprex package (v1.0.0)
I don't understand why he doesn't find tidy (I took the formula from a lesson)
Furthermore I would like to comput a lot of Anova (I want to see wich protein have a significative differential of expression)
I was thniking something like :
library(dplyr)
#>
#> Attachement du package : 'dplyr'
#> The following objects are masked from 'package:stats':
#>
#> filter, lag
#> The following objects are masked from 'package:base':
#>
#> intersect, setdiff, setequal, union
library(tidyverse)
#Building an AOV for every protein
data%>%filter(Protein=="TraesCS1A01G002200.1")->data_premiere_prot
#> Error in UseMethod("filter"): pas de méthode pour 'filter' applicable pour un objet de classe "function"
data_premiere_prot%>%aov(formula=formula("Norm~Strain*Cultivar*rep"))->aovprot
#> Error in terms.formula(formula, "Error", data = data): objet 'data_premiere_prot' introuvable
summary(aovprot)
#> Error in summary(aovprot): objet 'aovprot' introuvable
Created on 2021-03-03 by the reprex package (v1.0.0)
How can I do for build a loop able to change automatically the name of my protein ?
ps : on my R I don't have any problem with filter, I don't know why the reprex found one ^^''.