similar to: building formula objects

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.

Hi, I am new to R. Anyone can explain the following from R-help or anyone can direct me how to calculate odds ratio from logistic model in R. Thank you very much. Guoya Stefano <stecalza at tiscalinet.it <https://stat.ethz.ch/mailman/listinfo/r-help> > writes: >Hi all. > >A simple question. >Is there a function to compute the Odds Ratio and its confidence intervall,

I am having trouble getting a loop to work for the following problem. Any help would be much appreciated. Thanks. I need to find the slope and intercept from the linear regression of Drug Level on Day by Participant. There are a total of 37 Participants. I need to store the Participant, Label, Slope, and Intercept in a new data frame. This data is ordered by Participant number 37 total

Hi all. A simple question. Is there a function to compute the Odds Ratio and its confidence intervall, from a logistic model (glm(.......,family=binomial....). I've written my own, but certainly someone did a better job. Thank you in advance, Stefano *********************************************** Stefano Calza Istituto di Statistica Medica e Biometria Universit? degli Studi di Milano Via

Hi, Let's play a little game. Forget all you know, or think you know. Without any context at all, what does the following pattern make you think of? <p> X -> X X -> X | ^ v | X <- X X <- X | ^ v | X X -> X X | ^ | ^ v | v | X -> X

please help me. since 2 days i have tried desperately to find desperately a solution for the following problem: i have an array of 3 different vectors like: > clustertransformationstabelle[,15:17] x1 x2 x3 1 0.6270023 0.2464946 0.5074781 2 0.5629291 0.1783861 0.4398626 3 0.6773455 0.3107641 0.5637289 4 0.0000000 0.0000000 0.0000000 i want to calculate

Hi All, I've hit a problem using nls. I think it may be a restriction in the applicability of nls and I may have found a fix, but I've been wrong before. This example is simplified to the essentials. My real application is much more complicated. Take a function of matrix 'x' with additional arguments: matrix 'aMat' whose values are _not_ to be determined by nls vector

Hello, I have a question regarding the candisc package. My data are: species three five 1 2.95 6.63 1 2.53 7.79 1 3.57 5.65 1 3.16 5.47 2 2.58 4.46 2 2.16 6.22 2 3.27 3.52 I put these in a table and then a linear model >newdata <- lm(cbind(three, five) ~ species, data=rawdata) and then do a candisc on them >candata<-candisc(newdata)

Hello, I'm having some issues getting a ggplot figure to show up in the knitr output, when placed in a loop. Specifically, I have a loop inside a knitr chunk : ```{r fitting, warning=FALSE, fig.width=10, fig.height=10, fig.keep='high'} for (t in 1:T) { # do a regression of tgt.vals ~ predictors and compute coeffs and fitted values (fit.vals / fit.adj.vals) plot(

Am I correct in assuming that the output below essentially translates to "Males have a mean time that is significantly lower than Females"? Is this the correct way to interpret the fact that the coefficient is negative? Assume the variale sex is treated as a factor with Female =0 and Male=1. survmodel<-survreg(survobj~sex,data=data1, dist="weibull")

Hello, I am trying to bootstrap a function that extracts the log-likelihood value and the nls coefficients from an nls object. I want to sample my dataset (pdd) with replacement and for each sampled dataset, I want to run nls and output the nls coefficients and the log-likelihood value. Code: x<-c(1,2,3,4,5,6,7,8,9,10) y<-c(10,11,12,15,19,23,26,28,28,30) pdd<-data.frame(x,y)

Good day! I create diagram of encoder process. Using it i create implementation of encoder in FPGA (Xilinx/Altera). Please critique it. Is there missing stages? Here is blog http://developer-fpga.blogspot.com/ Here is picture of encoding stage 1 https://lh4.googleusercontent.com/-NV8o9DG3jvE/TYjYXr-dYGI/AAAAAAAAAos/U06O-YvhSI0/s1600/stage1.jpg Here is picture of encoding stage 2

I'm a new user, so my apologies for what is likely a dumb question... I am having a hard time getting a table of regression results when using Koenker's code for quantile regression with fixed effects (http://www.econ.uiuc.edu/~roger/research/panel/rq.fit.panel.R). I use the example data parameters that Koenker provides (see below). m <- 3 n <- 10 s <- rep(1:n,rep(m,n)) x <-

Hi Folks, I'm using the 'norm' package (based on Shafer's NORM) on some data. In outline, (X,Y) are bivariate normal, var(X)=0.29, var(Y)=24.4, cov(X,Y)=-0.277, there are some 900 cases, and some 170 values of Y have been set "missing" (NA). The puzzle is that, repeating the multiple imputation starting from the same random seed, I get different answers from the repeats

Hello, list why not add the smart proposal by Greg Snow as a built-in function in {stats}, just changing the "x234" and "newc" lines to allow for more distributions to be generated ? Or do I miss an already existing function to do that ? Regards. Olivier # slight modification of the original code by Greg Snow [mailto:Greg.Snow at intermountainmail.org] # on April 04, 2007

Dear all, The standard call to ARIMA in the base package such as arima(y,c(5,0,0),include.mean=FALSE) gives a full 5th order lag polynomial model with for example coeffs Coefficients: ar1 ar2 ar3 ar4 ar5 0.4715 0.067 -0.1772 0.0256 -0.2550 s.e. 0.1421 0.158 0.1569 0.1602 0.1469 Is it possible (I doubt it but am

Hi again, I have improved precision of my codec (on the encoder and decoder).I have also improved the -h1 quality setting.-Source code and binaries at http://nhwcodec.blogspot.com/-. I am still trying to improve precision of my codec, with keeping my low-complexity (fast) approach. I do not totally use the reference (and impressive) block prediction with different modes + residual coding scheme,

Hi I have gone through the document available in the net regarding the Vorbis encoder /Decoder. Based on that i have prepared a understanding document on the encoder/decoder block. I would like to know whether my understanding of the coder is OK. If there are any other additional block /information pl. provide me with the same. Thanks and regards S.Padmashri <HR NOSHADE> <UL>

Hello R users, I have a puzzle with the VGAM package, on my first excursion into generalized additive models, in that this very nice package seems to want to do either more or less than what I want. Precisely, I have a 4-component outcome, y, and am fitting multinomial logistic regression with one predictor x. What I would like to find out is, is there a single nonlinear function f(x) which acts

Hi Thanks for your response. However I seem to be doing something wrong regarding the try block resulting in yet another error described below. I have a function that takes in a file name and does the fit for the data in that file. Hence based on your input, I tried try ( (fit = lm(y~x, data = data_fitting)), silent = T); I left the subsequent lines of my code unchanged. coeffs =