similar to: scoping problem?

*Dear R-help: My name is Rodrigo and I have a question with nlme package in R to fit a mixed beta regression model. The details of the model are: Suppose that:* *j in {1, ..., J}* *(level 1)* *i in {1, ..., n_j}* *(level 2)* *y_{ij} ~ Beta(mu_{ij} * phi_{ij}; (1 - mu_{ij}) * phi_{ij}) y_{ij} = mu_{ij} + w_{ij} * *with* *logit(mu_{ij}) = Beta_{0i} + Beta_{1i} * x1_{ij} + b2 * x2_{ij}

*Dear R-help: My name is Rodrigo and I have a question with nlme package in R to fit a mixed beta regression model. I'm so sorry. In the last email, I forgot to say that W is also a unknown parameter in the mixed beta regression model. In any case, here I send you the correct formulation. ** Suppose that:* *j in {1, ..., J}* *(level 1)* *i in {1, ..., n_j}* *(level 2)* *y_{ij} ~

n_10;p_0.5;jjx_0:n;qbinom(pbinom(jjx,n,p),n,p) # This one works as expected n_100;p_0.5;jjx_0:n;qbinom(pbinom(jjx,n,p),n,p) # This one causes severe problems I cannot interrupt using ESC and I finally have to resort to the Windows Task manager to kill the R session. A friend of mine told me that he faced similar problems under Unix. --please do not edit the information below-- Version:

Everyone, One of our developers was kind enough to insert some bug checking into the mod_auth_pam and mod_auth_sys_group so that we could see a little more of what was going on with our authentication failures. Here is what we just saw. Two of our users NA\connelmp and NA\guminssa both started getting messages that they were not part of the required group. Here is the log for

Is there facility in R to read binary file? In Splus scan() has 'width=' argument and AsciiToInt() command to translate the binary data. Is there something similar in R? Thanks, -Yudi- -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or

I've had a play with this and, due to my own short-comings, remain none the wiser. In particular, I'm not sure what value of 'spar' is consistent with the magic lambda=1/1600 for quarterly data. I initially interpreted spar as lambda and tried setting spar=1/1600. This results in almost no smoothing while spar=1600 causes an error. The smooth.spline function seems to want

Hi, I posted some days ago a question concerning the computation of lambda in the smooth.spline function (which I repreat at the bottom of the mail) given _df_ . Unfortunately the documentation is not clear to me. Maybee someone can help to answer in my view the basic question: If the penalized log likelihood is L = (y - f)' W (y - f) + lambda c' Sigma c how the _lambda_ in the above

This came up on r-help but indicates a bug. dnbinom(x,n,p) calls dbinom_raw(n-1,...) which returns 0 for n<1. -thomas ---------- Forwarded message ---------- Date: Thu, 08 Feb 2001 17:10:23 +0000 From: Yudi Pawitan <yudi@stat.ucc.ie> To: Mark Myatt <mark@myatt.demon.co.uk> Cc: R-Help <r-help@stat.math.ethz.ch> Subject: Re: [R] Goodness of fit to Poisson / NegBinomial

Everyone, Here is a challenge for all of you samba experts! Lately I have been seeing a problem where winbind is not correctly identifying a user as a member of a group he most certainly belong to. This is with a Domain Local group so I know samba should support it. Users access a HTTPS (SSL) webpage that is secured by a Domain Local group. Sometimes they get in,

Hi, I am having an issue with a custom panel for lattice. The problem comes when I try passing a groups argument. Here is the custom panel, a wrapper for smooth spline. I copied panel.loess and replaced the loess arguments with smooth.spline(). [Note: I would like to use the cross-validation fitting properties of smooth.spline.] library(lattice)

Hi is it possible to fit a trend line (or some other panel function) through each of multiple data series plotted on the same graph? Specifically, while one can do something like xyplot(a+b+c~x) which plots three series, a,b & c, but can one automatically fit lines through each of them? I suppose one could generate three more variables afit, bfit, and cfit with a model & predict and

Hello, I come from a non statistics background, but R is available to me, and I needed to test an implementation of smoothing spline that I have written in c++, so I would like to match the results with R (for my unit tests) I am following http://www.nabble.com/file/p25569553/SPLINES.PDF SPLINES.PDF where we have a list of points (xi, yi), the yi points are random such that: y_i = f(x_i) +

Dear list, if I do smooth.spline(tmpSec, tmpT, all.knots=T) with the attached data, I get this error-message: Error in smooth.spline(tmpSec, tmpT, all.knots = T) : smoothing parameter value too small If I do smooth.spline(tmpSec[-single arbitrary number], tmpT[-single arbitrary number], all.knots=T) it works! I just don't see it. It works for hundrets other datasets, but not for

Hi, I am trying to compile R-2.2.1 on Solaris 2.9 with a 64-bit build. Following the instructions in "R Installation and Adminstration", I changed the following settings in "config.site": CC="gcc -m64" F77="g77 -64" CXX="g++ -m64" LDFLAGS="-L/usr/local/lib/sparcv9 -L/usr/local/lib" But I got the following error messages:

Can someone please show me how to smooth time series data that I have in the form of a zoo object? I have a monthly economies series and all I really need is to see a less jagged line when I plot it. If I do something like s <- smooth.spline(d.zoo$Y, spar = 0.2) plot(predict(s,index(d.zoo)), xlab = "Year") # not defined for Date objects and if I do something like

I'm putting a title over different plots, where 'that' varies. How do I make the value appear? Example: that <- 1 plot(1:10) title(expression(paste(hat(theta),'= ',that))) Yudi Pawitan yudi at stat.ucc.ie Department of Statistics UCC Cork, Ireland Ph 353-21-490 2906 Fax 353-21-427 1040 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-

>I use WinNT. You have to "launch" R from a DOS-shell window; you 1st >change to the directory where you have RGui.exe, then type : > Rgui --vsize 15M --nsize 1000k > >and R opens with increased memore size. It works here with my NT box, and >will probably do with Win95 (though I think that the DOS-shells are not >strictly similar in both OSs but this may not be a

(Sorry, I had the previous emails on a different machine.) Thanks to Prof Ripley for diagnosing the problem with inconsistent capitals under Win98. It's 'interesting': * I had ...\acepack\Src\Ace.f Avas.f and Rlsmo.f (tar'ed in a Unix machine and Ws_FTP'ed to the PC) * when the directory name was changed to ...\acepack\src it still didn't work *

To whom it may concern, I don't know whether this is really a bug with the Pspline package or only a problem with my installation. Things work fine in Linux but not in Mac OS X (Darwin). Both system run the latest public versions of R and Pspline. predict.smooth.Pspline produces only NaN instead of predicted values when norder>2: > library (Pspline) > tt <- seq

Dear All, I know how to export graphics as pdf files and then how to include them in LaTeX documents. However, I do not know how to do in order to have the text of the graphics written with the font selected for the LaTeX document. Is that possible? Thanks in advance, Paul