search for: ilovebacon

...sage is in MIME format. The first part should be readable text, while the remaining parts are likely unreadable without MIME-aware tools. --1660387551-1458482416-1249639718=:2997 Content-Type: TEXT/PLAIN; CHARSET=US-ASCII; FORMAT=flowed Content-ID: <Pine.LNX.4.64.0908071039211.2997 at parser.ilovebacon.org> Hello, There appears to be a bug in the nlm function, which I discovered when trying to run an SEM. The SEM is not bizarre; the covariance matrix is solve-able and the eigenvalues are greater than zero, and I do not believe the "sem" package per se to be the problem (as nlm k...

...ed ones) has been that you may be looking at a memory problem on you machine. Given that you can't replicate on another platform (and the .csv file didn't come through), I would think it wise to start there. My 2c. And I love bacon too :) Jeff On Fri, Aug 7, 2009 at 1:10 PM, <adik at ilovebacon.org> wrote: > =A0This message is in MIME format. =A0The first part should be readable t= ext, > =A0while the remaining parts are likely unreadable without MIME-aware too= ls. > > --1660387551-1458482416-1249639718=3D:2997 > Content-Type: TEXT/PLAIN; CHARSET=3DUS-ASCII; FORMAT=3Dfl...

...at a memory > problem on you machine. Given that you can't replicate on another > platform (and the .csv file didn't come through), I would think it > wise to start there. > > My 2c. And I love bacon too :) > > Jeff > > On Fri, Aug 7, 2009 at 1:10 PM, <adik at ilovebacon.org> wrote: >> =A0This message is in MIME format. =A0The first part should be readable = text, >> =A0while the remaining parts are likely unreadable without MIME-aware to= ols. >> >> --1660387551-1458482416-1249639718=3D:2997 >> Content-Type: TEXT/PLAIN; CHARSET=3DU...

Hello, I'm using a very large data set (n > 100,000 for 7 columns), for which I'm pretty happy dealing with pairwise-deleted correlations to populate my correlation table. E.g., a <- cor(cbind(col1, col2, col3),use="pairwise.complete.obs") ...however, I am interested in the number of cases used to compute each cell of the correlation table. I am unable to find such a

Dear R community I have stored the results of arrays in a list consist of J-components (say 200 components). Each component containing same no of columns but may be different no of rows. e.g [[1]] [,1] [,2] [,3] [,4] [,5] [1,] 4 0 0 0 0 [2,] 4 3 4 0 0 [3,] 4 3 4 0 0 [4,] 4 3 0 0 0 [[2]] [,1] [,2] [,3] [,4] [,5]

Dear Colleagues, Occasionally I deal with computer-generated (i.e., websurvey) data files that haven't quite worked correctly. When I try to read the data into R, I get something like this: Error in scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, : line 26 did not have 648 elements ...is there any way to get R to tell me how many elements line 26 *did* have? That

Hello, I have searched and failed for a program or script or method to conduct a power analysis for a MANOVA. My interest is a fairly simple case of 5 dependent variables and a single two-level categorical predictor (though the categories aren't balanced). If anybody happens to know of a script that will do this in R, I'd love to know of it! Otherwise, I'll see about writing one

Dear colleagues, I run R on a few different machines, and view graphs and the like by tunnelling X through SSH to my local machine. This is useful for me because my local machine can't easily handle some of the data sets I work with. However, when an ssh connection dies, the tunnelled X session also dies, which breaks R's device connection, generating this error: > Error: X11

Hi, I'm attempting to compile R 2.7.0 on my G5, which is running OSX 10.5.2. Yes, I know there is a precompiled binary, but my experience has shown that if I hand-compile ATLAS on my machine and then link R against it when compiling, R runs faster, hence me compiling it myself. Anyway, I have this set of errors: checking X11/Intrinsic.h usability... no checking X11/Intrinsic.h presence...

Hello, I need to perform a mixed-model (with nesting) MANCOVA, using Type III sums of squares. I know how to perform each of these types of tests individually, but I am not sure if performing a mixed-model MANCOVA is possible. Please let me know. Erika <>< <>< <>< <>< <>< <>< <>< Erika Crispo, PhD candidate

Hello, I frequently have to export a large quantity of data from some source (for example, a database, or a hand-written perl script) and then read it into R. This occasionally takes a lot of time; I'm usually using read.table("filename",comment.char="",quote="") to read the data once it is written to disk. However, I *know* that the program that generates

Dear Colleagues, Apologies for a long email to ask what I feel may be a very simple question; I figure it's better to overspecify my situation. I was asked a question, recently, by a colleague in my department about pre-aggregating variables, i.e., computing the mean of defined subsets of a data frame. Naturally, I thought of the 'by' and 'tapply' functions, as

Hello, I am trying to use self created functions in other scripts than the one where they are stored. For the moment I am using the following structure of commands to do that: 1. Load the text file with the functions in the current script: x=parse("path") 2. transform the tex in a function: f1=eval(x[1]), f2=eval(x[2]) if more than one function is stored in the text file 3. use the

I have data that looks like this: > df1 group id 1 red A 2 red B 3 red C 4 blue D 5 blue E 6 blue F I want a list of the groups containing vectors with the ids. I am avoiding subset(), as it is only recommended for interactive use. Here's what I have so far: df1 <- data.frame(group=c("red", "red", "red", "blue",

Hi All, I've been having a little trouble with creating a loop that will run a a series of t.tests for inspection, Below is the code i've tried, and some checks i've looked at. I've used the get(paste()) idea as i was told previously that the use of the eval should try and be avoided. I've run a single syntax to check that my systax is correct and works without any problems

Hello, I have an HLM-like situation: a data frame with rows representing individuals, and columns representing some statistics for the individual (e.g., a mean and variance for that individual, plus the number of observations for that person). I'm interested in computing a confidence interval around each individual's mean, but to do that, I would need to know the critical t-value for

Hello, I'm having the following questionable behavior: > summary(m) Min. 1st Qu. Median Mean 3rd Qu. Max. 1 13000 26280 25890 38550 50910 > max(m) [1] 50912 > typeof(m) [1] "integer" > class(m) [1] "integer" ...it seems to me like max() and summary(m)[6] ought to return the same number. Am I doing something wrong? I'm

Hi, I use R in a terminal environment on a linux box, using X11 for graphics, often tunnelled to the terminal I'm using. When an internet connection dies (say, if the ssh connection dies when forwarding), R usually gives a scary warning like "Error: X11 fatal IO error: please save work and shut down R." That's all well and good. Then, if I ignore it and try to plot()

Hello, Just upgraded to 2.5.0, and found that R now includes an rparen (right parentheses) or rbracket whenever I enter in an lparen. While I can see the use of this function, it doesn't mesh well with my personal style of using R (e.g., using the up arrow, adding an rparen, jumping to the beginning of the line, and then wrapping a summary, for instance). Some 10 minutes of google

Hello, I've been making a decent amount of use of the lme function recently, and I am aware that I can extract the variance and correlation components with the VarCorr(model) function. However, I am interested in testing these components as well...is there a module or function available for testing these? I understand there's some debate as to which test is best for a given situation