similar to: Problem Matching Exact Values

Dear all, I get a strange error when I find polychoric correlations with the ML method, which I have been able to reproduce using randomly-generated data. What is wrong? I realize that the data that I generated randomly is a bit strange, but it is the only way that I duplicate the error message. > n<-100 > test.x<-rnorm(n, mean=0, sd=1) > test.c<-test.x + rnorm(n, mean=0,

Hello all, Trying to get this piece of code to work on my data set. It is from http://www.itc.nl/personal/rossiter. logit.roc <- function(model, steps=100) { field.name <- attr(attr(terms(formula(model)), "factors"), "dimnames")[[1]][1] eval(parse(text=paste("tmp <- ", ifelse(class(model$data) == "data.frame", "model$data$",

Hello, I am writing out to a file and have two quick questions that I can't seem to track down the correct answers for. Luckily, I *think* they are both simple enough that someone might be able to point me in the right direction on them without too much trouble. Both questions relate to the process below where CompleteFrame is a data frame containing what I want printed to a file. filename

I'm not trying to start a Windows vs. Linux debate, but I've been using R on a Windows machine for a while, and was recently wondering if R's performance would be faster on a Linux machine. And similarly, if any incremental increase in processing speed would be worth the time it would take me to migrate my entire system to Linux (including a database that I access via an R package.)

hi i have gone through the achieves but still could not get my bridge to suvive a reboot . please can anyone help me am using fedora core 2 -----Original message----- From: bridge-request@lists.osdl.org Date: Fri, 28 Jan 2005 03:08:06 +0100 To: bridge@lists.osdl.org Subject: Bridge Digest, Vol 17, Issue 25 > Send Bridge mailing list submissions to > bridge@lists.osdl.org > >

Hi: I have encountered problems with imodwt and universal.thresh.modwt and cannot find any reference in R Search. Hope someone can give me some ideas: Starting with modwt.la8 <- modwt(xdata, "la8", n.level=6) <-- this seems to work fine (1) ydata <- imodwt(modwt.la8) will always give ydata as numeric(0) (no values) instead of being a time series data with

Dear Community, I hope that I have the right category selected because I am relatively new to the "R" world. I come with a relatively challenging problem in the luggage. I would like to realize, that "R" reads text files (there are several hundred pieces in my folder) sequentially, and screens for specific terms. If the term is found, the program should write a 1, if not a

This is hopefully a quick question on decimal accuracy. Is any decimal accuracy lost when casting a numeric vector as a matrix? And then again casting the result back to a numeric? I'm finding that my calculation values are different when I run for loops that manually calculate matrix multiplication as compared to when I cast the vectors as matrices and multiply them using "%*%".

At useR 2006 I mentioned that it would be nice to have a way to specify binomial links that involved free parameters and described some experience with a Gosset link involving a free degrees of freedom parameter, and a Tukey-lambda link with two free parameters. My implementation of this involved some rather kludgey modifications of binomial, make.link and glm that (essentially) added a

Hi, I'm trying to run an R script using Rcmd Batch from the command line on a Windows Vista machine. I am using R version 2.8.1. I installed the batch files 4-3 found at http://cran.r-project.org/contrib/extra/batchfiles/ and added them to my path. I also had to install the latest version of perl (it's Strawberry perl if that makes a difference) and have added this to my path. Now when

Hi, Whenever I issue this commands (as in the help for the threshold function), # # Do a 1D decomposition # tdecomp <- wd(test.data) # # Threshold it # tdecomp.thresh <- threshold(tdecomp) # # Reconstruct from the thresholded coefficients # trecons <- wr(tdecomp.thresh) # #

I fit the model fit<-lm(thresh~cond*Ne) where thresh is the reponse cond is a factor with levels a, b, and c Ne is a continuous indep var I think of this full model as having three lines: thresh as a function of Ne for each condition. Thus we have slopea, slopeb, slopec, inta, intb, intc. lm output my params ------------------------- (Intercept) inta condb intb - inta condc

Airon, I don't think you have to find an English computer 'cause the following must work in your Chinese one :-) Let me explain. First of all, change your lines to xdata <- ckhdat$Adj..Close[1:1447] #names(ckhdwt.la8) <- c("w1", "w2", "w3", "w4", "w5","w6", "v6") note the # sign, i.e., DO NOT change the names

I have a problem with Windows 2000 (SP2) printing to an Epson C60 on a Linux server (running Mandrake 8.1 and Samba 2.2.2 with Cups 1.1.12). I've scoured the net and posts on Usenet have not yielded any helpful response. I have 2 Win2k boxes and they both respond differently, but neither is printing. I've installed the printer on both boxes. On one I get a "unable to

Is there a way to start multiple instances of R in an automated manner? Since I'm not sure that question makes tons of sense, here's my scenario: I have a number of data updates that need to be completed on an ongoing basis with the data pulled from and then stored to another location. The updates are manually triggered, so can build up in terms of the volume that needs to be updated.

I am still having some difficulties with connecting R to a C++ function. I am able to call the function as expected after compiling the shared library and such. However, the call to the function is via .C; parameters from the .C call are not being passed correctly to the function. As an example, I have attached a GDB run of the code. I set a breakpoint on entry to the function I am calling

I'm working on a problem where I'm introducing random error and have been using the built in function runif to provide that random error. However, I realized that I seem to be getting some unexpected behavior out of the function and was hoping someone could share some insight. I don't know the runif algorithm at all, but from the behavior I'm seeing, it seems that whenever I open

Hi All, I'm generating 10 different data sets with 1 and 0 in a matrix form and writing the output in separate files. Now I need to stack all these data sets in one vector and I know that stack only operates on list or data frame however I got these data sets by converting list to a matrix so can't go backwards now. Is there a way i can still use Stack? Please see the program:

Dear all, I'm trying to use "scatterplot" from package "car" to get a scatterplot with row-names as labels and I want my data to be divided into two groups defined by the column "tresh". Leaving away the groups-command gives me the right labels but they are not distinguished by colors according to "tresh". Does anybody know why the following example

Hi, I'm trying to use mars function in R to interpolate nonlinear multivariate functions. However, it seems that mars gives me a fit which uses only very few basis function and it underfits very badly. For example, I have tried the following code to test mars: require("mda") f <- function(x,y) { x^2-y^2 }; #f <- function(x,y) { x+2*y }; # Grid x <-