Exp Regression Slope Annualized with R-squared TRADINGVIEW

Vinks_o_7 · 2021-03-26T21:01:50+00:00

Hi I was wondering if it was possible to convert such indicator for PRT: this is a useful tool to build a portfolio by ranking stocks by strength based on Andreas Clenow’s book “Stocks on the Move” . Here’s the strategy: 1. Determine the 90 day annualized exponential regression slope for each stock 2. Determine the 90 day R-squared value for each of the stock 3. Multiply #1 * #2 and then sort from highest to lowest —> screener PRT 4. Buy the top stocks based on the formula ((AccountValue * 0.001)/ATR(20)) Thanks. //@version=4 // This source code is subject to the terms of the Mozilla Public License 2.0 at https://mozilla.org/MPL/2.0/ // Exponential Regression Slope Annualized with R-squared Histogram // by SparkyFlary study(title="Exponential Regression Slope Annualized with R-squared Histogram", shorttitle="Exponential Regression Slope Annualized with R-squared", format=format.percent, overlay=false) length = input(90, title="Slope length") rSquaredThreshold = input(0.0, title="R-squared threshold", minval=0, maxval=1) withRsquared = input(true, title="Multiply by R-squared?") withAnnualize = input(true, title="Annualize slope?") daysInYear = input(252, title="Number of bars in a year(if slope is annualized)") src = input(close, title="Source") //natural log of price nlPrice = log(src) //alexgrover's regression slope code x = bar_index y = nlPrice x_ = sma(x,length) y_ = sma(y,length) mx = stdev(x,length) my = stdev(y,length) c = correlation(x,y,length) slope = c * (my/mx) //inter = y_ - slope*x_ //reg = x*slope + inter exp_slope = slope annualized_slope = (pow(exp(slope),daysInYear) - 1) * 100 rSquared = pow(correlation(cum(1), nlPrice, length), 2) slopeAnnualized_with_Rsquared = withAnnualize ? annualized_slope * (withRsquared ? rSquared : 1) : exp_slope * (withRsquared ? rSquared : 1) trend = rSquared>rSquaredThreshold rising = slopeAnnualized_with_Rsquared>slopeAnnualized_with_Rsquared[1] falling = slopeAnnualized_with_Rsquared 0 down = slopeAnnualized_with_Rsquared<0 up_rising_color = color.new(color.green, 0) up_falling_color = color.new(color.green, 50) down_falling_color = color.new(color.red, 0) down_rising_color = color.new(color.red, 50) col = trend and up and rising?up_rising_color:trend and up and falling?up_falling_color:trend and down and falling?down_falling_color:trend and down and rising?down_rising_color : color.gray plot(slopeAnnualized_with_Rsquared, title="Annualized slope x R-squared", style=plot.style_histogram, color = col, linewidth=3)

Viewing 6 posts - 16 through 21 (of 21 total)

← 1 2

Author

Posts
09/17/2021 at 7:57 AM #177770 quote

RP-Markets
Participant

Junior

It seems that the condition that the price should be above the GD100 (pink line) does not work

EMA100.jpg

09/17/2021 at 8:52 AM #177782 quote

Nicolas
Keymaster

Legend

What is GD100 please? An MA with period 100?

09/17/2021 at 9:01 AM #177784 quote

RP-Markets
Participant

Junior

Oh sorry yes it was the German expression. 🙂

a simple moving average with 100 days period

09/17/2021 at 9:22 AM #177785 quote

jerome148
Participant

Junior

RP-Markets wrote:

Sorry now i see val = Momentum

but I actually forgot one point.

Can you still exclude stocks that have made a gap of more than 15% in the last few days?

Thank you very much for your work

This is what those 2 lines are for:

testVol = abs(close/close[1]–1) > 0.15

volatile = summation[90](testVol)

09/17/2021 at 9:36 AM #177788 quote

jerome148
Participant

Junior

RP-Markets wrote:

It seems that the condition that the price should be above the GD100 (pink line) does not work

Could you please check if you’re using a daily timeframe for your screener? There is a condition on the SMA100 so this is the only thing i can think of.

bullish = close > Average[100](Close)

09/17/2021 at 10:20 AM #177790 quote

RP-Markets
Participant

Junior

OK

i try it again and now i think it works

dont know what was the problem
Author

Posts