[R] vglm help

2010-09-03 Thread choonhong ang
Hi All, I am using vglm & the fitted values is NA. As in the R documentation that for Pareto1, if the estimate of k is less than or equal to unity then the fitted values will be NA. what is NA means ? how to solve it ? how to get the k estimate ? (in the R document the estimate of alpha is f.

[R] vglm

2010-08-31 Thread choonhong ang
Hi All, could anybody help me to understand what is this error means ? mydata=read.table("C:/Documents and Settings/angieb/Desktop/CommercialGL/cl_ilf_claimdata.csv",header=TRUE,sep=",") > names(mydata) [1] "ILFTable""liabLimit" "AnnAggLimit" "DedAmt" "Loss" "TIL" > fit=vglm(Loss~1,pa

[R] help - read SAS into R

2010-08-31 Thread choonhong ang
HI All, How to read SAS data directly into R ? Thank you [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-g

[R] A Gamma-GLM with log link

2009-12-07 Thread choonhong ang
Hi, I have a set of data (total number of record = 144,122), and I would like to use gamma-glm with log link to set up a model. IC is number of records IL is paid amount The table below shows that I have 30.578% of the data in the level of "1 - 1000" paid amount 20.320% of the data in the level

[R] detect outliers and high levarage points

2009-03-03 Thread choonhong ang
Hi friends, How to detect outliers and high leverage points for GLM ? Could I use plot(model) (i) "Residuals vs Fitted" graph to detect the outliers ? (ii) "Residuals vs Leverage" graph to detect the high leverage points ? And then remove those points from the data and re-run the model ?

[R] comment on this book "A Handbook of Statistical Analyses Using R by Brian S. Everitt (Author), Torsten Hothorn (Author)"

2009-03-02 Thread choonhong ang
Is this book a good reference to learn R for statistical analysis ? A Handbook of Statistical Analyses Using R by Brian S. Everitt(Author), Torste

[R] select Intercept coefficients only

2009-02-27 Thread choonhong ang
Hi friends, Is there a function to select intercept coefficients only ? When I use "coeficients" it shows me all the coefficients, but I only want a specific coefficients. [[alternative HTML version deleted]] __ R-help@r-project.org mailing li

[R] R question - combine values

2009-02-25 Thread choonhong ang
The district a is the baseline and we observe the difference between District a & b is not significant, we can choose to combine these 2 values. How to write code to combine these 2 value ? > m1=glm(Claims~District+Group+Age+log(Holders),fami ly=poisson,data=mydata) > summary(m1) Call: glm(formul

Re: [R] Insurance data in library(MASS)

2009-02-24 Thread choonhong ang
Hi, In the result shown, the District 1 is used as the base category. How to change to make District 4 as a base category ? On Mon, Feb 23, 2009 at 11:05 AM, choonhong ang wrote: > I have used the insurance data from R library and I have 2 questions: > I use the following: > >

[R] Insurance data in library(MASS)

2009-02-23 Thread choonhong ang
I have used the insurance data from R library and I have 2 questions: I use the following: >library(MASS) >data(Insurance) > m1=glm(Claims ~ District + Group + Age + offset(log(Holders)),data = Insurance, family = poisson) >summary(m1) Call: glm(formula = Claims ~ District + Group + Age + offset(l