[R] Inflated Array

2006-07-18 Thread Hadassa Brunschwig
Hi R-users! I am trying to create a what I call inflated array (maybe there is already some other name for that). It is an array that changes dinamically its dimensions, e.g. the higher the number of third dimensions, the more rows in the array. So for example the array could look like the

[R] Negative Binomial: Simulation

2006-07-14 Thread Hadassa Brunschwig
Hi R-Users! I fitted a negative binomial distribution to my count data using the function glm.nb() and obtained the calculated parameters theta (dispersion) and mu. I would like to simulate values from this negative binomial distribution. Looking at the function rnbinom() I was looking at the

[R] Negative Binomial: Simulation

2006-07-14 Thread Hadassa Brunschwig
Hi R-Users! (sorry about the last email) I fitted a negative binomial distribution to my count data using the function glm.nb() and obtained the calculated parameters theta (dispersion) and mu. I would like to simulate values from this negative binomial distribution. Looking at the function

Re: [R] Negative Binomial: Simulation

2006-07-14 Thread Hadassa Brunschwig
Addition to last email: The problem basically arises also when I want to calculate (plot) the cumulative distribution of the negative binomial with the calculated parameters mu and theta! On 7/14/06, Hadassa Brunschwig [EMAIL PROTECTED] wrote: Hi R-Users! (sorry about the last email) I

[R] tick mark intervals

2006-03-19 Thread Hadassa Brunschwig
Hi everyone! This must be a common question but I have not found an answer to it in the archives yet. I am producing four plots with par(mfrow=c(2,2)). The x-axis is the same for all of them but the y-axis is different. What I would like to do is to have a different range of the y-axis and

[R] R2WinBUGS: Comparison to WinBUGS

2005-11-01 Thread Hadassa Brunschwig
... Thanks. Hadassa -- Hadassa Brunschwig Birmannsgasse 10A CH-4055 Basel Switzerland Phone: +41 78 797 6065 Email: [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] R2WinBUGS: Comparison to WinBUGS

2005-09-28 Thread Hadassa Brunschwig
- read.bugs(mcmcA[1]) plot(codaA1) THANKS A LOT!! -- Hadassa Brunschwig Birmannsgasse 10A CH-4055 Basel Switzerland Phone: +41 78 797 6065 Email: [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] [Fwd: R2WinBUGS: Comparison to WinBUGS]

2005-09-28 Thread Hadassa Brunschwig
should rely on. Thanks a lot. -- Hadassa Brunschwig Birmannsgasse 10A CH-4055 Basel Switzerland Phone: +41 78 797 6065 Email: [EMAIL PROTECTED] Quoting Sibylle Sturtz [EMAIL PROTECTED]: This is due to the following: In bugs(), the default for thinning is n.thin = max(1, floor(n.chains

Re: [R] R2WinBUGS: Data loading error

2005-09-22 Thread Hadassa Brunschwig
Thank you for your answer. I found out the mistake and it finally works!! Hadassa -- Hadassa Brunschwig Birmannsgasse 10A CH-4055 Basel Switzerland Phone: +41 78 797 6065 Email: [EMAIL PROTECTED] Quoting Uwe Ligges [EMAIL PROTECTED]: Hadassa Brunschwig wrote: Hi R-Help! I am

[R] survreg: fitting different location parameters

2005-01-24 Thread Hadassa Brunschwig
Hi R-Help! My question: I have lifetime/failure data of machines with different stress levels and i think an weibull/extreme value distribution would fit this data. So I did: model1 - survreg(Surv(lfailure)~stress,data=steel,dist=extreme) (where lfailure=log(failure)) Now I would like to

[R] adding to Trellis plot

2003-08-14 Thread Hadassa Brunschwig
Hello everyone (Im not sure if i already posted this problem) I have grouped data and i used a Trellis plot to show the curve of a fitted model for each group. Additionally i have a vector whose values are the highest points the curve can reach (the curve rises exponetially). Each value of the

[R] function lme

2003-08-12 Thread Hadassa Brunschwig
Thanks to everyone who replied to my first problem. Here is the second one: Is the function lme (Mixed Model) also in a foreign library(shouldnt, as this is a basic statistical function)? I am trying to use lme in R and it replies that it couldnt find the function lme...weird... Thanks for

[R] Importing Data

2003-08-08 Thread Hadassa Brunschwig
Im trying to import data from an excel sheet or a sas file to R...im not succeeding. Apparently the function read.xport for reading a SAS file doesnt exist. What do i have to type in EXACTLY to read from an excel sheet(i guess i would be using read.table?)? Thanks in advance for an answer