Yaacov,
I cannot really help you, I don't even understand what you're trying to do!
What's the link between your function and all the information you've
just sent (by the way, you forgot to send it to the list)?
By the way, attachments are not really advised, you'd better use dput()
I've never used these packages and functions, that could explain why I
don't understand.
Maybe someone else will be able to help you
Ivan
-------- Message original --------
Sujet: RE: [R] dims error
Date : Tue, 17 Aug 2010 09:11:18 -0400
De : Yaacov Petscher <ypetsc...@fcrr.org>
Pour : ivan.calan...@uni-hamburg.de <ivan.calan...@uni-hamburg.de>
Thanks for your response Ivan. Below is the call and attached are the
referenced tables. My sense is that it's related to the table, but I'm not
quite sure why.
library(statmod)
library(rgenoud)
aE<- read.table("Z:\\Measurement Grant\\TCC linking\\f1p1.dat",
header=T,sep="\t")$mean[1:20]
bE<- read.table("Z:\\Measurement Grant\\TCC linking\\f1p1.dat",
header=T,sep="\t")$mean[31:50]
cE<- rep(0,25)
aB<- read.table("Z:\\Measurement Grant\\TCC linking\\f1p28.dat",
header=T,sep="\t")$mean[1:20]
bB<- read.table("Z:\\Measurement Grant\\TCC linking\\f1p28.dat",
header=T,sep="\t")$mean[31:50]
cB<- rep(0,25)
theta<-seq(-4, 4, length=20)
nxi<- 40
std<- sqrt(c(0.306))
item<- c(9,8)
ax<-aB
cx<-cB
bx<-bB
ay<-aE
cy<-cE
by<-bE
testlet.linking.coefficients(ax, cx, bx, theta, nxi, std, item, ay,
cy, by,"Z:\\Measurement Grant\\TCC linking\\f1P28p1OUT.txt")
Yaacov Petscher, Ph.D.
Director of Research, Florida Center for Reading Research
Associate in Research, Florida State University
2010 Levy Ave
Suite 100
Tallahassee, Florida, 32310
850-645-8963 (phone Tuesday/Thursday)
850-645-9477 (phone Monday/Wednesday/Friday)
850-644-9085 (f)
-----Original Message-----
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Ivan Calandra
Sent: Tuesday, August 17, 2010 9:08 AM
To: r-help@r-project.org
Subject: Re: [R] dims error
Hi,
Maybe knowing what you provided to the function call as a, b, c, and xi would
help!
Since you're at it, think about using the dput() function to provide
copy/pastable objects ;)
HTH,
Ivan
Le 8/17/2010 15:00, Yaacov Petscher a écrit :
Greetings,
I am a very novice user with R, and in the course of running a linking
procedure :
P.old<- function (a, c, b, xi){
#a contains a parameters
#c contains c parameters
#b contains b parameters
#xi is a one column vector containing quadrature
points for xi(=theta-gamma)for one item
tmp<- a*(xi-b)
p<- c+(1-c)*pnorm(tmp)
return(p)
}
I receive an error which states:
Error in P.old(a[i], c[i], b[i], node) :
dims [product 40] do not match the length of object [0]
I'm not exactly sure to what I should attend in changing the code. Any help
would be greatly appreciated!
Yaacov
Yaacov Petscher, Ph.D.
Director of Research, Florida Center for Reading Research Associate in
Research, Florida State University 2010 Levy Ave Suite 100
Tallahassee, Florida, 32310
850-645-8963 (phone Tuesday/Thursday)
850-645-9477 (phone Monday/Wednesday/Friday)
850-644-9085 (f)
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
--
Ivan CALANDRA
PhD Student
University of Hamburg
Biozentrum Grindel und Zoologisches Museum
Abt. Säugetiere
Martin-Luther-King-Platz 3
D-20146 Hamburg, GERMANY
+49(0)40 42838 6231
ivan.calan...@uni-hamburg.de
**********
http://www.for771.uni-bonn.de
http://webapp5.rrz.uni-hamburg.de/mammals/eng/mitarbeiter.php
[[alternative HTML version deleted]]
3.016 -5.804
2.569 -5.053
2.519 -4.333
1.331 -2.452
2.927 -5.183
1.149 -2.37
0.629 0.341
0.911 -0.951
1.437 -2.428
3.16 -3.811
1.024 -1.137
1.853 -2.176
3.365 -4.084
3.157 -3.162
1.491 -1.276
1.592 -1.102
2.724 -3.626
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.