similar to: Padding in lattice plots

Dear list, I'm trying to fit a glm model using family=poisson(link = "identity"). The problem is that the glm function fits a model with a negative intercept, which sounds like a nonsense to me, being the response a Poisson variable. >From a previous discussion on this list I've understood that the glm function uses IRLS for the fitting without any constraint so it is

Hi all, Is it possible to have a test value for assessing the normality of residuals from a linear regression model, instead of simply relying on qqplots? I've tried to use fitdistr to try and fit the residuals with a normal distribution, but fitdsitr only returns the parameters of the distribution and the standard errors, not the p-value. Am I missing something? Cheers, Federico

Full_Name: Federico Gherardini Version: 1.8.1 OS: Gentoo linux Submission from: (NULL) (81.208.36.89) When I use fix() to edit a matrix R immediately quits with a segmentation fault if I use the copy button in the graphical device window. Clicking on the name of the variables completely messes up the screen too. In this case you have to close the grahpic device and call fix() again. I've

Hi, you can use this simple function: add.col<-function(df, new.col) {n.row<-dim(df)[1] length(new.col)<-n.row cbind(df, new.col) } see this example: > x<-cbind(c(1,2,3),c(4,5,6)) > x [,1] [,2] [1,] 1 4 [2,] 2 5 [3,] 3 6 > y<-c(7,8) > y [1] 7 8 > add.col<-function(df, new.col)

Hi all, I'm trying to read a 1220 * 20000 table in R but I'm having lot of problems. Basically what it happens is that R.bin starts eating all my memory until it gets about 90%. At that point it locks itself in a uninterruptible sleep status (at least that's what top says) where it just sits there barely using the cpu at all but keeping its tons of memory. I've tried with

Hi all I'm trying to use the function glm from the MASS package to do the following fit. fit <- glm(FP ~ rand, data = tab, family = poisson(link = "identity"), subset = rand >= 1) (FP is >= 0) but I get the following error Error: no valid set of coefficients has been found:please supply starting values In addition: Warning message: NaNs produced in: log(x) in contrast

Dear R-users, I am trying to understand what the different padding arguments in trellis.par.set are exactly controlling the space around lattice plots. I have used the following code as a basis for testing but it did not really help me to visualize how the value of each argument changes the margins and the plotting area. I guess a better way to visualize the effects of these padding items

Dear Federico, see: ? shapiro.test(stats) Shapiro-Wilk Normality Test and ? jarque.bera.test(tseries) Jarque-Bera Test They are the most common tests used for normality testing. Ciao Vito Federico Gherardini wrote on Fri Oct 15 14:44:18 CEST 2004: Hi all, Is it possible to have a test value for assessing the normality of residuals from a linear regression model,

Dear All, is any way I can change the colours of the box and umbrella of a bwplot without having to go to the length of: box.rectangle<-trellis.par.get("box.rectangle") box.rectangle$col<-"black" trellis.par.set("box.rectangle", box.rectangle) etc... but straight from the call: bwplot(y ~ x | z, mydata) ? Regards, Federico Calboli

Hi Everyone, I would like to add errorbars to a lattice plot. I already have a function that adds error bars to a plot (originally written by George Gilchrist) which takes as arguments the coordinates of the point and then the size of the error bar. I would like to integrate it in my lattice code but I seem to be stumped. I can put the error bar function inside my code as if it were a panel

Hi, R is very new for me, so excuse if my questions are too basic... BTW, are there any forum where new R users could get help without annoying this huge mailing list ? In following code, I'd like to choose the color for each of the curve diplayed. png(filename = filename, width = 950, height = 600, pointsize = 10, bg = "white") xyplot(HITS+MS1*3+FREQ~TIME

Hi You could try to find an equivalent representation as a string and try to match those. > (A <- cbind(sample(1:2, 10, rep=TRUE), sample(1:2, 10, rep=TRUE))) [,1] [,2] [1,] 1 2 [2,] 1 2 [3,] 1 2 [4,] 2 2 [5,] 1 1 [6,] 1 2 [7,] 1 2 [8,] 1 1 [9,] 1 2 [10,] 1 1 > (B <- unique(A)) [,1] [,2] [1,] 1 2

Dear list, I need to align two plots on top of each other for comparison (they only have the x-axis in common). When the y-labels have a different extent, I cannot find a way to align the x-axes, as illustrated below, > library(grid) > library(lattice) > x <- seq(0, 10, length=100) > y <- sin(x) > y2 <- 10*sin(x) > f <- rep(c("1", "2"),

Dear members of the list, it is likely a stupid question but I cannot find the information neither in R manuals nor in google. I am generating a plot (from hclust results) but I cannot see properly the labels because the default font size is too large. How can I change it? Thanks! Federico

On Wed, 28 Jul 2004 09:53:08 +0200 Uwe Ligges <ligges at statistik.uni-dortmund.de> wrote: > > If your data is numeric, you will need roughly > > 1220 * 20000 * 8 / 1024 / 1024 ~~ 200 MB > > just to store one copy in memory. If you need more than two copies, your > machine with its 512MB will start to use swap space ..... > Hence either use a machine with more

Hi All, I am trying to build a new R package to submit, but it's failing to create a tex manual: R CMD check Biodem * checking for working latex ... OK * using log directory '/home/greatsage/Fede/R-packages/temp/Biodem.Rcheck' * checking for file 'Biodem/DESCRIPTION' ... OK * checking if this is a source package ... OK * Installing *source* package 'Biodem' ... ** R

[ This was posted to the R-packages list (which I moderate), but definitely doesn't belong there. Martin Maechler ] ------- start of forwarded message ------- From: Federico Calboli <f.calboli at ucl.ac.uk> To: r-packages at stat.math.ethz.ch Subject: how to upload Date: 15 Jul 2004 18:30:26 +0100 Dear All, I just finished a add on lib called Malmig, to calculate Malecot

On Thu, Oct 26, 2017 at 7:30 AM, Aki Tuomi <aki.tuomi at dovecot.fi> wrote: > > > On October 26, 2017 at 4:30 PM Federico Bartolucci <federico at aruba.it> > wrote: > > > > > > Hello, > > > > it's the first time for me writing to the list, I'm trying to change the > > location into which the Dovecot's locks are done reserving

Hi All, I have a list of matrices: > x [,1] [,2] [1,] 1 4 [2,] 2 5 [3,] 3 6 > y [,1] [,2] [,3] [,4] [,5] [,6] [1,] 18 21 24 27 30 33 [2,] 19 22 25 28 31 34 [3,] 20 23 26 29 32 35 > z =list(x,y) I want to create a second list that is has a subset each matrix in the list subsetting so I get the 2nd and 3rd row of each (and

Hi, yes sure, we already use proxy with mysql auythentication. On the same server we have the proxy working with mysql and it regularly uses the backends servers, but when we configure the proxy to use LDAP instead it always uses the default localhost and not the backend IPs. Regards -federico Il 21/05/2019 07:25, Aki Tuomi ha scritto: > > > On 20.5.2019 20.06, Federico Bartolucci via