similar to: estimate inverse gaussian in R

Dear R Helpers Which package is available for estimatine the parameters of three parameter FRECHET distribution. Also, how to generate the random numbers for Frechet using these three estimated parameters. Thanking in advance Maithili

Hi, I am traying to implement an STAR model, but I have some problems. I am following the instruction of the model, that they are in: http://bm2.genes.nig.ac.jp/RGM2/R_current/library/tsDyn/man/star.html that they are from: http://bm2.genes.nig.ac.jp/RGM2/pkg.php?p=tsDyn The model is: star(x, m=2, noRegimes, d = 1, steps = d, series, rob = FALSE, mTh, thDelay, thVar, sig=0.05, trace=TRUE,

Hi All, I want to produce a similar "Forest Plot" as it is on the following link, but my data would be having only two columns (one for "Estimate" and other for "Std. Dev"). Can anyone suggest some function() {Package} which can take such file as an input and give following forest plot:

Hi, I hope somebody can help me with this issue: I am doing population pyramids using the barplot command, so in the left side I have male age structure and in the right side the female age structure. To plot the male age structure I put the data in negative numbers. Now, I want to change the sign in the bar plot in such way that I have no-sign numbers, both in left and right side of the graph. I

Hi! I'm trying to compare maximum likelihood estimator with maximum spacing estimator (in Generalized extreme value distribution). My problem is the following: i know how to determinate maximum likelihood estimator on R, but i dont find a way to do the maximum spacing estimator...:-( Any help would be great:-) -- View this message in context:

Odette > Wondering how can I generate "Akaike weight" with R? I know the description, > but is there any function to generate by R on the web-site or R library? > I am using GLM or GLMM (family=binomial), so would be appreciated if you > help me. You could have a look at this. http://bm2.genes.nig.ac.jp/RGM2/R_current/library/aod/man/summary.aic.html Which is in the OAD

Hello I'm running into a variety of limitations with my computer complaining that it doesn't have sufficient memory. It is a Dell Latitude D630 Duo Core CPU T7700 @ 2.40 GHz, with 3.50 GB RAM. I'm working with relatively large dataframes > dim(AllYears) returns 4071840 records with 13variables sometimes this complains that it can not allocate vectors of 16.0 Mb and will

Hi. I have data which i would want to assess the degree of agreement between two assays, e.g., to evaluate reproducibility or for inter-rater reliability. I have used the Pearson product-moment correlation coefficient. It looks good ranginging between 0.90 to 0.998. Though this looks good. I am told the Concordance correlation coefficient will give a better picture of how reproducible the assay

Two questions: 1. inverse.gaussian is up there as one of the glm families, but do people ever use it? There is no inverse.gaussian in the R distribution family, and when I checked McCullagh & Nelder, it only appeared twice in the book (according to subject index), once in the table on p. 30 and once on p. 38 in a passing sentence. Is there a good reference on this distribution? 2. When I

Hello, I'm Emanuela, I'm implemented a survival analysis and I'm trying to use a frailty model with inverse gaussian distribution, but I'm not able to find the right code, because it seems to be only a gamma and a gaussian distribution. Is there also the inverse gaussian distribution? Thanks a lot Emanuela Rossi

Dear list, ? I am trying to program semi-random movement within a circle, with no particles leaving the circle. I would like them to bounce back when they come to close to the wall, but?I don't seem to be able to get this right. ?Would somebody?be able to?give me a hint ? This is my code so far, the particle starts at some point and moves towards the wall, but I don't get the

Hi, What is the function for inverse normal transformation? Thanks, Carol [[alternative HTML version deleted]]

Dear, I want to fit an inverse gaussion distribution to a data set. The predictor variables are gender, area and agecategory. For each of these variables I've defined a baseline e.g. #agecat: baseline is 3 data<-transform(data, agecat=C(factor(agecat,ordered=TRUE), contr.treatment(n=6,base=3))) The variable 'area' goes from A to F (6 areas: A,B,C,D,E,F) How can i

Dear R-people Does someone have a routine to ngenerate inverse-gaussian random numbers. I am thinking of something similar to rinvgauss, pinvgauss etc. in S-plus. best regards Helgi -- Helgi Tomasson FAX: 354-552-6806 University of Iceland PHONE:354-525-4571 Faculty of Economics and Business Administration

I needed to generate some data from the Inverse Gaussian distribution but it is not in R. A google search took me to http://www.maths.uq.edu.au/~gks/s/invgauss.html which contains the documentation for the d, p, r, q functions with a download link of the functions. The functions work fine in R. I am curious why the facilities are not included in R. Is it for legal reasons or just the

Dear all, I'd like tweaking the ... arguments that one user can pass in my function for fitting a model. More precisely, my objective function is (really) problematic to optimize using the "optim" function. Consequently, I'd like to add in the "control" argument of the latter function a "ndeps = rep(something, #par)" and/or "parscale =

Problem on flexmix when trying to apply signature developed in one model to a new sample. Dear R Users, R Core Team, I have a problem when trying to know the classification of the tested cases using two variables with the function of flexmix: After importing the database and creating a matrix: BM<-cbind(Data$var1,Data$var2) I see that the best model has 2 groups and use: ex2

I have a 64-bit windows box - Intel Xeon CPU E7340 @ 2.4GHz 31.9GB of RAM I have R 2.11.1 (64bit) running on it. My csv data is 3.6 GB (with about 15 million obs, 120 variables.) ------------------------------------------------ I have successfully imported the data above into R. No problem. Now I am trying to run 'rpart' on my data. But I got the following error : Error: cannot

Hello, When I put in the following script line: gap.boxplot(CLI3, CLI4, CLI5, CLI6, CLI7, gap=list(top=c(8000,280000), bottom=c(0,250)), range=50, outline=TRUE, par(ask=FALSE) I get a '+' telling me I am missing something. I have tried adding ')', 'width=NULL', etc and then I get this error: Error: unexpected symbol in: "gap.boxplot(CLI3, CLI4, CLI5, CLI6, CLI7,

Hi, I have the two vectors mean and sd of individual columns, but I am unsure how to generate bar charts with the standard deviation, even after looking the help of barplot and barplot.2. [[alternative HTML version deleted]]