Displaying 5 results from an estimated 5 matches for "phihat".
2002 Jul 03
1
complex numbers in library function (PR#1737)
Full_Name: Alex Deckmyn
Version: 1.5.1
OS: Linux Mandrake 7.0
Submission from: (NULL) (193.190.63.62)
A short function in a library contains the complex number -2i :
phihat <- function(x)
{
x+ (-2i)
}
If i read the source file with source("phihat") this works fine. If I include it
as
a function in a library and read it with library(test), the function is still
read corectly, but when I print it out it shows as
> phihat
function (x)
{
x + (...
2012 Feb 10
3
Help needed please
...4
v1 <- cbind(simtimeseries[1:1020])
v2 <- cbind(simtimeseries[5:1024])
v3 <- t(v1)%*%v2
s4 <- v3/1204
i0 <- c(s0,s1,s2,s3)
i1 <- c(s1,s0,s1,s2)
i2 <- c(s2,s1,s0,s1)
i3 <- c(s3,s2,s1,s0)
gamma <- cbind(i0,i1,i2,i3)
eta <-c(s1,s2,s3,s4)
inversegamma <- solve(gamma)
phihat <- inversegamma%*%eta
phihat
Phihat <- cbind(phihat)
s <- c(s1,s2,s3,s4)
S <- cbind(s)
sigmasquaredyule <- s0 - (t(Phihat)%*%S)
sigmasquaredyule
I did a yule walker estimate on the simulated data and wanted to work out phi hat which is a vector of 4 values and sigmasquaredyule wh...
2009 Feb 13
1
need help with errors in betareg analysis
...g to fit a model in betareg and I'm getting errors, but have no
idea what they mean or how to solve them. Does anyone have experience with
this?
> model <- betareg(ACT ~ ST*SoilT, data = actDL_F)
Warning messages:
1: In sqrt(W) : NaNs produced
2: In sqrt(W) : NaNs produced
3: In sqrt(1 + phihat) : NaNs produced
data summaries don't give any na's or problems I can see
> summary(actDL_F$ACT)
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.0006227 0.0505600 0.2163000 0.2643000 0.4084000 0.9994000
> summary(actDL_F$ST)
Min. 1st Qu. Median Mean 3rd Qu....
2008 Oct 21
2
Question about glm using R
Good morning,
I am using R to try to model the proportion of burned area in Portugal.
The dependent variable is the proportion. The family used is binomial
and the epsilon would be binary.
I am not able to find the package to be used when the proportion (%) has
to be used in glm. Could someone help me?
I am using normal commands of glm.. for example:
glm_5<- glm(formula=p~Precipitation,
2011 Mar 12
3
betareg help
Dear R users,
I'm trying to do betareg on my dataset.
Dependent variable is not normally distributed and is proportion (of condom
use (0,1)).
But I'm having problems:
gyl<-betareg(cond ~ alcoh + drug, data=results)
Error in optim(par = start, fn = loglikfun, gr = gradfun, method = method, :
initial value in 'vmmin' is not finite
Why is R returning me error in optim()?
What