Displaying 7 results from an estimated 7 matches for "imatrix".
Did you mean:
matrix
2011 Jan 27
3
agnes clustering and NAs
...onents:
[1] "order" "height" "ac" "merge" "diss" "call" "method" "data"
But I have a large matrix (23371 rows, 50 columns) with some NAs in it and it runs for about a minute, then gives an error :
> agnes(iMatrix)
Error in agnes(iMatrix) :
No clustering performed, NA-values in the dissimilarity matrix.
I've also tried getting rid of rows with all NAs in them, and it still gave me the same error. Is this a bug in agnes() ? It doesn't seem to fulfil the claim made by its documentation.
The matrix...
2007 Jun 15
0
Question with nlm
...;- as.matrix(object at data$output[1:totalTime,1]);
x <- as.matrix(object at data$input[1:totalTime,]);
# compute necessary matrices
M <- as.matrix(object at model$M);
P <- diag(sigma*sigma);
A <- AMatrix(totalTime, M, object at data$input[1:totalTime,]);
Q <- IMatrix(totalTime)+A %*% outerM(IMatrix(totalTime-1),P) %*% t(A);
invQ <- solve(Q,IMatrix(dim(Q)[1]));
xM <- matrix(rep(0, dim(M)[2]*totalTime), ncol=dim(M)[2], nrow=totalTime);
for (i in 1:totalTime) {
xM[i,] <- x[i,] %*% powerM(M, -totalTime+i);
}
tmp <- solve((t(xM...
2001 Nov 05
1
Why doesn't outer work?
...;
M<-matrix(0,n,n);
P<-matrix(0,n,n);
B<-matrix(0,n,n);
P<-t(outer(y,y,pxy))
B<-t(outer(y,y,fxy))
M<-P+B
K<-M;
M<-(U-L)*M/n;
P<-(U-L)*P/n;
B<-(U-L)*B/n;
return(list(matrix=M,kernel=K,meshpts=y,Pmatrix=P,Bmatrix=B,Imatrix=I));
}
bigmatrix<-function(n) {
# upper and lower integration limits
L<-minsize; U<-maxsize;
# boundary points b and mesh points y
b<-L+c(0:n)*(U-L)/n;
y<-0.5*(b[1:n]+b[2:(n+1)]);
# loop to construct the matrix
I<-diag(n);
M<-matrix(0,n,n);
P<...
2012 Jul 21
2
car::Anova - Can it be used for ANCOVA with repeated-measures factors.
...ted (i.e., does not interact with the other factors), but when calling Anova on the model, I don't know how I can specify the between-within design (i.e., which parts of the model should interact with the repeated measures factors).
As far as I understand it, neither the idesign, icontrasts or imatrix arguments, nor the linearHypothesis function can specify the within-between design (as far as I get it they all specify the within or intra-subject design, see John Fox's slides from User 2011: http://web.warwick.ac.uk/statsdept/useR-2011/TalkSlides/Contributed/17Aug_1705_FocusV_4-Multivariate_...
2010 Jul 29
0
[R-pkgs] heplots 0.9-3 and candisc 0.5-18 released to CRAN
...th icontrasts
o Added markH0 argument to heplot
o Added label argument to mark.H0()
o Added arrow3d()
o Fixed buglet in pairs.mlm (type)
Version 0.9-0 (2009-12-25)
o Implemented 2D HE plots in heplot.mlm for repeated measures designs
by adding
idata=, idesign=, icontrasts= arguments.
o imatrix= argument for doubly-multivariate designs implemented, but
depends on car 2.0+
o Updated example(OBK)
o Updated example(WeightLoss)
o Implemented 3D HE plots in heplot3d.mlm for repeated measures designs
(OK)
o Implemented pairs() for repeated measures designs (OK)
== candisc
candisc includ...
2010 Aug 23
3
extracting p-values from Anova objects (from the car library)
Dear all,
is there anyone who can help me extracting p-values from an Anova object
from the car library? I can't seem to locate the p-values using
str(result) or str(summary(result)) in the example below
> A <- factor( rep(1:2,each=3) )
> B <- factor( rep(1:3,times=2) )
> idata <- data.frame(A,B)
> fit <- lm( cbind(a1_b1,a1_b2,a1_b3,a2_b1,a2_b2,a2_b3) ? sex,
2010 Aug 24
0
mlm for within subject design
...s "1","2","3": 1 2 3 1 2 3
.. ..- attr(*, "contrasts")= chr "contr.sum"
$ idesign :Class 'formula' length 2 ~A * B
.. ..- attr(*, ".Environment")=
$ icontrasts: chr [1:2] "contr.sum" "contr.poly"
$ imatrix : NULL
- attr(*, "class")= chr "Anova.mlm"
> str(summary(result))
Type III Repeated Measures MANOVA Tests:
------------------------------------------
Term: (Intercept)
Response transformation matrix:
(Intercept)
a1_b1 1
a1_b2 1
a1_b3...