search for: rsquare

Displaying 20 results from an estimated 87 matches for "rsquare".

Did you mean: square
2013 May 16
2
R looping help
...having trouble with my R code. I am trying to get the p-values, R^2s etc for a number of different groups of variables that are all in one dataset. This is the code: #Stand counter st<-1 #Collections stands<-numeric(67) slopes<-numeric(67) intercepts<-numeric(67) mses<-numeric(67) rsquares<-numeric(67) pValues<-numeric(67) #Start lists for X and Y values within each stand xi<-numeric(0) yi<-numeric(0) #Set the first element to the starting X and Y values xi[1]=X[1] yi[1]=Y[1] #Start looping working through your data, record by record for (i in 2:length(X)) { #If you ar...
2007 Sep 13
5
statistics - hypothesis testing question
I estimate two competing simple regression models, A and B where the LHS is the same in both cases but the predictor is different ( I handle the intercept issue based on other postings I have seen ). I estimate the two models on a weekly basis over 24 weeks. So, I end up with 24 RSquaredAs and 24 RsquaredBs, so essentally 2 time series of Rsquareds. This doesn't have to be necessarily thought of as a time series problem but, is there a usual way, given the Rsquared data, to test H0 : Rsquared B = Rsquared A versus H1 : Rsquared B > Rsquared A so that I can map the 24 R...
2005 Feb 08
3
logistic regression
Hi, I'm using glm function to do logistic regression and now I want to know if it exists a kind of R-squared with this function in order to check the model. Thank you.
2003 Jul 21
1
how to calculate Rsquare
...y Df Sum Sq Mean Sq F value Pr(>F) x 1 4062.9 4062.9 400.04 2.990e-13 *** f 1 4421.5 4421.5 435.35 1.496e-13 *** Residuals 17 172.7 10.2 --- Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1 > #rsquare of model > (4062.9+4421.5)/(4062.9+4421.5+172.7) [1] 0.980051 In this way I calculate the model rsquare, but how to calculate the rsquare of each levels "a" and "b"? This is only an example, the model maybe glm, lme etc. Thanks Ronaldo -- NOWPRINT. NOWPRINT. Clemclone,...
2011 Jul 13
1
max possible rsquare
...- Signif. codes: 0 ?***? 0.001 ?**? 0.01 ?*? 0.05 ?.? 0.1 ? ? 1 exp(coef) exp(-coef) lower .95 upper .95 smocc_zyban 0.6451 1.5502 0.1008 4.130 self 3.2728 0.3055 2.4666 4.343 smocc_zyban:self 1.1529 0.8673 0.7229 1.839 Rsquare= 0.116 (max possible= 0.821 ) Likelihood ratio test= 137 on 3 df, p=0 Wald test = 114.3 on 3 df, p=0 Score (logrank) test = 141.3 on 3 df, p=0, Robust = 70.98 p=2.665e-15 -- View this message in context: http://r.789695.n4.nabble.com/max-possible-rsquare-tp3664925p3664925....
2011 Apr 15
3
Rsquared for anova
...#39;,'tt1','tt2','tt3','path','iters','type','length','tos','tws','time') fit<-(aov(tos~nh1*nh2*nh3*randsize*aweights*tt1*tt2*tt3*iters*length, data=expdata)) summary(fit) I want to check the fit of the model with Rsquared, but I cannot seem to find a way to do this... Any help is much appreciated... Thanks, Dorien [[alternative HTML version deleted]]
2004 Nov 20
1
annotation problems (conditional text())
...dge into the right direction? Thanks, Joh if(m >= 0){ text( 4,0.19, cex=0.75, adj=0, paste( substitute( y == m*x+b, list( m=round(m,digits=5), b=round(b,digits=5) ) ), "\n", substitute( R^2==rsquared, list( rsquared=round(summary(fit)$r.squared,digits=3) ) ) ) ); text( 4,0.19, cex=0.75, adj=0, paste( substitute( y == m*x-b, list( m=round(m,digits=5), b=round(sqrt(b^2),digits=5) ) ),...
2012 Jan 17
2
result numeric(0) when using variable1[which(variable2="max(variable2)"]
Dear all, I have a question about the knowing for which row I have the max value of one of my variables. I calculated the Rsquared for different columns and made a list to gather them. I unlisted this list to create a vector with this values. I want to know for which column I have the max value of Rsquared. The columns were always named in the same way. They always start with results4$depth_ following by the number. The num...
2005 Jun 22
0
Rsquare from glmmPQL or another GLMM?
Hi, I know that Rsquare in glm or in non-linear models is "wrong", but some people like this. How I make to estimate the Rsquare from a model ajusted with glmmPQL or another GLMM? Thanks for all Ronaldo -- A simplicidade ?? o ??ltimo degrau da sabedoria. -- Victor Hugo -- |> // | \\ [...
2003 Aug 27
1
how to calculate Rsquare
...ntation of R, comparison and compatibility with S and S-plus, and for the posting of nice examples and benchmarks." -----Original Message----- From: Ronaldo Reis Jr. [mailto:chrysopa at insecta.ufv.br] Sent: Wednesday, August 27, 2003 10:04 AM To: R-Help Subject: Re: [R] how to calculate Rsquare Can anybody send these articles for me? > NagelKerke, N. J. D. (1991) "A note on a general definition of the > coefficient of determination", Biometrika 78: 691-2. > > Cox, D. R. and Wermuth, N. (1992) "A comment on the coefficient of > determination for bi...
2011 Jul 12
0
Adding a correlation value (like Rsquared) to a 4 parameter logistic fit model.
...d, data = SC, fct = LL.4()) plot(FourP, main = "LTB4 Standard Curve Zi Phase 7", xlab = "Expected (pg/mL)", ylab = "Response (%Bound)") Thanks for any help. Kevin McEnroy -- View this message in context: http://r.789695.n4.nabble.com/Adding-a-correlation-value-like-Rsquared-to-a-4-parameter-logistic-fit-model-tp3662480p3662480.html Sent from the R help mailing list archive at Nabble.com.
2012 Feb 16
1
how to get r-squared for a predefined curve or function with "other" data points
...{ x <- subset(curvedata,Group == i)$MPa y <- subset(curvedata,Group == i)$PLC name <- subset(curvedata,Group == i)$Sample points(x,y) vlc <- nls(y ~ 100/(1+exp(a*(x-b))), start=c(a=1, b=-3), data=list(x,y)) curve(100/(1+exp(coef(vlc)[1]*(x-coef(vlc)[2]))), col=1, add = TRUE) Rsquared <- 1 - var(residuals(vlc))/var(y) summarizeall[i ,"Run"] <- i summarizeall[i ,"Sample"] <- name[1] summarizeall[i ,"a"] <- coef(vlc)[1] summarizeall[i ,"b"] <- coef(vlc)[2] summarizeall[i ,"R2"] <- Rsquared listnow <- d...
2002 May 09
4
Rsquared in summary(lm)
Hello, I'm doing some linear regression: >lm<-lm(osas~alp,data) >summary(lm) However, the Rsquared in the output of summary() is not the same as the "standard" Rsquared calculated by spreadsheets, and outlined in statistical guidebooks, being SSR/SSTO. The output says "multiple Rsquared", but it is no multiple regression... What's the difference? Thanks, Wouter Buytae...
2007 Dec 08
1
lm: how to calculate rsquared of the predicted values?
Hi, I've built a linear model using multiple linear regression which leads me a R-squared value of 73.58%. After that, I used this model to predicted some values based on the test data. Now I'm wondering how: 1. can I measure de R-squared value between the predicted(by the model) and real (observed) values.? 2. Measure the RMSE error . Example: suppose my data its below: REAL
2013 Jun 18
2
offset en bucle
...Summ <- summary(GAM) Var.p.values <- Summ$s.table[,"p-value"] MaxSignif <-c(MaxSignif,max(Var.p.values)) Significance <- all(Var.p.values<0.1) ### Especificar el p values que queremos dentro del modelo GlobalSignificance <- c(GlobalSignificance,Significance) RSquare <- c(RSquare,Summ$r.sq) } Si teneis alguna sugerencia os estaría muy agradecido -- *Juan Diego Alcaraz Hernández * ** *http://juandiegoalcaraz.wordpress.com/* * Institut d'Investigació per a la Gestió Integrada de Zones Costaneres (IGIC) Universitat Politècnica de València C/ Paranimf, 1...
2001 Apr 19
0
All subsets weighted regression in R ?
To everyone, I am currently doing analysis through Splus but am having severe problems. I am trying to fit all subset weighted regression to a multivariate data set and get PRESS statistics, along with Rsquare, Rsquare adjusted and Mallows Cp for all possible models. Splus, will give me these statistics but the leaps procedure is highly sensitive to nearly singular data matrices (which mine is) and the Rsquare and adjusted Rsquare values it gives me are biased high (seems to be the case whenever you use...
2012 Aug 06
1
cannot find function "simpleRDA2"
...s function implements FS using the F-test described ## by Miller and Farr (1971). This test requires that ## -- the Y variables be standardized, ## -- the error in the response variables be normally distributed (to be verified by the user). ## ## This function uses 'simpleRDA2' and 'RsquareAdj' developed for 'varpart' in 'vegan'. ## ## Pierre Legendre & Guillaume Blanchet, May 2007 ## ## Arguments -- ## ## Y Response data matrix with n rows and m columns containing quantitative variables. ## X Explanatory data matrix with n rows a...
2007 Aug 28
1
FW: How to fit an linear model withou intercept
...It's actually in a time series framework also but that's not relevant > for my question. My question has 2 parts : > > 1) I was leaning towards using the R squared as the decision criteria ( > I will be Regressing monthly and over a couple of years so I will have > about 24 rsquareds. I have tons of data For one monthly regression so I > don't have to just do one big regression over the whole time period ) > but I noticed in your previous example that the model with intercept ( > compared to the model forced to have zero intercept ) had a lower R^2 > and a low...
2007 Dec 11
1
postResample R² and lm() R²
Hello, I'm with a conceptual doubt regarding Rsquared of both lm() and postResample(library caret). I've got a multiple regression linear model (lets say mlr) with anR² value of 67.52%. Then I use this model pro make predictions with predict() function using the same data as input , that is, use the generated model to predict the value associate...
2013 Mar 02
2
caret pls model statistics
...ies) inTrain1=createDataPartition(y=iris$Species, p=.75, list=FALSE) training1=iris[inTrain1,] testing1=iris[-inTrain1,] ctrl1=trainControl(method="cv", number=10) plsFit2=train(Species~., data=training1, method="pls", trControl=ctrl1, metric="Rsquared", preProc=c("scale")) data(iris) training1=iris[inTrain1,] datvars=training1[,1:4] dat.sc=scale(datvars) n=nrow(dat.sc) dat.indices=seq(1,n) timematrix=with(training1, classvec2classmat(Species[dat.indices])) pls.dat=plsr(timematrix ~ dat.sc, ncomp=3, method=&quo...