search for: inman

...", i, sep="") <- input.value[i] } Of course this loop will not work since the paste function returns a character string that cannot be a variable. How can I construct a loop to create sequential variables and assign a value to them in the same spirit of the loop above? Brant Inman Mayo Clinic [[alternative HTML version deleted]]

...milar to Faraway (Regression), Pinheiro&Bates (Mixed Effects) and Wood (GAMs) and would be very modern (i.e. published within the last couple of years). Any ideas? If such a book does not exist, one of the experts on this help list should write it! (I will gladly buy the first copy.) Brant Inman Mayo Clinic [[alternative HTML version deleted]]

...#39;l' but I don't want to, I find that the plots look less cluttered with disjoint axes. However, I was intrigued enough by her comments that I decided to solicit the opinions of others on this issue. Are there principled reasons why one should prefer joined axes or disjoint axes? Brant Inman

...ollowing website: http://gnuwin32.sourceforge.net/packages/tiff.htm I would like to ask the R experts for help with the following things: 1) Where do I get the tiffio.h file? 2) Where do I install or relocate the tiffio.h and TIFF library to so that rtiff will work? Thanks for your help. Brant Inman Mayo Clinic

R-helpers: I would like to measure the correlation coefficient between the repeated measures of a single variable that is measured over time and is unbalanced. As an example, consider the Orthodont dataset from package nlme, where the model is: fit <- lmer(distance ~ age + (1 | Subject), data=Orthodont) I would like to measure the correlation b/t the variable "distance" at

...use the symbol 'createLLVMOptimizer' (or anything with 'lto', etc). The document at http://llvm.org/releases/2.1/docs/LinkTimeOptimization.html doesn't seem to be of much help for getting started with this. How do I provoke LTO? Any suggestions? Thanks, Jeff Jeff Inman -- jti at lanl.gov Advanced Architectures for Computation Los Alamos National Laboratory

...- read.S(file='H:\\Research\\data.sdd') Error in read.S(file = "H:\\Research\\data.sdd") : not an S object The dataset is fairly large, roughly 13000 rows and 100 columns, but within S-Plus 7.0 it is stored as a normal dataframe (not a bd dataframe). Any ideas? Brant Inman [[alternative HTML version deleted]]

...- read.S(file='H:\\Research\\data.sdd') Error in read.S(file = "H:\\Research\\data.sdd") : not an S object The dataset is fairly large, roughly 13000 rows and 100 columns, but within S-Plus 7.0 it is stored as a normal dataframe (not a bd dataframe). Any ideas? Brant Inman [[alternative HTML version deleted]]

...rcepts: mild|normal normal|severe -1.630612 2.563055 Residual Deviance: 503.7251 AIC: 509.7251 ------ Can someone explain what I am doing wrong when using vglm and polr with the modified dataset? I do not understand why these two formulations should give different results. Brant Inman

...plies to his posting how he managed to solve the problem, so I am reposting a similar note again in hopes of guidance. In particular, I would like to know if and how I need to modify the confint function call so that I can access it from within another function that I have written. Thanks, Brant Inman

R-Experts: Does anyone know if there are R functions to perform the Freeman-Tukey double arcsine transformation and then backtransform it? Thanks, Brant Inman Mayo Clinic

...*binary* package ‘pcvsuite’ ... * DONE (pcvsuite) > library(pcvsuite) Error: package 'pcvsuite' was built before R 2.10.0: please re-install it > Anybody have an idea of how I can use the pcvsuite package running on R 2.12.0 for Mac? Thank you in advance for any solutions. Brant Inman Duke University [[alternative HTML version deleted]]

...ote the problems in rows 15, 21, 37 and 43. They should read [0.5, 45.5], [14.5, 21.5], etc... The matrix function seems to be rounding the second column up to the next integer. Why would this occur? Can I do something to prevent this? I would very much appreciate any comments. Thanks, Brant Inman Mayo Clinic [[alternative HTML version deleted]]

...rm is the size of the usnews dataset or the amount of missing data. -------------------- > dim(usnews) [1] 1302 7 -------------------- Does anyone have any ideas? If not, are there alternatives to norm for implementing the MLE and EM methods of dealing with missing data? Thanks, Brant Inman Mayo Clinic

...ain site (http://bugs.r-project.org/cgi-bin/R) , I get an error page with the following message: The system encountered a fatal error cannot open config file /home/sfe/r-bugs/jitterbug/R : No such file or directory The last error code was: No such file or directory uid/gid=30/8 Thanks, Brant Inman

...There appears to be a problem with the predict.drc function. This code previous generated confidence and prediction intervals in columns 2 and 3 of the prediction matrices but now fails for reasons unknown to me. Anyone have an idea of what is going on and how I can remedy the situation? Brant Inman [[alternative HTML version deleted]]

...ror in Surv(time, 1 - cens) : Time variable is not numeric In addition: Warning message: is.na() applied to non-(list or vector) in: is.na(time) 2) Can the linear predictor (linp) be used in sbrier in the same way as a probability might? 3) Has anyone implemented Royston/Sauerbrei's D? Brant Inman Mayo Clinic [[alternative HTML version deleted]]

...two.sided') } } Error in result[i, j] <- power.t.test(delta = DELTA[i], sd = SD[j], sig.level = 0.05, : number of items to replace is not a multiple of replacement length ----------------- Can some one tell me what I am doing wrong here? Thanks in advance for your help, Brant Inman

...0.016, which suggests that x1 is indeed a significant predictor of y. How do I reconcile these two differing messages? I do recognize that the upper bound of the confidence interval is pretty close to 1, but I am certain that some journal reviewer will point out the problem as inconsistent. Brant Inman

R plotting experts: I have a bivariate dataset composed of 300 (x,y) continuous datapoints. 297 of these points are located within the y range of [0,10], while 2 are located at 20 and one at 55. No coding errors, real outliers. When plotting these data with a scatterplot, I obviously have a problem. If I plot the full dataset with ylim = c(0,55), then I cannot see the structure in the data in