Skip to content

Commit

Permalink
No commit message
Browse files Browse the repository at this point in the history
  • Loading branch information
@PasaOpasen
PasaOpasen committed Feb 13, 2020
1 parent d413de7 commit 9c5cec7
Show file tree
Hide file tree
Showing 6 changed files with 225 additions and 179 deletions.
Binary file renamed ... по пауэрлифтингу и не только/циклы/9.xls → ...лько/циклы/Пиковый цикл Верхошанского.xls
View file
Binary file not shown.
Binary file renamed ...по пауэрлифтингу и не только/циклы/26.xls → ... не только/циклы/гусеница Суровецкого.xls
View file
Binary file not shown.
336 changes: 168 additions & 168 deletions Уравнение для повторного максимума/.Rhistory
View file
Original file line number Diff line number Diff line change
@@ -1,171 +1,3 @@
Show(predict(md,data%>% select(Val,Count,Weight,High,Body)),allrows)
md=lm(SM~Val:Count+Val*Body-1+I(Val*Weight/(High-100))+I((Val-100)/Val)),data)
summary(md)
Show(predict(md,data%>% select(Val,Count,Weight,High,Body)),allrows)
md=lm(SM~Val:Count+Val*Body-1+I(Val*Weight/(High-100))+I((Val-100)/Val),data)
summary(md)
Show(predict(md,data%>% select(Val,Count,Weight,High,Body)),allrows)
#Загрузка данных####
library(tidyverse)
library(magrittr)
data=read_tsv("data.tsv",
skip=1,col_names = F,na="",
col_types = "fddnfffnnnff"
) %>% tbl_df()
colnames(data)=c("Date","SM","Val","Count","Type","Sex","Experience","Age","Weight","High","Body","Mail")
data %<>% mutate(CountGroup=cut(Count,breaks = c(1,3,6,9,12,20,30)))
allrows=1:nrow(data)
#Разведочный анализ и описание выборки####
summary(data)
Error=function(target,weight) (target-weight)^2 %>% sum()
Show=function(vals,rows){
cbind(value=vals,Target=data$SM[rows],
ERROR=abs(data$SM[rows]-vals),ERRORproc=abs(data$SM[rows]-vals)/data$SM[rows]*100,
data[rows,c(3:11,13)]) %>% tbl_df() %>% print()
cat("---------------------> Error is",Error(vals,data$SM[rows]),"\n")
}
Show(data$Val*(1+0.0333*data$Count),allrows)
Error=function(target,weight) (target-weight)^2 %>% sum()
Show=function(vals,rows){
cbind(value=vals,Target=data$SM[rows],
ERROR=abs(data$SM[rows]-vals),ErrorPercent=abs(data$SM[rows]-vals)/data$SM[rows]*100,
data[rows,c(3:11,13)]) %>% tbl_df() %>% print()
cat("---------------------> Error is",Error(vals,data$SM[rows]),"\n")
}
Show(data$Val*(1+0.0333*data$Count),allrows)
md=lm(SM~Val:Count+Val*Body-1+I(Val*Weight/(High-100))+I((Val-100)/Val),data)
summary(md)
Show(predict(md,data%>% select(Val,Count,Weight,High,Body)),allrows)
#Загрузка данных####
library(tidyverse)
library(magrittr)
data=read_tsv("data.tsv",
skip=1,col_names = F,na="",
col_types = "fddnfffnnnff"
) %>% tbl_df()
colnames(data)=c("Date","SM","Val","Count","Type","Sex","Experience","Age","Weight","High","Body","Mail")
data %<>% mutate(CountGroup=cut(Count,breaks = c(1,3,6,9,12,20,30)),ValCoef=Val/Weight,SMcoef=SM/Weight)
allrows=1:nrow(data)
data
md=lm(SM~Val:Count+Val-1,data)
summary(md)
Show(predict(md,data %>% select(Val,Count)),allrows)
md=lm(SMcoef~ValCoef:Count+ValCoef-1,data)
summary(md)
Show(predict(md,data %>% select(ValCoef,Count))*data$Weight,allrows)
md=lm(SM~Val:Count+Val*Body-1+I(Val*Weight/(High-100))+I((Val-100)/Val),data)
summary(md)
Show(predict(md,data%>% select(Val,Count,Weight,High,Body)),allrows)
pr=predict(md,data[c(3:4,9,13)], interval = "prediction", level = 0.95)
obj+
#geom_ribbon(aes(x=Val,ymin = pr[,2], ymax = pr[,3]), fill = "grey70") +
geom_line(aes(x=Val,y=pr[,1]),size=1)+
geom_point(aes(x=Val,y=pr[,1]),size=3)+
geom_point(aes(x=Val,y=SM,col=Body,shape=Type),size=4)+theme(legend.position = c(0.85,0.3))
obj=ggplot(data %>% select(-Date,-Mail))+theme_bw()
pr=predict(md,data[c(3:4,9,13)], interval = "prediction", level = 0.95)
obj+
#geom_ribbon(aes(x=Val,ymin = pr[,2], ymax = pr[,3]), fill = "grey70") +
geom_line(aes(x=Val,y=pr[,1]),size=1)+
geom_point(aes(x=Val,y=pr[,1]),size=3)+
geom_point(aes(x=Val,y=SM,col=Body,shape=Type),size=4)+theme(legend.position = c(0.85,0.3))
pr=predict(md,data%>% select(Val,Count,Weight,High,Body), interval = "prediction", level = 0.95)
obj+
#geom_ribbon(aes(x=Val,ymin = pr[,2], ymax = pr[,3]), fill = "grey70") +
geom_line(aes(x=Val,y=pr[,1]),size=1)+
geom_point(aes(x=Val,y=pr[,1]),size=3)+
geom_point(aes(x=Val,y=SM,col=Body,shape=Type),size=4)+theme(legend.position = c(0.85,0.3))
getPIE=function(vec,main=""){
ln=length(levels(vec))
x=numeric(ln)
ns=character(ln)
for(i in seq(ln)){
x[i]=sum(vec==levels(vec)[i])/length(vec)
ns[i]=paste0(levels(vec)[i]," (",round(x[i]*100,2),"%)")
}
pie(x=x,labels=ns,main=main)
}
getPIE(data$Body)
getPIE(data$CountGroup)
getPIE(data$Experience)
getPIE(factor(data$Count))
obj=ggplot(data %>% select(-Date,-Mail))+theme_bw()
obj+geom_bar(aes(x=CountGroup))
obj+geom_bar(aes(x=Body))
obj+geom_boxplot(aes(x=Type,y=SM))
GGally::ggpairs(data %>% select(-Date,-Mail))
#Загрузка данных####
library(tidyverse)
library(magrittr)
data=read_tsv("data.tsv",
skip=1,col_names = F,na="",
col_types = "fddnfffnnnff"
) %>% tbl_df()
colnames(data)=c("Date","SM","Val","Count","Type","Sex","Experience","Age","Weight","High","Body","Mail")
data %<>% mutate(CountGroup=cut(Count,breaks = c(1,3,6,9,12,20,30)),ValCoef=Val/Weight,SMcoef=SM/Weight)
allrows=1:nrow(data)
obj=ggplot(data %>% select(-Date,-Mail))+theme_bw()
md=lm(SM~Val:Count+Val*Body-1+I(Val*Weight/(High-100))+I((Val-100)/Val),data)
summary(md)
Show(predict(md,data%>% select(Val,Count,Weight,High,Body)),allrows)
Error=function(target,weight) (target-weight)^2 %>% sum()
Show=function(vals,rows){
cbind(value=vals,Target=data$SM[rows],
ERROR=abs(data$SM[rows]-vals),ErrorPercent=abs(data$SM[rows]-vals)/data$SM[rows]*100,
data[rows,c(3:11,13)]) %>% tbl_df() %>% print()
cat("---------------------> Error is",Error(vals,data$SM[rows]),"\n")
}
md=lm(SM~Val:Count+Val*Body-1+I(Val*Weight/(High-100))+I((Val-100)/Val),data)
summary(md)
Show(predict(md,data%>% select(Val,Count,Weight,High,Body)),allrows)
p=prcomp(~Val+Count+High+Weight,data)
summary(p)
pr=predict(md,data%>% select(Val,Count,Weight,High,Body), interval = "prediction", level = 0.95)
obj+
#geom_ribbon(aes(x=Val,ymin = pr[,2], ymax = pr[,3]), fill = "grey70") +
geom_line(aes(x=Val,y=pr[,1]),size=1)+
geom_point(aes(x=Val,y=pr[,1]),size=3)+
geom_point(aes(x=Val,y=SM,col=Body,shape=Type),size=4)+theme(legend.position = c(0.85,0.3))
install.packages("ggformula")
library(ggformula)
pr=predict(md,data%>% select(Val,Count,Weight,High,Body), interval = "prediction", level = 0.95)
obj+
#geom_ribbon(aes(x=Val,ymin = pr[,2], ymax = pr[,3]), fill = "grey70") +
geom_spline(aes(x=Val,y=pr[,1]),size=1)+
geom_point(aes(x=Val,y=pr[,1]),size=3)+
geom_point(aes(x=Val,y=SM,col=Body,shape=Type),size=4)+theme(legend.position = c(0.85,0.3))
pr=predict(md,data%>% select(Val,Count,Weight,High,Body), interval = "prediction", level = 0.95)
obj+
#geom_ribbon(aes(x=Val,ymin = pr[,2], ymax = pr[,3]), fill = "grey70") +
geom_spline(aes(x=Val,y=pr[,1]),size=1,spar=0.2)+
geom_point(aes(x=Val,y=pr[,1]),size=3)+
geom_point(aes(x=Val,y=SM,col=Body,shape=Type),size=4)+theme(legend.position = c(0.85,0.3))
pr=predict(md,data%>% select(Val,Count,Weight,High,Body), interval = "prediction", level = 0.95)
obj+
#geom_ribbon(aes(x=Val,ymin = pr[,2], ymax = pr[,3]), fill = "grey70") +
geom_spline(aes(x=Val,y=pr[,1]),size=1,spar=0.4)+
geom_point(aes(x=Val,y=pr[,1]),size=3)+
geom_point(aes(x=Val,y=SM,col=Body,shape=Type),size=4)+theme(legend.position = c(0.85,0.3))
pr=predict(md,data%>% select(Val,Count,Weight,High,Body), interval = "prediction", level = 0.95)
obj+
#geom_ribbon(aes(x=Val,ymin = pr[,2], ymax = pr[,3]), fill = "grey70") +
geom_spline(aes(x=Val,y=pr[,1]),size=1,spar=0.02)+
geom_point(aes(x=Val,y=pr[,1]),size=3)+
geom_point(aes(x=Val,y=SM,col=Body,shape=Type),size=4)+theme(legend.position = c(0.85,0.3))
pr=predict(md,data%>% select(Val,Count,Weight,High,Body), interval = "prediction", level = 0.95)
obj+
#geom_ribbon(aes(x=Val,ymin = pr[,2], ymax = pr[,3]), fill = "grey70") +
geom_spline(aes(x=Val,y=pr[,1]),size=1,spar=0.02,cv=T)+
geom_point(aes(x=Val,y=pr[,1]),size=3)+
geom_point(aes(x=Val,y=SM,col=Body,shape=Type),size=4)+theme(legend.position = c(0.85,0.3))
pr=predict(md,data%>% select(Val,Count,Weight,High,Body), interval = "prediction", level = 0.95)
obj+
#geom_ribbon(aes(x=Val,ymin = pr[,2], ymax = pr[,3]), fill = "grey70") +
geom_spline(aes(x=Val,y=pr[,1]),size=1,spar=0.02,cv=F)+
geom_point(aes(x=Val,y=pr[,1]),size=3)+
geom_point(aes(x=Val,y=SM,col=Body,shape=Type),size=4)+theme(legend.position = c(0.85,0.3))
pr=predict(md,data%>% select(Val,Count,Weight,High,Body), interval = "prediction", level = 0.95)
obj+
#geom_ribbon(aes(x=Val,ymin = pr[,2], ymax = pr[,3]), fill = "grey70") +
geom_line(aes(x=Val,y=pr[,1]),size=1,col="grey")+
geom_point(aes(x=Val,y=pr[,1]),size=3)+
Expand Down Expand Up @@ -510,3 +342,171 @@ tb=tibble(d=1:length(vc),val=vc) %>% mutate(day=factor(ifelse(d%%(length(x)+1)==
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1)+
geom_point(aes(col=day),size=4)
library(tidyverse)
library(magrittr)
library(ggformula)
library(ggalt)
#простой цикл####
mx=function(vals,count){
t=vals*(1+0.0333*count)
t[count==1]=vals[count==1]
s=t[length(t)]*(1+runif(1,0.,0.02))
return(c(t,s))
}
x=c(100,120,115,105,130,120,125)
y=c(10,4,6,10,3,6,5)
vc=c(mx(x,y),mx(x+4,y),mx(x+9,y))
tb=tibble(d=1:length(vc),val=vc) %>% mutate(day=factor(ifelse(d%%(length(x)+1)==0,"проходка","тренировка")))
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1)+
geom_point(aes(col=day),size=4)
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1)+
geom_point(aes(col=day),size=4)+theme_bw()
mx=function(vals,count){
t=vals*(1+0.0333*count)
t[count==1]=vals[count==1]
s=t[length(t)]*(1+runif(1,0.,0.01))
return(c(t,s))
}
x=c(100,120,115,105,130,120,125)
y=c(10,4,6,10,3,6,5)
vc=c(mx(x,y),mx(x+4,y),mx(x+9,y))
tb=tibble(d=1:length(vc),val=vc) %>% mutate(day=factor(ifelse(d%%(length(x)+1)==0,"проходка","тренировка")))
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1)+
geom_point(aes(col=day),size=4)+theme_bw()
#простой цикл####
mx=function(vals,count){
t=vals*(1+0.0333*count)
t[count==1]=vals[count==1]
s=t[length(t)]*(1+runif(1,0.,0.005))
return(c(t,s))
}
x=c(100,120,115,105,130,120,125)
y=c(10,4,6,10,3,6,5)
vc=c(mx(x,y),mx(x+4,y),mx(x+9,y))
tb=tibble(d=1:length(vc),val=vc) %>% mutate(day=factor(ifelse(d%%(length(x)+1)==0,"проходка","тренировка")))
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1)+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1,spline_shape = 1)+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1,spline_shape = 0.2)+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1,spline_shape = 0.4,line_type="dash")+
geom_point(aes(col=day),size=4)+theme_bw()
vignette("ggplot2-specs")
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1,spline_shape = 0.6,linetype="dashed")+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1,spline_shape = 0.5,linetype="dashed")+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1,spline_shape = 0.7,linetype="dashed")+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1,spline_shape = 0.8,linetype="dashed")+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1,spline_shape = -0.8,linetype="dashed")+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1,spline_shape = -0.5,linetype="dashed")+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1,spline_shape = -0.3,linetype="dashed")+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1,spline_shape = -0.3,linetype="dotdash",alpha=0.5)+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1,spline_shape = -0.3,linetype="dotdash",alpha=0.7)+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1,spline_shape = -0.3,linetype="dotdash",col="orange",alpha=0.7)+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1,spline_shape = -0.3,linetype="dotdash",col="gold",alpha=0.7)+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1,spline_shape = -0.3,linetype="dotdash",col="orange",alpha=0.7)+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1,spline_shape = -0.3,linetype="dotdash",col="orange",alpha=0.8)+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1,spline_shape = -0.3,linetype="dotdash",col="orange",alpha=0.9)+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1,spline_shape = -0.3,linetype="dotdash",col="red",alpha=0.9)+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_xspline(size=1,spline_shape = -0.3,linetype="dotdash",col="green",alpha=0.9)+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_hline(yintercept = 140)+
geom_xspline(size=1,spline_shape = -0.3,linetype="dotdash",col="green",alpha=0.9)+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_hline(yintercept = 142)+
geom_hline(yintercept = vc[length(vc)])+
geom_xspline(size=1,spline_shape = -0.3,linetype="dotdash",col="green",alpha=0.9)+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_hline(yintercept = 142,show.legend = T)+
geom_hline(yintercept = vc[length(vc)])+
geom_xspline(size=1,spline_shape = -0.3,linetype="dotdash",col="green",alpha=0.9)+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_hline(yintercept = 142,size=1.2)+
geom_hline(yintercept = vc[length(vc)])+
geom_xspline(size=1,spline_shape = -0.3,linetype="dotdash",col="green",alpha=0.9)+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_hline(yintercept = 142,size=1.2,linetype="dashed")+
geom_hline(yintercept = vc[length(vc)])+
geom_xspline(size=1,spline_shape = -0.3,linetype="dotdash",col="green",alpha=0.9)+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_hline(yintercept = 142,size=1.2,linetype="dashed")+
geom_hline(yintercept = vc[length(vc)],size=1.2,linetype="dashed")+
geom_xspline(size=1,spline_shape = -0.3,linetype="dotdash",col="green",alpha=0.9)+
geom_point(aes(col=day),size=4)+theme_bw()
ggplot(tb,aes(x=d,y=val))+
geom_hline(yintercept = 142,size=1.1,linetype="dashed")+
geom_hline(yintercept = vc[length(vc)],size=1.1,linetype="dashed")+
geom_xspline(size=1,spline_shape = -0.3,linetype="dotdash",col="green",alpha=0.9)+
geom_point(aes(col=day),size=4)+theme_bw()
mx2=function(vals,count){
t=vals*(1+0.0333*count)
t[count==1]=vals[count==1]
return(t)
}
#155
x=c(100,125,107.5,132.5,120,140,127.5,147.5,117.5,132.5,125,162.5)
y=c(6,5,5,4,3,3,3,2,5,2,2,1)
vc=c(mx2(x,y),mx2(x+7.5,y),mx2(x+15,y))
plot(vc)
plot(vc,type = "b")
vc=c(mx2(x,y),mx2(x+8,y),mx2(x+15,y))
plot(vc,type = "b")
x=c(112,112,117,122,126,131,136,141,145,150,155,160,166,177,122.5,122.5,128,133,138,143,149,154,159,165,170,175,182,195,136,136,142,148,154,160,165,171,177,183,190,195,203,216)
y=rep(c(10,10,8,8,5,5,5,5,3,3,2,2,1,1),3)
vc=c(mx2(x,y))
plot(vc,type = "b")
mx=function(vals,count){
t=vals*(1+0.0333*count)
t[count==1]=vals[count==1]
s=t[length(t)]*(1+runif(1,0.,0.005))
return(c(t,s))
}
ggplot(tb,aes(x=d,y=val))+
geom_hline(yintercept = 141,size=1.1,linetype="dashed")+
geom_hline(yintercept = vc[length(vc)],size=1.1,linetype="dashed")+
geom_xspline(size=1,spline_shape = -0.3,linetype="dotdash",col="green",alpha=0.9)+
geom_point(aes(col=day),size=4)+theme_bw()
Loading

0 comments on commit 9c5cec7

Please sign in to comment.