[R] One-dimensional point processes

2007-06-05 Thread dhinds
I'm interested in characterizing data from a one-dimensional point process... say, testing for stationarity, or measuring the degree of clustering at different length scales. This sounds like a spatial statistics problem but the available tools all seem focused on 2+ dimensional data. Time

[R] Goodness of fit for hclust?

2007-05-21 Thread dhinds
I'd like to get a measure of goodness of fit for a heirarchical clustering result from hclust. Something that would indicate the extent to which the dendrogram accurately represents the original dissimilarity matrix. Is there an easy way to do this? Or, does anyone have code for computing

Re: [R] Object attributes in R

2006-10-12 Thread dhinds
Michael Toews [EMAIL PROTECTED] wrote: Is there any way of keeping the attributes when subsetted from primitive classes, like a fictional attr.drop option within the [ braces? The best alternative I have found is to make a new object, and copy the attributes: tm2 - tm[3:5]

Re: [R] Object attributes in R

2006-10-12 Thread dhinds
[EMAIL PROTECTED] wrote: What I've done is to define a subclass that keeps attributes, that can be added to any object, shown below. The keep.attr() function is supposed to return just user attributes but I'm not sure if my list of special ones is complete. Sorry, the following version

Re: [R] Numerical print format and loading with RMySQL

2006-06-13 Thread dhinds
Nathan Johnson [EMAIL PROTECTED] wrote: I constructed a data frame with columns corresponding to the table fields, and tried with and without a null column for the auto- increment internal id column. I also tried renaming all the column names to match those of the table fields. All of

Re: [R] Fastest way to do HWE.exact test on 100K SNP data?

2006-06-06 Thread dhinds
Anna Pluzhnikov [EMAIL PROTECTED] wrote: Hi everyone, I'm using the function 'HWE.exact' of 'genetics' package to compute p-values of the HWE test. My data set consists of ~600 subjects (cases and controls) typed at ~ 10K SNP markers; the test is applied separately to cases and controls.

Re: [R] Reshaping genetic data from long to wide

2006-04-12 Thread dhinds
Farrel Buchinsky [EMAIL PROTECTED] wrote: 2)Storing all the SNP data as a string seems quite clever and a space-saving way of doing it. However, if you were to analyze a whole chromosome at a time you would still be creating one almighty big table albeit only temporarily. Do you use R to run

Re: [R] Reshaping genetic data from long to wide

2006-04-08 Thread dhinds
Farrel Buchinsky [EMAIL PROTECTED] wrote: Bottom Line Up Front: How does one reshape genetic data from long to wide? I avoid both your long and wide formats because they are awkward and inefficient for large data sets. The long format wastes a huge amount of space with redundant column values,

Re: [R] hypothesis testing for rank-deficient linear models

2006-01-12 Thread dhinds
[EMAIL PROTECTED] wrote: Take the following example: a - rnorm(100) b - trunc(3*runif(100)) g - factor(trunc(4*runif(100)),labels=c('A','B','C','D')) y - rnorm(100) + a + (b+1) * (unclass(g)+2) ... Here's a cleaned-up function to compute estimable within-group effects for

[R] hypothesis testing for rank-deficient linear models

2006-01-10 Thread dhinds
Take the following example: a - rnorm(100) b - trunc(3*runif(100)) g - factor(trunc(4*runif(100)),labels=c('A','B','C','D')) y - rnorm(100) + a + (b+1) * (unclass(g)+2) m - lm(y~a+b*g) summary(m) Here b is discrete but not treated as a factor. I am interested in

Re: [R] Problem installing ROracle package under R

2005-11-04 Thread dhinds
In gmane.comp.lang.r.general Vittorio [EMAIL PROTECTED] wrote: configuration warning: Oracle pre-compiler proc not in /usr/local/oracle8-client/bin/proc you may not be able to compile ROracle You don't have the Oracle Pro*C precompiler installed. You can't build ROracle without that.

Re: [R] problem with the assignment function

2005-11-04 Thread dhinds
Wang, Yan [EMAIL PROTECTED] wrote: I run into the most weird problem I have ever met. I wrote a function rhopair, which calls a .C function. I cannot assign its value to a variable using either = nor -. After I did the assignment, rhopair cannot reproduce the same result as before with the

Re: [R] Problem with memory footprint of qq plot generated with lattice

2005-09-29 Thread dhinds
Witold Eryk Wolski [EMAIL PROTECTED] wrote: Adding f.value=fn as argument to qqmath reduces the size of the image, but neither the axis (absicissae) nor the line added by panel.qqmathline are right. Adding f.value=fn as argument to panel.qqmathline and panel.qqmath generates the right

Re: [R] Problem with memory footprint of qq plot generated with lattice

2005-09-28 Thread dhinds
nwew [EMAIL PROTECTED] wrote: Dear R helpers, I generate a qq plot using the following function call. ... dim(xx) [1] 680237 2 How about doing something like this: fn - function(n,cut=0.001,m=1000) { p - ppoints(n) p - p[pmin(p, 1-p) cut] q -

Re: [R] ROracle install problem

2005-09-08 Thread dhinds
Ariel Chernomoretz [EMAIL PROTECTED] wrote: Hi, I am trying to install the ROracle package in a Linux-64 machine. I downloaded from Oracle's site their Instant Client bundle but it seems that ROracle needs some stuff not included in that kit in order to compile (in particuar, the 'proc'

Re: [R] ROracle and select query empty

2005-08-29 Thread dhinds
Mathieu Drapeau [EMAIL PROTECTED] wrote: I am using Linux Oracle Client 9i and I am running my R scripts on a Linux box. Does /usr/include/sqlca.h exist? This is a Postgres file. My hazy memory is that this conflicts with an Oracle header; and that if ROracle sees this one at compile time,

Re: [R] ROracle and select query empty

2005-08-29 Thread dhinds
Mathieu Drapeau [EMAIL PROTECTED] wrote: Thanks!!! It is working now. I'll look into whether ROracle can be tweaked to prevent this problem. -- Dave __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] ROracle and select query empty

2005-08-27 Thread dhinds
Mathieu Drapeau [EMAIL PROTECTED] wrote: Hi, I just installed ROracle and RDBI. The connection to the database seems to work also. My problem is when I am selection rows that really exist in the database, it is returning nothing. Where should I look to see what could be my problem? What

[R] Raw data type transformations

2005-08-17 Thread dhinds
Not sure if this belongs here or on the devel list: I've needed a more efficient way to manipulate raw binary data in R, with more than the minimal raw transformation functions in the base package. So I've written a small package in C that effectively lets me cast back and forth between raw

[R] New trellis settings

2004-11-29 Thread dhinds
I've just upgraded to 2.0.1 and was taken by surprise by the changes in graphical parameter handling for lattice plots. I'd previously been using 1.9.1. The old settings seem to have been replaced by a daunting number of new options. I've poked around a bit and have not seen any discussion of