similar to: Heatmap Colors

I have a matrix of, say, M and N dimensions: my_matrix=matrix(c(1:60),nrow=6,ncol=10) I have two "id" vectors corresponding to the rows and columns, e.g.: id_m=seq(10,60,by=10) id_n=seq(100,1000,by=100) I would like to create a "proper" database (let's say a data.frame for this example -- i'm going to be loading these into an SQLite database, but we'll leave

I'm working on some heatmaps, and the person I'm working with would prefer a red-black-green color palette (red denoting gene induction and green denoting gene repression). Does such a palette exist already? If not, is there an easy way to create one? Thanks, Jake

Hi R-users I can not change the name of one column only of my matrix. my_matrix <- matrix (1:12,ncol=3) colnames(my_matrix)[1] <- 'myname' Error in dimnames(x) <- dn : length of 'dimnames' [2] not equal to array extent thank you for your help Lorenzo [[alternative HTML version deleted]]

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,

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

Hello, I saw this nice trick I want to replicate but I lost the source and I hope one of you can point me to the solution. My problem is that I don't know the correct words to query this. When I print to screen a matrix or data.frame the columns are split and printed below the previous ones; even though I have plenty of screen left. E.g., > my_matrix = matrix(runif(30),nrow=3,ncol=10)

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,

Hello, I have 2 functions (a and b) a = function(n) { matrix (runif(n*2,0.0,1), n) } > > > b = function (m, matrix) { > n=nrow (matrix) > p=ceiling (n/m) > lapply (1:p, function (l,n,m) { > inf = ((l-1)*m)+1 > if (l<p) sup=((l-1)*m)+m > else sup=n >

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: > >

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

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

just reinstalled @home and i have a one of those 100 cards, anyways when i call from the pstn the box picks up but i hear nothing, then it clicks a couple times, then nothing again, i am trying to get the digital receptionist to work but it won't save my wav file to the @home box and all the radio buttons under incoming calls are greyed out. the greyed out thing seems to be my biggest problem

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

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")

Hello Everyone, I am using R for creating heatmap from a square matrix. Below is my script to do so my_map <- read.csv("Desktop/input.csv", sep=",") my_matrix <- data.matrix(my_map) my_heatmap <- heatmap(my_matrix, Rowv=NA, Colv=NA, col = cm.colors(256), scale="column", margins=c(5,10)) I get a beautiful plot as output, but I could not save this plot

Does anyone know of any sip wifi phones? Only one i can find that is redily availiable is the zyxel prestige 2000w and from what i hear it is flaky. Kurt Fankhauser WaveLinc HYPERLINK "http://www.wavelinc.com/"www.wavelinc.com 114 S. Walnut St. Bucyrus, OH 44820 419-562-6405 -- No virus found in this outgoing message. Checked by AVG Anti-Virus. Version: 7.0.300 / Virus Database:

Hi, I have a comma separated file with element names in first column like shown below : Name_1,0 Name_2,0.8878,0 Name_3,0.6777,0.7643,0 Name_4,0.9844,0.1234,0.1414,0 Original data is a 10000x10000 symmetric matrix (600 MB). To reduce file size, I have minimized matrix to only lower triangle. Is there a (memory) efficient way to 1) read file 2) compute first and second principal components and

Hi, I just joined the list, anyways i am trying to setup an @home box with a x100p card and so far i can't even get the box to pickup the incoming call and in the amp management under the section "send calls from PSTN too" page all the radio buttons are blank and i want to use the digital receptionist, also when i try to setup digital receptionist via uploading wav file and save, it

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

I'm using functions to return a matrix of all permutations of a specified size from a larger list for predictor selection. For each predictor size I use a seperate function like this: bag2 <- function(n) { rl <- c() for (i1 in seq(n)) { for (i2 in seq(n)) { if (length(unique(c(i1,i2)))==1) {next} rl <- cbind(rl,matrix(c(i1,i2))) } }