similar to: Strange results for Beta Distribution

Hello, I am planning a study with the main point to evaluate the interaction of two treatments, but for ethical reasons one cell is empty, that with patients receaving no treatment at all Treatment B

Dear All, I am looking for an implementation in R of the Appell Hypergeometric function. Any suggestions will be more than appreciated! GP -- dr. Giovanni Parrinello External Lecturer Medical Statistics Unit Department of Biomedical Sciences Viale Europa, 11 - 25123 Brescia Italy Tel: +390303717528 Fax: +390303717488 email: parrinel at med.unibs.it

Dear All, discussing with a statistician of a pharmaceutical company I received this answer about the statistical package that I have planned to use: As R is not a validated software package, we would like to ask if it would rather be possible for you to use SAS, SPSS or another approved statistical software system. Could someone suggest me a 'polite' answer? TIA Giovanni -- dr.

Dear R Fellow-Travellers: What is your recommended way of dealing with a left-censored response (non-detects) in (linear Gaussian) mixed effects models? Specifics: Response is a numeric positive measurement (of volume, actually); but when it falls below some unknown and slightly random value (depending on how the sample is prepared and measured), it cannot be measured and is recorded as 0.

Dear All, from some days when I try to update the installed packages I receive this answer: update.packages(ask='graphics') --- Please select a CRAN mirror for use in this session --- Error in as.list(read.dcf(file = file)[1, ]) : indice fuori limite # index out of bound TIA Giovanni PS: I have already tried to solve the problem reinstalling R or changing the options of

Dear All, I would like to find something ( references, code,..) to implement a comparison of three treatments in an observational study using the 'Propensity Score'. Any help is much appreciated. Thanks! Giovanni -- dr. Giovanni Parrinello Department of Biotecnologies Medical Statistics Unit University of Brescia Viale Europa, 11 25123 Brescia email: parrinel at med.unibs.it Phone:

Dear All, can anyone explain this message Error in load("matched2.RData") : Value of SET_STRING_ELT() must be a 'CHARSXP' not a 'builtin'? Have I the possibility to retrieve the data? TIA Giovanni -- dr. Giovanni Parrinello Department of Biotecnologies Medical Statistics Unit University of Brescia Viale Europa, 11 25123 Brescia email: parrinel at med.unibs.it Phone:

latex(s6a2,title="",caption="Baseline vs zonr",file="",label="Base",long=TRUE,landscape=F, middle.bold=TRUE,here=T, + ,size="smaller[5]",outer.size="smaller",Nsize="smaller",midsize="smaller") "latex" is not reconized \as a internalor external command , un programma eseguibile o un file batch. Warning

Dear All, I am looking for the library QVF [ NONLINEAR MEASUREMENT ERROR PROGRAMS FOR SAS AND SPLUS ]. TIA Giovanni Parrinello dr. Giovanni Parrinello Medical Statistics Unit Department of Medical Sciences University of Brescia Via Valsabbina, 19 25127 Brescia Italy tel.: +39303717528 fax:+39303701157 E-mail: parrinel at med.unibs.it url: http://www.med.unibs.it/dip/stat/index.html

Last week Giovanni Parrinello posted a message asking why various packages were loaded when he loaded an .Rdata file. Brian Ripley replied saying he thought it was because the saved workspace contained a reference to the namespace of ipred. (Correspondence copied below). This begs the question: how did the reference to the namespace of ipred come to be in the .Rdata file? Brian did say it is

Hello to All. I'd want to use a one-way ANOVA. This means that I have only one factor, with, lets say, 5 levels. I made a dataframe, called "DATA", with two Columns: A, that is my response, and it is "class numerical". B, that defines the different levels of my factor, and it is "class factor". If I want to use a fixed effect model, I know that the formula I have

Dear All, I am trying to apply the function 'cox.zph' of the library survival, but I receive this error message: not found the object 'residuals.coxph'. I have re-installed the library 'survival' without any change and also a search with RSiteSearch was unsuccessful.. Any suggestion? TIA Giovanni PS: R 2.2.0 Windows XP HE dr. Giovanni Parrinello Section of Medical

Dear All, simetimes when I load an Rdata I get this message ####### Code: load('bladder1.RData') Carico il pacchetto richiesto: rpart ( Bad traslastion: Load required package-...) Carico il pacchetto richiesto: MASS Carico il pacchetto richiesto: mlbench Carico il pacchetto richiesto: survival Carico il pacchetto richiesto: splines Carico il pacchetto richiesto: 'survival'

Hi all How can I change the limits (xlim or ylim) in a plot that has been already created? For example, consider this naive example curve(dbeta(x,2,4)) curve(dbeta(x,8,13),add=T,col=2) When adding the second curve, it goes off the original limits computed by R for the first graph, which are roughly, c(0,2.1) I know two obvious solutions for this, which are: 1) passing a sufficiently large

Hi, sorry for repeating the question but this is kind of important to me and i don't know whom should i ask. So as noted before when I do a parameter fit to the beta distr i get: fitdist(vectNorm,"beta"); Fitting of the distribution ' beta ' by maximum likelihood Parameters: estimate Std. Error shape1 2.148779 0.1458042 shape2 810.067515 61.8608126 Warning

Full_Name: Hans Peter Wolf Version: 1.2.1 OS: hpux10.20 Submission from: (NULL) (129.70.84.25) dbeta computes a wrong result with parameters (1.3,1) > version platform hppa2.0-hp-hpux10.20 arch hppa2.0 os hpux10.20 system hppa2.0, hpux10.20 status major 1 minor 2.1 year

Dear all, we received a bug report for betareg, that in some cases the optim call in betareg.fit would hang the R session and the command cannot be interrupted by Ctrl-C? We narrowed down the problem to the dbeta function which is used for the log likelihood evaluation in betareg.fit. Particularly, the following command hangs the R session to a 100% CPU usage in all systems we tried it (OS X

Hi, Interpretation problem ! so what i did is by using the: >fit1 <- fitdist(vectNorm,"beta") Warning messages: 1: In dbeta(x, shape1, shape2, log) : NaNs produced 2: In dbeta(x, shape1, shape2, log) : NaNs produced 3: In dbeta(x, shape1, shape2, log) : NaNs produced 4: In dbeta(x, shape1, shape2, log) : NaNs produced 5: In dbeta(x, shape1, shape2, log) : NaNs produced 6: In

Hi all, I compiled 64-bit R on an Apple Mac G5 running OS X, but it failed make check. Simon Urbanek suggested I post results to R-devel. > On Nov 6, 2007, at 10:23 PM, Steven McKinney wrote: > > > Hi Simon, > > > > Would you be able to give more guidance on how to compile 64-bit > > libiconv for Tiger, > > You can get the sources from Apple and compile

Dear All, There is one problem I encountered when I used ADAPT to compute some 2-D integral w.r.t beta density. For example, when I try to run the following comments: fun2<-function(theta){return(dbeta(theta[1],0.005,0.005)*dbeta(theta[2],0.005,0.005))} int.fun2<-adapt(ndim=2,lo = c(0,0), up = c(1,1),functn = fun2,eps = 1e-4) It seems it will take very long time to run. Acturally, I