Displaying 20 results from an estimated 6000 matches similar to: "Finding distance matrix for categorical data"
2010 May 26
3
cluster analysis and supervised classification: an alternative to knn1?
Hi,
I have a 1.000 observations with 10 attributes (of different types: numeric,
dicotomic, categorical ecc..) and a measure M.
I need to cluster these observations in order to assign a new observation
(with the same 10 attributes but not the measure) to a cluster.
I want to calculate for the new observation a measure as the average of the
meausures M of the observations in the cluster
2010 Jun 22
1
Mahalanobis distance
I am a new R user. i have a question about Mahalanobis distance.actually i have 300 rows and 7 columns. columns are different measurements, 300 rows are genes. since genes can
classify into 4 categories. i used dist() with euclidean distance and cmdscale to do MDS plot. but find out Mahalanobis distance may be
better. how do i use Mahalanobis() to generate similar dist object which i can use
2010 Nov 17
2
Bug in agrep computing edit distance?
I posted this yesterday to r-help and Ben Bolker suggested reposting it
here...
Dickison, Daniel <ddickison <at> carnegielearning.com> writes:
>
> The documentation for agrep says it uses the Levenshtein edit distance,
> but it seems to get this wrong in certain cases when there is a
> combination of deletions and substitutions. For example:
>
> >
2009 Aug 23
3
[LLVMdev] RFC: Supporting ELF symbol aliases via GlobalAlias GEPs
Hi Everyone,
Chris suggested[1] I should ask for feedback as to whether this is a
desired feature before I put too much effort into it, so here goes:
I would like to be able to export a symbol that is inside an LLVM
structure. This is possible on ELF targets[2], and the attached proof-
of-concept patch to AsmWriter makes it work (although in a hackish way
that I am NOT suggesting be
2009 Aug 25
0
[LLVMdev] RFC: Supporting ELF symbol aliases via GlobalAlias GEPs
I've attached a less-hackish implementation of this. This includes
the following modifications:
- getSupportsOverlappingAliases() method on TargetMachine which
returns whether the target supports multiple symbols to the same
object. This returns false in the superclass and needs to be
explicitly overridden for each target to enable it.
- An implementation of this method in
2010 Dec 21
2
Warning message when items of Hmisc are masked by loading a package.
I've noticed that I get a warning message every time a package masks
some functions from Hmisc. The warning message says :
Warning message:
In identical(get(., i), get(., lib.pos)) : ignoring non-pairlist attributes
This happens with eg:
library(plyr)
library(xtable)
I think I've seen this passing by before, but I'm not sure any more.
Just thought I'd mention it.
Cheers
Joris
2010 Jun 18
3
Non-procedural access to columns of a matrix
Hi,
I would like to have an index for a column in a matrix encoded in a
cell of the same matrix.
For example:
x = matrix(c(11,12,13,1,
21,22,23,3,
31,32,33,2),byrow=T,ncol=4)
In this case, column 4 is the index. I then access the column
specified in the index by:
> for (i in 1:3) print(x[i,x[i,4]])
[1] 11
[1] 23
[1] 32
>
> for (i in 1:3) {x[i,x[i,4]] <- x[i,x[i,4]] + 5}
> x
2017 Mar 28
2
`[` not recognized as a primitive in certain cases.
?typeof? is your friend here:
> typeof(`[`)
[1] "special"
> typeof(mc[[1]])
[1] "symbol"
> typeof(mc2[[1]])
[1] "special"
so mc[[1]] is a symbol, and thus not a primitive.
- Lukas
> On 28 Mar 2017, at 14:46, Michael Lawrence <lawrence.michael at gene.com> wrote:
>
> There is a difference between the symbol and the function (primitive
>
2009 Aug 25
1
[LLVMdev] RFC: Supporting ELF symbol aliases via GlobalAlias GEPs
Hi David,
Even if this works on Linux/ELF, do know whether it is "officially" supported?
Using aliases to point to the interior of objects seems like something
that could be very likely to break, but I don't know anything about
how ELF encodes aliases.
- Daniel
On Tue, Aug 25, 2009 at 9:07 AM, David Chisnall<csdavec at swansea.ac.uk> wrote:
> I've attached a
2016 Sep 06
2
The use of match.fun
Dear gurus,
I was utterly surprised to learn that one of my examples illustrating the
need of match.fun() doesn't give me the expected result.
center <- function(x,FUN) FUN(x)
center(1:10, mean)
mean <- 4
center(1:10, mean)
Used to give me the error message "could not find function FUN". Now it
just works, even though I didn't expect it to. I believe this is at least
2018 Jan 31
3
Best practices in developing package: From a single file
On 31/01/2018 6:33 AM, Joris Meys wrote:
> 3. given your criticism, I'd like your opinion on where I can improve
> the documentation of https://github.com/CenterForStatistics-UGent/pim.
> I'm currently busy updating the help files for a next release on CRAN,
> so your input is more than welcome.
After this invitation I sent some private comments to Joris. I would
say his
2010 Jun 24
2
count data with a specific range
I would like to prepare the data for barplot. But I only have the data frame
now.
x1=rnorm(10,mean=2)
x2=rnorm(20,mean=-1)
x3=rnorm(15,mean=3)
data=data.frame(x1,x2,x3)
If there a way to put data within a specific range? The expected result is
as follows:
range x1 x2 x3
-10-0 2 5 1 (# points in
this
2010 Jun 08
2
Please help me
Dear Mr. or Ms.,
I used the R-software to run the zero-inflatoin negative binomial model (zeroinfl()) .
Firstly, I introduced one dummy variable to the model as an independent variable, and I got the estimators of parameters. But the results are not satisfied to me. So I introduced three dummy variables to the model. but I could not get the results. And the error message is
2015 Apr 01
4
evaluation in transform versus within
On 01/04/2015 1:35 PM, Gabriel Becker wrote:
> Joris,
>
>
> The second argument to evalq is envir, so that line says, roughly, "call
> environment() to generate me a new environment within the environment
> defined by data".
I think that's not quite right. environment() returns the current
environment, it doesn't create a new one. It is evalq() that created
2017 May 31
4
stats::line() does not produce correct Tukey line when n mod 6 is 2 or 3
Seriously, if a method gives a wrong result, it's wrong. line() does NOT
implement the algorithm of Tukey, even not after the patch. We're not
discussing Excel here, are we?
The method of Tukey is rather clear, and it is NOT using the default
quantile definition from the quantile function. Actually, it doesn't even
use quantiles to define the groups. It just says that the groups
2014 Oct 03
2
How I() works in a formula
Dear all,
I'm updating a package regarding a new type of models, and I'm looking to
extend the formula interface with two functions (L() and R() ) for
construction of these models. I want to use as much of the formula
interface as possible, and hoped to do something similarly to I().
I know the I() function does nothing more than add the class "AsIs". I've
been browsing the
2011 Feb 04
2
terribly annoying bug with POSIXlt : one o'clock is midnight?
Apparently, as.POSIXlt takes one o'clock as the start of the day :
> as.POSIXlt(0,origin="1970-01-01")
[1] "1970-01-01 01:00:00 CET"
> as.POSIXlt(0,origin="1970-01-01 00:00:00")
[1] "1970-01-01 01:00:00 CET"
> as.POSIXlt(0,origin="1970-01-01 23:59:59")
[1] "1970-01-02 00:59:59 CET"
Cheers
--
Joris Meys
Statistical
2010 Mar 30
2
weighted.median function from package R.basic
Dear all,
I want to apply a weighted median on a huge dataset, and I remember a
function from the package R.basic that could do this using an internal
sorting algorithm qsort. This speeded things up quite a bit. Alas, I can't
find that package anywhere anymore. There is a weighted.median function in
the package limma too, but I didn't use that before.
Anybody who knows what happened to
2010 May 25
2
summary of arima model in R
Hi,
I want to give a summary or anova for "arima" model in R, as
"summary", and "anova" for "lm".
As including various intervention factors in arima(xreg = ) part, I
want to assess the significancy of thse factors.
I can do it using interrupted analysis of time series by linear
regression, but want to see whether arima model works for the data
first.
2014 Apr 19
1
lag() not returning a time series object
Dear all,
Before I file this as a bug, I wanted to check if I didn't miss something.
The help page of lag() says that the function returns a time series object.
It actually does return something that looks like a ts object (the
attribute tsp is set). But when using a vector, the class "ts" is not added
to the result:
> avec <- 1:10
> lag(avec)
[1] 1 2 3 4 5 6 7 8