similar to: negative binomial: expected number of events?

Dear experts, I just noticed that I get different results conducting type III anova using drop1 or the Anova command from the car package. I suppose I made a mistake and hope you can offer me some help. I have no idea where I got wrong and would be very grateful for explaination as R is new terrain for me. If I run the commands in line, they produce the same results. But if I run them in

Hello, Is there any R plug-in available for Eclipse other than StatET? StatET doesn't seem to work with the latest release of Eclipse properly, i.e. syntax highlighting isn't enabled for R commands. I already contacted the author some time ago, but have not yet received a response. Would also be great if somebody could tell me whether they have the same problem. I already tried

Hi WINE-people! If you are from germany and are interested in Vista, then have a look at http://www.mediaonline.de/ There you can "buy" Microsoft Windows Vista Intallation DVDs 32bit and 64bit, german version, for 0 ? ! You only pay for shipment. This is already annonced on some places here in the internet. But if you still don't know it, I only want to tell it. The only

Hallo, das Subject sagt schon fast alles. Zum besseren Verst?ndnis der Hintergrund. nachdem ich vor vielen jahren Pathworks-Server durch ein OS/2-Serversystem abgel?st habe,haben wir es auch noch geschafft, die n?chste Evolutionsstufe bis jetzt herauszuziehen. Wir nutzen File- und Printservices. Das ist alles. Die Anzahl der User tendiert gegen 150. Wir sind an zwei Standorten, aber ?ber VPN

Hi, I installed the package Hmisc with the command install.packages("Hmisc") without errors. When I try to load the library with command library(Hmisc) I get the error > library(Hmisc) Error in library(Hmisc) : there is no package called 'Hmisc' > version _ platform i386-pc-linux-gnu arch i386 os linux-gnu system i386,

Dear all, I tried to use "step" function to do model selection, but I got an error massage. What I don't understand is that data as data.frame worked well for my other programs, how come I cannot make it run this time. Could you please tell me how I can fix it? ***************************************************************************************************

Hello, I was wondering if I can plot two curves I get from "density(data)" into one plot. I want to compare both. With the following commad, I just get one curve plotted: plot( density(mydata) ) Sorry for this stupid question but I could not find a solution until now... Antje

Hello, I am interested in plotting my regression analysis results(regression coefficients and standard errors obtained through OLS and Tobit models) in the form of graphs.What is the best way to accomplish this? Thanks. Murali Kuchibhotla [[alternative HTML version deleted]]

I''m looking for a R function which can maximise this logliklihood function, under the constraits a>0 e b>0 f<-function(param){ a<-param[1] b <-param[2] log(prod)-(a*s2)-(b*s)-n*log(1-((0.5*b/sqrt(a))*(exp((b^2)/(4*a)))*((sqrt(pi ))*(1-pnorm(-b/(2*sqrt(a)), mean=0, sd=1)))))} I''ve tried maxlik constrOptim e donlp2 but without success. Thanks so

Hello all, In the course of trying to implement the CDF of an InverseGammaDistribution, I have run across the need for an igamma() function. Several others have needed this function but the answers I have found so far are not totally clear to me. I'm writing for three reasons: 1) to present a small error in the docs 2) to clarify the approach we are expected to take 3) to request,for the

Dear all, I want to design a histogram and I need to have the frequency at certain points. For example I have the following 2 columns: *X Y* 0.1 25 0.4 22 0.45 11 0.55 21 I want the chart to have 4 columns. First column is from 0.0-0.1 (on X) and frequency is 25. Next colum is wider and form 0.1-0.4 with 22 frequency. Next column is narrow with 11 frequency and the last column is

Hi All, Please pardon me if I am missing something obvious here. How do I get the Kaplan-Meier estimate function that is created by survfit and plotted by the code. fit <- survfit(Surv(time, status) , data=aml) plot(fit) That is, I need a function that will give me the survival estimate at a given time: \hat{S}(t). Thanks in advance. Ritwik Sinha ritwik.sinha at gmail.com | +12033042111 |

Hi, all, Does anybody got the most general expression of a unitary matrix? I found one in the book, four entries of the matrix are: (cos\theta) exp(j\alpha); -(sin\theta)exp(j(\alpha-\Omega)); (sin\theta)exp(j(\beta+\Omega)); (cos\theta) exp(j\beta); where "j" is for complex. However, since for any two unitary matrices, their product should also be a unitary matrix. When I

I was reading about the project drawer feature in Textmate, which is Mac only. Is there a similar feature in a Windows based text editor that works with R. This feature sounds really useful. Thanks, Graham [[alternative HTML version deleted]]

Hi All, I have a data output as below.I want to display them in an interactive html report using shiny but the data table is not rendering properly and instead giving NA values. max_usage_hours_per_region<-setNames(aggregate(df3_machine_region$sum_as_hours~df3_machine_region$Region,df3_machine_region,max),c("Region","Sum_as_Hours")) Region Sum_as_Hours 1 Africa 1156.0833

Hi All, I have a data output as below.I want to display them in an interactive html report using shiny but the data table is not rendering properly and instead giving NA values. max_usage_hours_per_region<-setNames(aggregate(df3_machine_region$sum_as_hours~df3_machine_region$Region,df3_machine_region,max),c("Region","Sum_as_Hours")) Region Sum_as_Hours 1 Africa 1156.0833

Hi, Suppose I have a multivariate response Y (n x k) obtained at a set of predictors X (n x p). I would like to perform a linear regression taking into consideration the covariance structure of Y within each unit - this would be represented by a specified matrix V (k x k), assumed to be the same across units. How do I use "lm" to do this? One approach that I was thinking of

Do any of the R libraries have an implementation of the Inverse Error Function (Inverse ERF)? ref: http://mathworld.wolfram.com/InverseErf.html http://functions.wolfram.com/GammaBetaErf/InverseErf/ Thanks, Nathan [[alternative HTML version deleted]]

To R-helpers, R offers the polyroot function for solving mentioned equations iteratively. However, Dr Math and Mathworld (and other places) show in detail how to solve mentioned equations non-iteratively. Do implementations for R that are non-iterative and that solve mentioned equations exists? Regards, Mads Jeppe

I'm having trouble with the standard deviation distribution as shown on http://mathworld.wolfram.com/StandardDeviationDistribution.html . (Eric Weisstein references Kenney and Keeping 1951, which I can't check.) I believe the graphs they show, but when I code the function in R, according to the listed formula, I get very different graphs. Would someone please point out my error or tell