Displaying 20 results from an estimated 1000 matches similar to: "how to call out the estimated parameter values"
2006 Oct 27
0
VGAM package released on CRAN
Dear useRs,
upon request, the VGAM package (currently version 0.7-1) has been
officially released on CRAN (the package has been at my website
http://www.stat.auckland.ac.nz/~yee/VGAM for a number of years now).
VGAM implements a general framework for several classes of
regression models using iteratively reweighted least squares
(IRLS). The key ideas are Fisher scoring, generalized linear
and
2007 Nov 29
2
convert an S plus file to R?
hi!
i send again my question because there was a problem earlier that
someone did not see my attached file. If you really can't download it,
this is the attached file.
Please help me how to convert this S plus file to R.
Is there a quick method to do it? I don't have an S plus
installer here.
---------------------------------------------------
# Computes a possible choice for
2007 Jul 27
1
R codes for g-and-h distribution
hi!
I would like to ask help how to generate numbers from g-and-h distribution. This distribution is like normal distribution but span more of the kurtosis and skewness plane. Has R any package on how to generate them?
Any help will be greatly appreciated. Thank you so much!
Form,
Filame Uyaco
---------------------------------
[[alternative HTML version deleted]]
2012 Jan 10
1
S4 summary method not being called (VGAM)
The symptom triggering this email is that an S4 summary method sometimes refuses to be invoked, even when a package is explicitly loaded, if the first load of the package is implicit. It may or may not be specific to 'summary' methods and/or the 'VGAM' package. I've sent to R-devel because (i) it looks like some kind of bug to me, but I'm not sure; (ii) it's not
2013 Apr 17
1
Bug in VGAM z value and coefficient ?
Dear,
When i multiply the y of a regression by 10, I would expect that the
coefficient would be multiply by 10 and the z value to stay constant. Here
some reproducible code to support the case.
*Ex 1*
library(mvtnorm)
library(VGAM)
set.seed(1)
x=rmvnorm(1000,sigma=matrix(c(1,0.75,0.75,1),2,2))
2011 Aug 23
1
P values for vglm(zibinomial) function in VGAM
Hi ,
I know this question has been asked twice in the past but to my knowldege,
it still hasn't been solved.
I am doing a zero inflated binomial model using the VGAM package, I need to
obtain p values for my Tvalues in the vglm output. code is as follows
> mod2=vglm(dmat~Season+Diel+Tidal.phase+Tidal.cycle,zibinomial, data=mp1)
> summary(mod2)
Call:
vglm(formula = dmat ~ Season +
2007 May 31
0
VGAM package
Hi, R-users
Could someone help me to understand this following error. I'm using vglm
function in VGAM package
Best regards and thank you for your ehlp
########
mydata <- read.table("Data2_overruns.csv", sep =";", header = T,
row.names=NULL)
> attach(mydata)
>
> y <- mydata$cat.event
> phase.vol <-mydata$phase.vol
> pilote <- mydata$pilote
2008 Apr 18
2
rzinb (VGAM) and dnbinom in optim
Dear R-help gurus (and T.Yee, the VGAM maintainer) -
I've been banging my head against the keyboard for too long now, hopefully someone can pick up on the errors of my ways...
I am trying to use optim to fit a zero-inflated negative binomial distribution. No matter what I try I can't get optim to recognize my initial parameters. I think the problem is that dnbinom allows either
2009 Jun 05
2
p-values from VGAM function vglm
Anyone know how to get p-values for the t-values from the coefficients
produced in vglm?
Attached is the code and output ? see comment added to output to show
where I need p-values
+ print(paste("********** Using VGAM function gamma2 **********"))
+ modl2<-
vglm(MidPoint~Count,gamma2,data=modl.subset,trace=TRUE,crit="c")
+ print(coef(modl2,matrix=TRUE))
2009 Jun 16
0
Generation from COX PH with gamma frailty
Hello,
I want to generate data set from Cox PH model with gamma frailty effects.
theta(parameter for frailty distribution)=2
beta=1.5
n=300
cluster size=30
number of clusters=10
I think I should first generate u from Gamma(Theta,theta) and then using
this theta I could not decide how I should generate the survival times?
Is there any package for this? or any document you could suggest?
Any
2012 Apr 04
0
multivariate ordered probit regression---use standard bivariate normal distribution?
Hello.
I have yet to receive a response to my previous post, so I may have
done a poor job asking the question. So, here is the general question:
how can I run a run a multivariate (more than one non-independent,
response variables) ordered probit regression model? I've had success
doing this in the univariate case using the vglm() function in the
VGAM package. For example:
2007 Oct 29
1
VGAM and vglm
Hi Folks,
I wonderif someone who is familiar with the details
of vglm in the VGAM package can assist me. I'm new
to using it, and there doesn;t seem much in the
documentation that's relevant to the question below.
Say I have a vector x of 0/1 responses and another
vector y of 0/1 responses, these in fact being a
bivariate set of 0/1 responses equivalent to
cbind(x,y).
E.g.
2007 Jan 06
2
Using VGAM's vglm function for ordinal logistic regression
R-Experts:
I am using the vglm function of the VGAM library to perform proportional
odds ordinal logistic regression. The issue that I would like help with
concerns the format in which the response variable must be provided for
this function to work correctly. Consider the following example:
------
library(VGAM)
library(MASS)
attach(pneumo)
pneumo # Inspect the format of the original dataset
2011 Jun 22
2
VGAM constraints-related puzzle
Hello R users,
I have a puzzle with the VGAM package, on my first excursion into
generalized additive models, in that this very nice package seems to
want to do either more or less than what I want.
Precisely, I have a 4-component outcome, y, and am fitting multinomial
logistic regression with one predictor x. What I would like to find
out is, is there a single nonlinear function f(x) which acts
2007 Jul 16
2
Error while fitting Partial Proportional Odds model using vglm
Dear R developers:
I am trying to fit a PPO model using vglm from the library VGAM, and get an
error while executing the code. Here is the data, code, and error:
Data: first row is the column names. a = age, and 1,2,3, 4 and 5 are
condition grades.
a 1 2 3 4 5
1 1 0 0 0 0
2 84 2 7 10 2
3 16 0 6 6 2
4 13 0 3 4 0
5 0 0 0 1 0
Library(VGAM)
2010 Jul 05
1
question concerning VGAM
Hello everyone,
using the VGAM package and the following code
library(VGAM)
bp1 <- vglm(cbind(daten$anzahl_b, daten$deckung_b) ~ ., binom2.rho,
data=daten1)
summary(bp1)
coef(bp1, matrix=TRUE)
produced this error message:
error in object$coefficients : $ operator not defined for this S4 class
I am bit confused because some day ago this error message did not show up
and
2009 Mar 26
0
(Interpretation) VGAM - FRECHET 3 parameters by maximum likelihood estimation for
Dear R Helpers
This is the R code (which I have slightly changed) I got in VGAM package for estimating the parameters of FRECHET.
_________________________________________________________________
y = rfrechet(n <- 100, shape=exp(exp(0))) # (A)
fit3 = vglm(y ~ 1, frechet3(ilocation=0), trace=TRUE, maxit=155) # (B)
coef(fit3, matrix=TRUE) # (C)
2009 Nov 04
1
vglm(), t values and p values
Hi All,
I'm fitting an proportional odds model using vglm() from VGAM.
My response variable is the severity of diseases, going from 0 to 5 (the
severity is actually an ordered factor).
The independent variables are: 1 genetic marker, time of medical observation,
age, sex. What I *need* is a p-value for the genetic marker. Because I have ~1.5
million markers I'd rather not faffing
2009 Jan 29
0
Problem VGAM and Predict
Hello,
since I installed the package VGAM I have problems useing the predict for othere methods.
for example I have a model from glm and polr the command predict(model) I get the error: unable to find an inherited method for function "predict", for signature "polr".
Has perhaps anybody a solution, because Iwould need vglm and also other methods like tree in a loop.
Thanks a
2007 Jun 13
1
VGAM Pareto
I would like to fit a Pareto Distribution and I am using the following codes
fit=vglm(ycf1 ~ 1, pareto1(location=alpha), trace=TRUE, crit="c")
fitted(fit)
But the fitted values turn out to be the same for each observation. I guess
the problem is with "ycf1 ~ 1",
I would be grateful if anyone can give me some advice on how to define the
formula.
Many thanks
--
View this