Displaying 20 results from an estimated 600 matches similar to: "imporving performance of slicing on matrices and S4 their derivatives"
2014 Mar 31
1
C API to get numrow of data frame
Hi ,
Is there any C API to the R API nrow of dataframe ?
x<- data.frame()
n<- nrow(x)
print(n)
0
Example :
My C function which deals with data frame looks like and I don't to send
the number of rows of data frame .I want to detect it from the function
itself, my function take data frame as argument and do some on it. I want
API equivalent to nrow. I tried Rf_nrows,Rf_ncols . No
2003 Jun 25
3
joining columns as in a relational database
In our recent workshop on "Multilevel Modeling in R" we discussed
handling data for multilevel modeling. An classic example of such
data are test scores of students grouped into schools. We may wish to
model the scores as functions of both student-level covariates and
school-level covariates.
Such data are often organized in a multi-table format with a separate table
for each level of
2009 Dec 16
2
Flexclust barchart issue when mcol=NULL (PR#14150)
Full_Name: Chris Hane
Version: 2.10.1
OS: Windows
Submission from: (NULL) (198.203.181.181)
When using barchart in the flexcust package, setting mcol=NULL to avoid the
lollipops causes an error. Each panel shows the text message "Error using
packet n replacement has length zero." where n is the panel number.
> data(iris)
> cl <- cclust(iris[,-5], k=3)
> barplot(cl,
2023 Jul 08
1
Getting an error calling MASS::boxcox in a function
Dear Ron and Bert,
First (and without considering why one would want to do this, e.g.,
adding a start of 1 to the data), the following works for me:
------ snip ------
> library(MASS)
> BoxCoxLambda <- function(z){
+ b <- boxcox(z + 1 ~ 1,
+ lambda = seq(-5, 5, length.out = 101),
+ plotit = FALSE)
+ b$x[which.max(b$y)]
+ }
> mrow <- 500
2023 Jul 08
1
Getting an error calling MASS::boxcox in a function
Thanks John.
?boxcox says:
*************************
Arguments
object
a formula or fitted model object. Currently only lm and aov objects are handled.
*************************
I read that as saying that
boxcox(lm(z+1 ~ 1),...)
should run without error. But it didn't. And perhaps here's why:
BoxCoxLambda <- function(z){
b <- MASS:::boxcox.lm(lm(z+1 ~ 1), lambda = seq(-5, 5,
2013 Mar 28
1
scatterplot3d with densCols ?
Hi,
I was trying to make a 3D plot using densCols. The documentation for
densCols doesn't look like it'll work for 3D. For example:
-----------------------------------------
library(scatterplot3d)
v1 <- rnorm(10000)
v2 <- rnorm(10000)
v3 <- rnorm(10000)
## 2D with denscols
mat1 <- cbind(v1,v2)
mcols1 <- densCols(mat1)
plot(mat1,col=mcols1)
mat <- cbind(v1,v2,v3)
2023 Jul 08
1
Getting an error calling MASS::boxcox in a function
Hi Bert,
On 2023-07-08 3:42 p.m., Bert Gunter wrote:
> Caution: This email may have originated from outside the organization. Please exercise additional caution with any links and attachments.
>
>
> Thanks John.
>
> ?boxcox says:
>
> *************************
> Arguments
>
> object
>
> a formula or fitted model object. Currently only lm and aov objects
2023 Jul 08
1
Getting an error calling MASS::boxcox in a function
Hi,
Firstly, apologies as I have posted this on community.rstudio.com too.
I want to optimise a Box-Cox transformation on columns of a matrix (ie, a unique lambda for each column). So I wrote a function that includes the call to MASS::boxcox in order that it can be applied to each column easily. Except that I'm getting an error when calling the function. If I just extract a column of the
2023 Jul 08
1
Getting an error calling MASS::boxcox in a function
Try this for your function:
BoxCoxLambda <- function(z){
y <- z
b <- boxcox(y + 1 ~ 1,lambda = seq(-5, 5, length.out = 61), plotit =
FALSE)
b$x[which.max(b$y)] # best lambda
}
***I think*** (corrections and clarification strongly welcomed!) that `~`
(the formula function) is looking for 'z' in the GlobalEnv, the caller of
apply(), and not finding it. It finds
2023 Jul 08
1
Getting an error calling MASS::boxcox in a function
No, I'm afraid I'm wrong. Something went wrong with my R session and gave
me incorrect answers. After restarting, I continued to get the same error
as you did with my supposed "fix." So just ignore what I said and sorry for
the noise.
-- Bert
On Sat, Jul 8, 2023 at 8:28?AM Bert Gunter <bgunter.4567 at gmail.com> wrote:
> Try this for your function:
>
>
2012 Sep 28
1
Heatmap Colors
Hello R-Users!
I'm using a heatmap to visualize a matrix of values between -1 and 3.
How can I set the colors so that white is zero, below zero is blue of
increasing intensity towards -1 and above zero is red of increasing
intensity towards red?
I tried like this (using the marray and gplots packages from bioconductor):
mcol <- maPalette(low="blue", mid="white",
2011 May 15
1
pls help: lattice graph with both negative and positive value, x and y cross at 0 and negative value bars are plotted just oppositive direction in contrast to positive
Dear R experts:
Here is my problem:
#Data 1
Y <- c(0.5, 0.1, 0.5, 1.3, 1.4, 1.6, 1.65, 2.4, 2.6, 3.4, 3.6, 4.3, 4.42,
4.8, 4.7, 3.4, 3.3, 2.8, 2.8, 1.2, 1.1, 0.5, 0.2, 0.1, -0.2, -1.5, -2.5,
-1.3, -0.5, -0.1)
X <- seq(1:30)
X1 <- c(rep("T1", 24), rep("T2", 6))
dat1 <- data.frame(Y, X, X1)
require(lattice)
mcol <- c("green", "red")
2013 Mar 06
1
red SNA
Estimados
Le consulto porque estoy intentando presentar unos datos de alguna forma un
poco distinta, salir un poco de la costumbre.
Primero describo los datos:
Suponiendo una secuencia:
golpe al auto -> torura de vidrio -> rotura de aciento
otro auto
golpe al auto -> -> rotura de aciento , a este por alguna causa no se le
rompio el vidio o esa información es desconocida.
Esta
2003 May 21
2
moving onto returning a data.frame?
I've been studying some of the code and I'm still a little shakey on the
proper method for returning a data.frame from a C function (which is my
ultimate goal here). I've started some code that I've "stolen" from the
archives and I'm running into crashes, etc. I've been trying to gleen some
insight from the src/main/scan.c file and didn't find many comments in
2013 Nov 22
2
[LLVMdev] new dragonegg 3.4 warning on darwin12
Duncan,
Is this expected for x86_64-apple-darwin12 under Xcode 5? Using current
llvm/compiler-rt/clang/polly/testsuite with dragonegg 3.4 branch, I see
the warnings...
% /sw/lib/gcc4.8/bin/gcc-4 -fplugin=/sw/lib/gcc4.8/lib/dragonegg.so -specs=/sw/lib/gcc4.8/lib/integrated-as.specs -fplugin-arg-dragonegg-enable-gcc-optzns -Ofast himenoBMTxpa.c
himenoBMTxpa.c: In function ‘main’:
2013 Nov 22
0
[LLVMdev] new dragonegg 3.4 warning on darwin12
On Thu, Nov 21, 2013 at 07:37:46PM -0500, Jack Howarth wrote:
> Duncan,
> Is this expected for x86_64-apple-darwin12 under Xcode 5? Using current
> llvm/compiler-rt/clang/polly/testsuite with dragonegg 3.4 branch, I see
> the warnings...
>
> % /sw/lib/gcc4.8/bin/gcc-4 -fplugin=/sw/lib/gcc4.8/lib/dragonegg.so -specs=/sw/lib/gcc4.8/lib/integrated-as.specs
1998 Nov 28
1
No subject
Dear Friends,
Yesterday I posed a question to the list concerning the possibility of doing
animation examples in R. Here is an example S-Plus (4.5 for Windows ) that
I wrote to illustrate my problem. If I try this in R (comment out the
graphsheet and guilocator calls) I don't see my results until after the
function has iterated through the 100 trials. I would like to plot each
iteration.
2020 Sep 12
0
Garbage collection of seemingly PROTECTed pairlist
Hi,
In line 5, you are allocating a vector of length nc. Then, in line 12, you
are using nr as a limit, so if nr goes beyond nc, which is happening in
line 39, you are in trouble.
I?aki
On Sat, 12 Sep 2020 at 03:30, Rory Nolan <rorynoolan at gmail.com> wrote:
> I want to write an R function using R's C interface that takes a 2-column
> matrix of increasing, non-overlapping
2017 Apr 28
1
pairwiseAlignment Improvements
Good day,
The location of indels can be retrieved from a PairwiseAlignmentsSingleSubject object by using indel. Determining any difference between the two sequences, including substitutions, is not quick nor easy. I suppose that summary displays details of the mismatches, but the variable is of class PairwiseAlignmentsSingleSubjectSummary which has no documented accessors. So, the code to access
2008 Aug 15
2
[LLVMdev] Interprocedural Program Slicing
Dear all,
Is anyone aware of any implementation of an inter procedural slicing
algorithm
for LLVM? I am mostly interested in a static algorithm but a dynamic one can
be helpful too.
Thanks,
Silviu Ganceanu
Research Student,
Ecole Polytechique Federale de Luasanne
-------------- next part --------------
An HTML attachment was scrubbed...
URL: