search for: ririzarr

hi! using: Version 1.8.0 Under development (unstable) (2003-09-02) using the Sweave function of the tools package i get this error when latexing with texi2dvi --pdf ! LaTeX Error: File `upquote.sty' not found. if i change the line \usepackage{/users/faculty/ririzarr/R-1.8.0/lib/R/share/texmf/Sweave} to \usepackage{/users/faculty/ririzarr/R-1.7.1/lib/R/share/texmf/Sweave} in the tex file, then the file compiles fine. this happens to me in both redhat 7.3 and solaris 8.

one more bit of information about this problem. If I start R as the user "postgres" i dont have the segmentation fault. ---------- Forwarded message ---------- Date: Wed, 5 Feb 2003 19:19:39 -0500 (EST) From: Rafael A. Irizarry <ririzarr at jhsph.edu> Reply-To: rafa at jhu.edu To: "R-Help (E-mail)" <r-help at r-project.org> Subject: rdbi segmentation fault hi! i am experiencing the same behaviour explained here: http://finzi.psych.upenn.edu/R/Rhelp02/archive/2482.html (i copied the message below) has anybody...

Using Sweave in the tools library (R version 1.8.0: sorry i havent upgraded), it seems i cant use if statements in R chunks that make graphs. i have this: <<fig=TRUE,echo=F>>= par(mfrow=c(1,1)) if(exists("x")) plot(x,x) else{ plot(1,1,type="n") text(1,1,"data not available.\n") } @ and I get this error: Error: chunk 6 Error in parse(file, n,

im using plclust and want the labels to be different colors. i took a look at getS3method("plot","hclust") and saw a call to .Internal. i looked at the help on .Internal and dont know where to go next. any help appreciated! thanks, rafael

...if help page: "For that reason, one (somewhat extreme) attitude of defensive programming uses braces always, e.g., for 'if' clauses. always enclose expressions in if-else statements even if you dont think you need them." thanks again, -r "Rafael A. Irizarry" <ririzarr at jhsph.edu> writes: > Using Sweave in the tools library (R version 1.8.0: sorry i havent > upgraded), it seems i cant use if statements in R chunks t...

Hi! I have used the Matrix package (Version: 0.995-10) successfully to obtain the OLS solution for a problem where the design matrix X is 44000x6000. X is very sparse (about 80000 non-zeros elements). Now I want to do WLS: (X'WX)^-1X'Wy I tried W=Diagonal(length(w),w) and wX=solve(X,W) but after various minutes R gives a not enough memory error (Im using a 64bit machine with 16Gigs

Hi. guys, I am trying to write my own Stochastic Gradient Ascent for logistic regression in R. But it seems that I am having convergence problem. Am I doing anything wrong, or just the data is off? Here is my code in R - lbw <- read.table("http://www.biostat.jhsph.edu/~ririzarr/Teaching/754/lbw.dat" , header=TRUE) attach(lbw) lbw[1:2,] low age lwt race smoke ptl ht ui ftv bwt 1 0 19 182 2 0 0 0 1 0 2523 2 0 33 155 3 0 0 0 0 3 2551 #-----R implementation of logistic regression : gradient descent ------ sigmoid<-function(z) { 1/(1 + exp(-1*z)) } X<-c...

Ive created a postscript files that looks just like I want it to look except instead of being 6'' x 6'' i want it to be 2''x 2''. i used the postscipt() function (R 1.5.1 on redhat). if i use width and height in par to change the scale the figure looks terrible cause cex,mar,etc.. dont appear to change. before i start tinkering with the cex's and mar... is

Hi, I have a data frame with 6563 observations. I can run a regression with loess using four explanatory variables. If I add a fifth, R crashes. There are no missings in the data, and if I run a regression with any four of the five explanatory variables, it works. Its only when I go from four to five that it crashes. This leads me to believe that it is not an obvious problem with the data,

im using R 1.5.0 on redhat. when i use ccf in the ts library i get, what i think is, an incorrect entry in the lag and acf components. as an example, look at the second entry in the lag and acf components below: > library(ts) > tmp <- ccf(ts(rnorm(3)),ts(rnorm(3)),plot=F) > tmp$lag , , 1 [,1] [1,] -1 [2,] 0 [3,] 0 [4,] 1 > tmp$acf , , 1 [,1] [1,]

if i write a funciton like this f <- function(x) Sweave("report.Rnw") and report.Rnw has an R chunk like this: <<>>= print(x) @ i get an error: 'Object "x" not found'. is there a way to avoid this without using assign to create a global variable? thanks, rafael

Im tryng to understand an error i get with integrate. this is 1.7.0 on solaris 2.8. ##i am trying to approximate an integral of this function, f<-function(b) exp(-(b-mu)^2/(2*tau2))/(p-exp(b))*10^6 ##with tau2 <- .005;mu <- 7.96;p <- 2000 ##from -inf to different upper limits. using integrate(f,-Inf,log(p-exp(1))) ##i get the following error: ##Error in integrate(f, -Inf, log(p -

hi! i want to change the "robustness weights" used by loess. these are described on page 316 of chambers and hastie's "statistical models in S" book as r_i = B(e_i,6m) where B is tukey's biweight function, e_i are the residulas, and m is the median average distance from 0 of the residuals. i want to change 6m to, say, 3m. is there a way to do this? i cant

hi! i am experiencing the same behaviour explained here: http://finzi.psych.upenn.edu/R/Rhelp02/archive/2482.html (i copied the message below) has anybody found a solution? here is the code that gives me a segmentation fault library(Rdbi) library(Rdbi.PgSQL) conn <- dbConnect(PgSQL(), dbname = "PGA") tmp <- "create table test ( expid int, name varchar(128));"