similar to: glm error message when using family Gamma(link="inverse")

I have written a function that contains runs lm() vif() and glm() When the glm() blows up with an error message, I don't get the output from either the lm() or vf() even thought neither lm() nor vif() have any problems . How can I force the function to print sequential results rather than wait for the entire function to complete before listing the functhion's output? Thanks, John

Forgive my resending this post. To data I have received only one response (thank you Bert Gunter), and I still do not have an answer to my question. Respectfully, John Windows XP R 2.12.1 contrast package. I am trying to understand how to create contrasts for a model that contatains an interaction. I can get contrasts to work for a model without interaction, but not after adding the

I am running my code in a loop and it does not work but when I run it outside the loop I get the values I want. n <- 1000; # Sample size fitglm <- function(sigma,tau){ x <- rnorm(n,0,sigma) intercept <- 0 beta <- 0 ystar <- intercept+beta*x z <- rbinom(n,1,plogis(ystar)) xerr <- x + rnorm(n,0,tau) model<-glm(z ~ xerr, family=binomial(logit))

I am new to R and I am trying to do a monte carlo simulation where I generate data and interject error then test various cut points; however, my output was garbage (at x equal zero, I did not get .50) I am basically testing the performance of classifiers. Here is the code: n <- 1000; # Sample size fitglm <- function(sigma,tau){ x <- rnorm(n,0,sigma) intercept <- 0 beta

library(ROCR) n <- 1000 fitglm <- function(iteration,intercept,sigma,tau,beta){ x <- rnorm(n,0,sigma) ystar <- intercept+beta*x z <- rbinom(n,1,plogis(ystar)) xerr <- x + rnorm(n,0,tau) model<-glm(z ~ xerr, family=binomial(logit)) *int*<-coef(model)[1] *slope*<-coef(model)[2] # when add error you are suppose to get slightly bias slope. However when I change

Dear r-community, Today I have completed my PhD.  I would like to take this opportunity to thank the r-community for helping me with the r-coding.  I use r to do data manipulation during my PhD and I benefit a lot through the discussion in the r-forum. I will continue using R and help the others too. Thank you so much. Regards, Roslina UniSA. [[alternative HTML version deleted]]

Windows XP R 2.8.1 Colleagues, I am trying to run a function testone() and if the function completes without error do one set of instructions, and if the function generates either a warning or an error run another set of instructions. I have read try, and tryCatch help screens at least 20 times, have tried to experiment with code, but I can't understand how to accomplish my desired task.

I have a variable which is roughly age categories in decades. In the original data, it came in coded: > str(xxx) 'data.frame': 58271 obs. of 29 variables: $ issuecat : Factor w/ 5 levels "0 - 39","40 - 49",..: 1 1 1 1... snip I then defined issuecat as ordered: > xxx$issuecat<-as.ordered(xxx$issuecat) When I include issuecat in a glm model, the result

I have just committed this code to CVS, and here is more information about it for those interested : - the spec that I went off of for the protocol can be found here -> http://www.oddsock.org/icecast2yp/spec.html - since icecast2 now needs to potentially connect out to yp servers for communication, libcurl has been added as a dependency...I used curl 7.9.8... - I also have currently

Dear Mrs Rachel Pearce, I am looking for a function "cutpt-coxph" in R - like you did some years ago. How have you solved the problem? Have you found it or a similar function? thank you, Sincerely, Friederike "The title says it all really; I am looking for a function along the lines of cutpt.coxph as described in "Survival Analysis Using S" (Tableman and Kim), Chapter

On Sun, 2 Feb 2003, oddsock wrote: | I have just committed this code to CVS, and here is more information about | it for those interested : This is great! But I have a few wonderings: I'm suprised that the HTTP/GET method of notifying YP servers is still preferred. According to this protocol, we're still limited to essentially title/author, listeners, etc. In other words, we're

At 01:26 PM 2/2/2003 -0500, you wrote: >I'm suprised that the HTTP/GET method of notifying YP servers is still >preferred. According to this protocol, we're still limited to essentially >title/author, listeners, etc. In other words, we're still constrained by >what information is accepted by the YP server. > >I've been thinking that a more XML-like approach could

Pre-fletching on the file and device level has been disabled yielding good results so far. We''ve lowered the number of concurrent ios from 35 to 1 causing the service times to go even lower (1 -> 8ms) but inflating actv (.4 -> 2ms). I''ve followed your recommendation in setting primarycache to metadata. I''ll have to check with our tester in the morning if it made

I''m running into a issue where there seems to be a high number of read iops hitting disks and physical free memory is fluctuating between 200MB -> 450MB out of 16GB total. We have the l2arc configured on a 32GB Intel X25-E ssd and slog on another32GB X25-E ssd. According to our tester, Oracle writes are extremely slow (high latency). Below is a snippet of iostat: r/s w/s

Hi, I have a question about estimating IRT models by using nlme, not just rasch model, but also other models. Behavior Research Methods <http://www.springerlink.com/content/1554-351x/> Volume 37, Number 2 <http://www.springerlink.com/content/1554-351x/37/2/>, 202-218, DOI: 10.3758/BF03192688 Using SAS PROC NLMIXED to fit item response theory models (2005). Ching-Fan