[R] Q: Mean, median and confidence intervals with functions summary boxplot.stats

2007-08-30 Thread Tom Willems
Dear R ussers, My question is, How can my mean be outside the confidence intervals ?! I think i have the answer for it, but i would like to hear some other ideas on it. First my data is not continuose but categorical, it is a titre calculated on a dilution serie. It is stored as a column of

[R] How do i print a main title on a win.graph with several plots?

2007-08-22 Thread Tom Willems
Good Mornig All, How R you today? ;o) I have lots of questions, but i l start with the simplest one, to wich i am shy to say, i did not find the answer. It is the following: When i make a summary plot like for example plot( summary(glm)), i get one window, one main title, and 4 graph's in

[R] Ans on: How do i print a main title on a win.graph with several plots?

2007-08-22 Thread Tom Willems
Thanks Mr.Pr. Ripley, For pointing me in the good direction, i use win.graph() so i get an overview of the different graphs. Until now, i v had no problems with it, hope it stays that way. for those whom have the same porblem, this is what i do: win.graph(); op - par(mfrow = c(1,2) , oma=

[R] Q: combine 2 data frames with missing values

2007-08-20 Thread Tom Willems
hello R ussers, i have the same problem with my data, for aal the different variables, i have the same number of cases, but the are often out of detectionlimits so they produce na's . so the data looks like this: casevar1var2var3var4 ... 1 9 9 13 11 2

Re: [R] Does anyone.... worth a warning?!? No warning at all

2007-08-20 Thread Tom Willems
dear Mathew mean is a Generic function mean(x...) in wich x is a data object, like a data frame a list a numeric vector... so in your example it only reads the first character and then reports it. try x = c(1,1,2) mean(x) kind regards, Tom Disclaimer: click here [[alternative

[R] Q: how to extract coefisients from one glm and implement them in to an other glm?

2007-08-13 Thread Tom Willems
Dear R ussers, I'm still working with the categorical models, and i can not do cross validation on them because odf their nature. not cv works perfectly on other models so i wonderd if it is possible to extract the beta's from a categorical model, then use their value in a 0/1 glm. Their

[R] Logistic Models, Cross Validation, error tables , Monte Carlo Simulation?

2007-08-08 Thread Tom Willems
Thanks Mr. Ellison, Your remark helped solve my error table problem. However, I found a new one. Now that I have my error tables, i realised that it is no good statistical practise to validate a model, based one one error table. So i should use a tool like K-fold CV. ex: binary_model -

[R] logistic models validation for categorical data

2007-08-07 Thread Tom Willems
, but also that it does not give reference to the number of test samples, wich is clear in table 1 ( 30 test samples). Kind regards, Tom. Tom Willems CODA-CERVA-VAR Department of Virology Epizootic Diseases Section Groeselenberg 99 B-1180 Ukkel Tel.: +32 2 3790522 Fax: +32 2 3790666 E-mail

[R] question about logistic models (AIC)

2007-08-03 Thread Tom Willems
as - AIC . R displays it as a positive number, does this mean that a large AIC gives a small - AIC , so the bigger the better? Kind regards, Tom. Tom Willems CODA-CERVA-VAR Department of Virology Epizootic Diseases Section Groeselenberg 99 B-1180 Ukkel Tel.: +32 2 3790522 Fax: +32 2 3790666 E