Displaying 5 results from an estimated 5 matches for "regmod".
Did you mean:
regmon
2008 Dec 13
2
weird pasting of ".value" when list is returned
...code below is run and temp is returned.
I've been trying to figure this out for too long. It doesn't matter when
I put the FPVAL in the return statement. It happens regardless of
whether it's first or last. Thanks.
f.lmmultenhanced <-
function(response, pred1, pred2)
{
regmod <- lm(response ~ pred1 + pred2)
lmsum <- summary(regmod)
imbcoef<-lmsum$coefficients[2,1]
retcoef<-lmsum$coefficients[3,1]
imbpval<-lmsum$coefficients[2,4]
retpval<-lmsum$coefficients[3,4]
Fstat<-lmsum$fstatistic[1]
fpv...
2012 Mar 26
1
Seeming failure of options(width=60)
...reating a vignette. Is there something other than e.g., options(width=60) I
need to set? The Sweave FAQ suggests this should work.
options(width=60)
pastured <- data.frame(
time=c(9, 14, 21, 28, 42, 57, 63, 70, 79),
yield= c(8.93, 10.8, 18.59, 22.33, 39.35,
56.11, 61.73, 64.62, 67.08))
regmod<-"yield ~ t1 - t2*exp(-exp(t3+t4*log(time)))"
huetstart<-c(t1=70, t2=60, t3=0, t4=1)
anlsx<-try(nls(regmod, start=huetstart, trace=TRUE, data=pastured))
print(anlsx)
John Nash
2008 Dec 04
1
Formula parsing and updating
...;= 65" is not parsed as an entity.
It does not look as a problem with the environment to me, as is outlined
in http://developer.r-project.org/nonstandard-eval.pdf, rather a parsing
problem.
By contrast, the following works like a charm:
> x <- runif(50)
> y <- rnorm(50)+x
> regmod <- lm(y ~ x)
> update(regmod, " . ~ . + (ifelse(x>.5, 0, x))")
Could anybody give me a hint on my wrongdoings?
--
Johannes H?sing There is something fascinating about science.
One gets such wholesale returns of conjecture
mailto:jo...
2007 Dec 19
0
leaps
...1 ' ' 1
Residual standard error: 0.08135 on 68 degrees of freedom
Multiple R-Squared: 0.9932, Adjusted R-squared: 0.9914
F-statistic: 549.8 on 18 and 68 DF, p-value: < 2.2e-16*
newmat <- cosmat[,-c(5,8,10,12,13,15,16,18,20:44)]
newmat <- cbind(newmat,sinmat[,-c(7:11,13:44)])
Regmod <- regsubsets(newmat, yy)
rs <- summary(Regmod)
which.max(rs$adjr)
*[1] 8
*rs$which[which.max(rs$adjr), ]
*(Intercept) cos1 cos2 cos3 cos4 cos6
TRUE TRUE TRUE FALSE TRUE FALSE
cos7 cos9 cos11...
2009 Jul 12
2
Heckman Selection MOdel Help in R
...he
> above probit as the control in my outcome equation, ?hence I need to get the
> IMR (Is there another direct way?)
>
> 3. How can this be done in R using my concept or otherwise does there exist
> another way of doing what I wish to achieve
>
>
>
> On trying
>
> regmod <- heckit(s ~ age + gender + gemedu + gemhinc + es_gdppc +
>
> ??? imf_pop + estbbo_m, ln_oy5_1 ~ age+ gender+fearfail+gemedu,
> adpopdata,method="2step")
>
>
>
> I get
>
> Error: could not find function "heckit"
>
>
>
> Error: could not...