similar to: Matrix Integral

Hello, I was curious if there was a complex valued matrix exponential function available for R? I have some Laplace transforms of occupation times for a hidden Markov model. The matrix exponential function in the msm package does not seem to handle complex values. For example > MatrixExp(diag(1i,2)) [,1] [,2] [1,] 1 0 [2,] 0 1 Warning message: imaginary parts

ABCD are four matrix. A * Inverse((Transpose(A)*Tranpose(B)*B*A+C)) * Transpose(A) * Transpose(B) * D how to write in R in an efficient way? --------------------------------- [[alternative HTML version deleted]]

Hi guys, I need your help. I would like to select a subset from a dataset. This is the dimension of the dataset. > dim(data1) [1] 72 36 1916 so, it's like.. there are 1916 of 72 * 36 matrix. ==> looks like 72 * ( 36*1916 ) ** *1)* And I would like to select the first 72*36 matrix. This is how I did: > two=data1[,1:36] Error in data1[, 1:36] : incorrect number of dimensions

Hi, Not sure if this is the best section to post it in. I'm doing lengthy runs on a 30GB data set - I would like to know if it is possible through any means for R to send an email / text message or anything when it has finished running a source file. Any suggestions would be appreciated. Nick. -- View this message in context:

Is the matrix exponential available in some package? The cannonical reference is "Nineteen dubious ways to take the exponential of a matrix". (Love that title) Terry T.

Guys, I've downloaded AsteriskNOW few days ago so I'm new to this product. The first issue is on service provider area. I've already used a VoIP account already configured with my ISP, it works fine! This configuration has been used until now with the client SJphone, Now I would use this profile as main VoIP service provider to setup in AsteriskNOW. Here are the profile detail as

http://idisk.mac.com/jdeleeuw-Public/utilities/apl/apl.R Dedicated to the IBM 2741. Implemented for general multidimensional arrays: drop, take, reshape, shape, rank, select, generalized inner product, generalized outer product, representation, base value, join, expand, reduce, scan, member of, ravel, compress, tranpose, rotate Basically, the APL-I part is complete, and after some testing

I have code and documentation that don't match, but R CMD check didn't flag it. in mspath.R mspath <- function(formula, # formula with observed Markov states ~ observation times (required) qmatrix, # matrix of 1s and 0s with indices of allowed transitions (diagonal is ignored) (required) misc = FALSE, ematrix = NULL, # matrix

x <- read.table("LittleGarvin.csv", sep=",", header=TRUE) y <- t(x) str(y) chr [1:193, 1:288] "oligocha" "0" " 0" " 0" " 0" "0" ... - attr(*, "dimnames")=List of 2 ..$ : chr [1:193] "X" "upwd1201" "upwd0502" "upwd0702" ... ..$ : NULL x is a data frame with

> -----Original Message----- > From: Federico Calboli > Sent: Tuesday, January 20, 2004 5:40 PM > To: r-help > Subject: [R] matrix exponential: M^0 > > I would like to ask why the zeroeth power of a matrix gives me a matrix > of ones rather than the identity matrix: > > > D<-rbind(c(0,0,0),c(0,0,0),c(0,0,0)) > > D<-as.matrix(D) > > D >

I have an array with 44800 columns and 24 rows I would like to compute the row average for the array 100 columns at a time, so I would like to end up with an array of 24 rows x 448 columns. I have tried using apply(dataset, 1, function(x) mean(x[])), but I am not sure how to get it to take the average 100 columns at a time. Any ideas would be welcomed. thanks, Spencer [[alternative HTML

Hello best helpers, I am a new user and I have been struggling for hours. So finally I decide to ask you: If I have a matrix P, and P.2 = P%*%P, and P.3=P.2%*%P is there a function to calculate the power of a matrix?? if not how can i do: for (i in 1:10) {P.i=P^i} after this I need to sum them up and my problem is to combine P and i to P.i can anyone help me please??? Thanks and have a nice

hi, everyone: i have a question on the reshape function. i have the following dataset : gene tissue patient1 patient2 patient3............. _________________________________________________ gene1 breast 10 20 50 gene2 breast 20 40 60 gene3 breast 100 200 300 which i hope to convert to the following format: gene patientID

Hi, Is there a function for raising a matrix to a power? For example if you like to compute A%*%A%*%A, is there an abbreviation similar to A^3? Atte Tenkanen > A=rbind(c(1,1),c(-1,-2)) > A [,1] [,2] [1,] 1 1 [2,] -1 -2 > A^3 [,1] [,2] [1,] 1 1 [2,] -1 -8 But: > A%*%A%*%A [,1] [,2] [1,] 1 2 [2,] -2 -5

Dear R-Users, Suppose I have data in the following format: CODE_NAME ZIP_CODE John 12345 John 23456 John 34567 Jane 13242 Jane 22123 I want to transpose / convert it into: CODE_NAME ZIP_CODE John 12345,23456,34567 Jane 13242,22123 Any idea/pointer is appreciated.

Hi, I need to transpose and aggregate the 4th and 5th column of this data set (below). The period is 1 to 100 years and some period has 1 entry as shown in the example below. Any ideas how to do this. Thanks in advance. Noli ~~~~~~~~~~ From: P_ID CROP Period Ini_Age Area_Cut 83 SORI 1 31 528.2465512 84 SORI 1 32 74.55179899 85 SORI 1 33

I have been looking at some 'table' examples in Peter Dalgaard's ISwR book, and I am confused by how to get right justification of my table headings when I use the tables() command. Compare the following: # Produces right justfified column names caff.marital=matrix(c(652,1537,598,242,36,46,38,21,218,327,106,67),nrow=3,byrow=T)

I am trying to solve an integral in R. However, I am getting an error when I am trying to solve for that integral. The equation that I am trying to solve is as follows: $$ C_m = \frac{{abs{x}}e^{2x}}{\pi^{1/2}}\int_0^t t^{-3/2}e^{-x^2/t-t}dt $$ [image: enter image description here] The code that I am using is as follows: a <- seq(from=-10, by=0.5,length=100) ## Create a function to compute

Dear r-users, Is there any command in R allowing to evaluate a double integral? for instance let say I want to evaluate the following integral: integrate[lo=(0,1),up=(2,3)] f(x,y)=x^2+y^2 where lo is the vector of lower bounds and up that of upper bounds. I thaught the function "adapt" would work but it did not. Many thanks, Dominique K.

a=[0.1,0.2,0.1,0.3,0.4] b=[0.2,0.3,0.1,0.2,0.5] c=[1,1,1,1,1] log(c+a-x*b) where x=unknown scale variable. int=$$log(c+a-x*b)dadb, where $ denotes integral sign. Actually, how could I calculate the integral's approximation? double summation? best, moohwan