similar to: imporving performance of slicing on matrices and S4 their derivatives

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: