similar to: glm with a distribution free family

Hello Mark, Did you ever solve this problem? I saw your post but no replies. I'm in the exact same situation. If you have any pointers on where to get this I'd appreciate it. - Paul Hess >From: "Mark Mestdagh" <mestdagh@dma.be> >To: "samba digest" <samba@samba.anu.edu.au> >Subject: TCP/IP >Date: Thu, 26 Mar 1998 21:40:59

How might you fit a generalized linear model (glm) with variance = mu+theta*mu^2 (where mu = mean of the exponential family random variable and theta is a parameter to be estimated)? This appears in Table 2.7 of Fahrmeir and Tutz (2001) Multivariate Statisticial Modeling Based on Generalized Linear Models, 2nd ed. (Springer, p. 60), where they compare "log-linear model fits to

I am unfamiliar with R and I’m trying to do few statistical things like GLM and GAM with it. I hope my following questions will be clear enough: My datas ( y(i,j ))are run off triangles for example : J=1 J=2 J=3 I=1 1 2 3 I=2 4 5 I=3 6 My model is : E[y(i,j)] =m(i,j) Var[y(i,j)] =constant *m(i,j) Log(m(i,j)) = eta (i,j) eta (i,j) = c + alpha(i)

Dear R-Users, without going into details I tried to prepare a simple example to show you where I would need help. In particular I prepare two examples-template for a study I'm conduction on discrete-time methods for survival analysis. Each of this example has two datasets which are basically equal, with the exception that in the former one has individual data and in the latter one aggregated

Hi! Does anybody know a string function that would calculate how many characters two strings share? I.e. ("Hello World","Hello Peter") would be 7. Thanks. Laetitia

Hello, Is there a simple and fast function that returns a vector of the number of replications for each object of a vector ? For example : I have a vector of IDs : ids <- c( "ID1", "ID2", "ID2", "ID3", "ID3","ID3", "ID5") I want the function returns the following vector where each term is the number of replicates for the

Hi all, I have this application that lists videogames, which i am rewriting to follow REST. I need to have this url : /games/xbox360/halo-3 ( map.connect "/ games/:platform_id/:permalink",:controller => "games",:action => "show" ) This is no problem with map.connect, but with map.resources i can not make this. map.resources :games do |game|

Hi, I've had a little problem for several weeks now. It is annoying and therefore I will ask for help now: When I write a script with several iterations, I make it write out a text file to save the data during the run. For example I write: if (i %% 25) write.table(output,"temporary_output.txt") Later on, when I read in this output and want to calculate things, R complains that

Hi R people! I am looking for some suggestions writing an if-else function. The idea is to characterize different plots containing counts of variables (here parasites). If a plot has a count equal or higher than 4 for any parasite the function should return a 1 else a 0. Later I can loop the function over all plots. Here I have a little subset of my data: VariablePAR Plot1

Help! Being a bit of a novice, please bear with me if this is a stupid question! I am trying to fit a saturated model to some count data that I have: model<-glm(COUNT~SP*LOC*COL*TIME*TREAT,poisson) but R keeps on crashing and coming up with (occasionally before crashing) an error that states: Algorithm did not converge in: (if(is.empty.model(mt)) glm.fit.null else glm.fit)(x = X, y = Y,

Hello everyone, How can I specify in tcltk file browser the initial directory to "MyComputer" in Windows where Drives and Partition are accessible? And just a little question if anyone knows, is there a way to use the function choose.files under windows to select a directory? Thanks a lot for your help. Laetitia.

Buenas. Sé que en R hay multitud de datasets y me haría falta alguno que trataran de variables relacionadas con salud, sobre todo para aprender más acerca de cómo realizar una regresión logística binomial o multinomial. Gracias..

Hi, I have R installed under a Mandrake linux system and I don't have shell utilities any more under my R console such as completion when writing a file path, back and forth in the history, bindkeys... Moreover when I quit R by saving, no .Rhistory file is created while the .Rdata is. I don't get how this work, I thought it was based on the user unix shell but it does not seem. Any

Hello World, I have a function that makes pairwise comparisons between two strings. I would like to apply this function to my data (which consists of columns with different strings) in the way that it compares the first with the second entry, and then the third with the fourth, and then the fifth with the sixth, and so on down each column... So (2x-1) and (2x) would be the different entries to be

Hi, Can anybody show me how to extract all columns in my dataset that are polymorphic? Or phrased in another way I would like to delete all columns that have no more than one letter in it (that are monomorphic). Thank you. Laetitia -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: seqCol.txt URL:

In the previous versions of R (2.6.1), when a vector of NA was given to the functions 'sd' or 'var' with parameter na.rm = TRUE, it used to return NA. Now (2.7.1) it returns an ERROR : Example in 2.6.1: > sd(c(NA, NA, NA, NA), na.rm = TRUE) [1] NA Example in 2.7.1: > sd(c(NA, NA, NA, NA), na.rm = TRUE) Error in var(x, na.rm = na.rm) : paires d'éléments

Dear Wu Gong and Peter Ehlers, thank you very much for your help debugging my script. Now I have a general following up question: Is there a straightforward way to rearrange the following dataset so that all first letters of each column will be combined in one column, all the second letters in a second column, all the third ones in a third column and so on, resulting in 7 columns, i.e. for

Hi! I need some help to finish my script. I have two tables that I combine randomly to produce a third table. This I do for hundreds of iterations. In the output file I get all the simulated tables after each other. It looks like this (in this case 3 iterations): output file: [[1]] [,1] [,2] [,3] [1,] "GM030005" "WI920024" "CCCC" [2,]

Hello all, I am wondering if there is a way to do a spatial error probit/logit model in R? I can't seem to find it in any of the packages. I can do it in MATLAB with Gibbs sampling, but would like to confirm the results. Ideally I would like to use this model to predict probability of parcel conversion in a future time period. This seems especially difficult in a binary outcome model

Hi everyone, I want to write a function that modify directly variables passed as parameters (the equivalent in C language of *ptr/&ptr) so that I don't have to return a list and to reaffect all my variables. Is it possible to do so in R? Thanks a lot. Laetitia Marisa.