Hi all,
I ran into a strange error: I am trying to create a package skeleton for a
new source package from within a function. Objects that are created in this
function are to be included in my package, but for some reason, I get an
error message saying that these objects cannot be found.
Here is
Christian Ritz pisze:
Hi Jarek,
an alternative approach is to provide more precise starting values!
It pays off to realise that it's possible to find quite good guesses
for some of the parameters in your model function:
t ~ tr+(ts-tr)/((1+(a*h)^n)^(1-(1/n)))
The parameters
Indeed, this works.
Thanks!
On Jan 28, 2008 9:30 AM, Prof Brian Ripley [EMAIL PROTECTED] wrote:
You need to set the 'environment' argument (the help file is incomplete).
e.g.
env - sys.frames()[[sys.nframe()]]
package.skeleton(name = pkgName, list=c(f,e, myenv), env=env)
On Mon, 28 Jan
If you type
example(barplot)
you will find an example.
Ciao,
domenico
CHENYS wrote:
Hi, I'm looking for a tool which can fill bar chart with dash, skewed line,
or grids, rather than pure color. Any one have the idea how to do that in R?
Or maybe in Matlab will also be helpful.
Thanks
Dear Rainer,
I think you'll have to recode the labels of the levels to get the output
that you want.
levels(ssq$gi) - c(gi = 1, gi = 2, gi = 3)
HTH,
Thierry
ir. Thierry Onkelinx
Instituut voor natuur- en
Groot, Philip de wrote:
Hello all,
I am not sure whether it actually is a bug, but it is not the behaviour I
would expect. Please consider this:
Sibships
[1] Patient_2400 Patient_2400 Patient_345 Patient_345 Patient_8901
[6] Patient_8901 Patient_4008 Patient_4008
hey all
greetings
hey all am an engineering student...and am trying to learn R
i am trying to automate reading a specific type of file...and perform
certain functions...but i want to omit lines in the end of the file..
there is an option for skiping the lines before begining...but how can i ask
Hello all,
I am not sure whether it actually is a bug, but it is not the behaviour I would
expect. Please consider this:
Sibships
[1] Patient_2400 Patient_2400 Patient_345 Patient_345 Patient_8901
[6] Patient_8901 Patient_4008 Patient_4008 Patient_7991 Patient_7991
[11] Patient_8353
Hello;
Please coud you advise me of a simple way to select points (x,y
coordinates) that fall within a polygon.
I've got a set of polygons, each one defined by an arbitrary number of
points, and several points inside each polygon.
I know this is simple with a GIS, but I'd rather do it inside R.
Dear useRs,
I've release a new version of the dynlm package to CRAN which adds two
new features:
o instrumental variables regression (two-stage least squares) via
formulas like
dynlm(y ~ x1 + x2 | z1 + z2 + z3, data = mydata)
where z1, z2, z3 are the instruments which can again
Perhaps:
data - read.table(textConnection(rev(rev(readLines('data.txt'))[-(1:2)])))
will skip the last two lines '-(1:2)'
On 28/01/2008, mrafi [EMAIL PROTECTED] wrote:
hey all
greetings
hey all am an engineering student...and am trying to learn R
i am trying to automate reading a specific
jgarcia at ija.csic.es writes:
Hello;
Please coud you advise me of a simple way to select points (x,y
coordinates) that fall within a polygon.
I've got a set of polygons, each one defined by an arbitrary number of
points, and several points inside each polygon.
I know this is simple with
pieterprovoost at gmail.com writes:
No, I get the same error message there...
Please do not needlesly delete the thread content. Your original question was:
I'm having problems reading a shapefile with read.shape (maptools). I'm
absolutely sure my file is there, but I get no such file. The
Professor Ripley,
do you have an idea why it works for me to save different types of
objects, but not for the environment object I construct in the
function:
### Example of a function that works without error:
fun2 - function(myname){
f - function(x,y) x+y
g - function(x,y) x-y
d -
Roger Bivand Roger.Bivand at nhh.no writes:
pieterprovoost at gmail.com writes:
No, I get the same error message there...
Your original question was:
I'm having problems reading a shapefile with read.shape (maptools). I'm
absolutely sure my file is there, but I get no such file.
hello Duncan, Gabor
Many thanks for your help!
I think that the line:
if( chartr(/,\\,getwd() )==R.home() ) setwd(C:\\)
can solve my problem in this moment!! :-)
Thanks again
Cleber
On 26/01/2008 9:03 PM, Cleber Nogueira Borges wrote:
hi Gabor and Duncan,
I make a test and I find
gallon li [EMAIL PROTECTED] wrote in
news:[EMAIL PROTECTED]:
Does anybody know if there is such a function to estimate the
distribution for interval censored data?
survfit doesn't work for this type of data as I tried various
references.
My reading of R-help posts by Therneau and Lumley
Hi,
I want to plot a graph and here is my code:
ec-rep(0,length(e))
fc-rep(0,length(f))
plot(e,ec,type=p,col=1,pch=19)
points(f,fc,col=2,pch=20)
legend(1.0e+08,1.0,c(dog, human),text.col=green4,pch=c(19,20),col=c(1,2))
my major problem here is the x-axis is too large in
And subsetting a factor retains the original factor levels. To drop
unused levels, just use factor(f[index]) or f[index, drop=TRUE]. The
opposite behaviour can be even more annoying/dangerous because it leads
to empty cells dropping out of tables and bars disappearing from barplots.
Of course
ggopt(strip.text = function(variable, value) paste(variable, value, sep==
))
That's exactly what I was looking for - thanks.
One thing that I should mention is that this is likely to change at
some point in the future. Eventually it will become:
+ facet_grid(strip.text = function(variable,
Dear list,
I have following data, where I want to check if any value in each
column is out of range. For example, column f1 can only take values
1-5, so if any values less than 1 or 5 will be
defined as missing value (i.e. NA), column f4 can only take values
of 1-3 and any values that
On Mon, 28 Jan 2008, Barry Rowlingson wrote:
Henrique Dallazuanna wrote:
Perhaps:
data - read.table(textConnection(rev(rev(readLines('data.txt'))[-(1:2)])))
Euurgh! Am I the only one whose sense of aesthetics is enraged by
this? To get rid of the last two items you reverse the vector,
Hello r-help,
I have a lengthy vector of data (with values anywhere from 1-200), and
another index vector of 'groups' representing values 0-2, 3-5, 6-8, ...
of length 67. The index vector has the structure (1, 4, 7, ... , 196,
199), where each value is the midpoint of each respective group.
I'm
Dear list,
I wonder if somebody has succeeded in installing R on an eeePC (Xandros
desktop). Searching via Rseek (term eeePC) and in eeePC forums (term
Cran) left me without proper hits.
Best wishes,
Walter.
__
R-help@r-project.org mailing
I'm a new user and am having trouble with loops.
In the following, I'm trying to add the results of test and the loops are
not working.
I've simplified the loop. What am I doing wrong?
Thanks!
test-numeric(20)
tot-numeric(20)
for(i in 1:20){test[i]-1}
for (i in 1:20){tot[i]-(test[i]+tot[i])}
Try:
transform(data, f1 = factor(f1, levels = 1:5), f4 = factor(f4, 1:3))
On Jan 28, 2008 8:12 AM, Tom Cohen [EMAIL PROTECTED] wrote:
Dear list,
I have following data, where I want to check if any value in each
column is out of range. For example, column f1 can only take values
1-5, so
If you don't mind reading it in twice its just:
DF - read.table(xy.dat, header = TRUE, nrow =
length(readLines(xy.dat)) - 3)
tail(DF)
# or
DF - read.table(xy.dat, header = TRUE, nrow =
length(count.fields(xy.dat)) - 3)
tail(DF)
# or
DF - read.table(xy.dat, header = TRUE, nrow =
Hi
[EMAIL PROTECTED] napsal dne 28.01.2008 13:06:16:
Hi,
I want to plot a graph and here is my code:
ec-rep(0,length(e))
fc-rep(0,length(f))
Strange. If I understand correctly, you repeat zero length(e) or length(f)
times and then you plot those zeroes at x axis = e values and
On Mon, 28 Jan 2008, mrafi wrote:
but then the number of levels would reamain the same...!!
Please explain: the levels of factors are taken from the data which is
actually read.
Prof Brian Ripley wrote:
On Mon, 28 Jan 2008, Barry Rowlingson wrote:
Henrique Dallazuanna wrote:
Perhaps:
Thank you, and also thanks to John Kane.
I did some tests last night, if you plot something like
height - t(t(c(1,-1,1)))
bardensity - t(t(c(10,10,0)))
barangle - t(t(c(45,135,0)))
barplot(height, density = bardensity, angle = barangle)
You can get grids. Otherwise, only slanted lines or
density ?barplot
aa - c(4,5,6)
barplot(aa, density=2, col='red', border=blue)
--- CHENYS [EMAIL PROTECTED] wrote:
Hi, I'm looking for a tool which can fill bar chart
with dash, skewed line,
or grids, rather than pure color. Any one have the
idea how to do that in R?
Or maybe in Matlab
but then the number of levels would reamain the same...!!
Prof Brian Ripley wrote:
On Mon, 28 Jan 2008, Barry Rowlingson wrote:
Henrique Dallazuanna wrote:
Perhaps:
data -
read.table(textConnection(rev(rev(readLines('data.txt'))[-(1:2)])))
Euurgh! Am I the only one whose sense of
The Icens package implements a number of different estimation paradigms.
Briefly, for interval censored data, the KM does not apply, but
Turnbull's method does. Relevant citations for Icens are in it.
For a Cox model approach, there was at one time a package, not at CRAN,
by Commenges (and a
Ok, I agree, but I assume that her don't know the number of rows of file.
In your example:
rev2 - function(x,n=100,m=5){
for(i in 1:n){
y=x[1:(length(x)-m)]
}
return(y)
}
is needed that open other textConnection - if use the example posted by me.
Is there other option?
On
Tobin, Jared wrote:
Hello r-help,
I have a lengthy vector of data (with values anywhere from 1-200), and
another index vector of 'groups' representing values 0-2, 3-5, 6-8, ...
of length 67. The index vector has the structure (1, 4, 7, ... , 196,
199), where each value is the midpoint of
Doing 'RSiteSearch(eee)' yields some hits. I knew that the ASUS eeePC
had come up on r-help.
-kevin
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
On Behalf Of Dr. Walter H. Schreiber
Sent: Monday, January 28, 2008 9:32 AM
To: r-help@r-project.org
Subject: [R] R on
On 1/28/2008 8:48 AM, cvandy wrote:
I'm a new user and am having trouble with loops.
In the following, I'm trying to add the results of test and the loops are
not working.
I've simplified the loop. What am I doing wrong?
Thanks!
test-numeric(20)
tot-numeric(20)
for(i in 1:20){test[i]-1}
Hi, all
I've written some R script to calculate the linear regression of a matrix.
Here below is my script:
x-matrix(scan(h:/data/xxx.dat,0),nrow=46,ncol=561,byrow=TRUE)
year - NULL
year - cbind(year,as.matrix(x[,1]))
lm.sol-lm(x~year)
xtrend-coef(lm.sol)[2,]# get the matrix of
hi Duncan,
I can't access file greater than 3 mb in my work! :-( :-(
((( unfortunately, is not a joke! I need a new job :-D )))
When I return home, I will do download the patch and
report the results to you! :-(
Thanks for your attention!
Cleber
It is now fixed in the version
RobG == Robert Gentleman [EMAIL PROTECTED]
on Mon, 28 Jan 2008 08:18:06 -0800 writes:
RobG The Icens package implements a number of different estimation
paradigms.
RobG Briefly, for interval censored data, the KM does not apply, but
RobG Turnbull's method does. Relevant
hits=-2.6 tests=BAYES_00
X-USF-Spam-Flag: NO
Hi Michelle,
You don't show your read.csv or read.table call, nor the output of
str(obj) where obj is the name of the object you read the data into.
I notice that you have explicit 0 and NA. Is there a chance that you
have entered NA into the cells
Hello,
I am trying to create multiple matrices (to run a PVA) but can't import all
of them from a .csv without the numbers treated as labels and not factors.
I can enter the matrix slowly:
Site05_96 - matrix(c(0.07,0,0.03,0.00,NA,0.00,
0.09,0.16667,0.31,0.42,NA,0.00,
See 'fig' argument in ?par
example:
x - rnorm(100)
plot(x)
par(fig=c(2/3,1,2/3,1), new=T)
hist(x, main=)
On 28/01/2008, lamack lamack [EMAIL PROTECTED] wrote:
From: [EMAIL PROTECTED]: [EMAIL PROTECTED]: sub-plot in a plotDate: Mon, 28
Jan 2008 15:43:40 +
Dear all, how can I do
Hi
[EMAIL PROTECTED] napsal dne 28.01.2008 17:16:37:
Hi, all
I've written some R script to calculate the linear regression of a
matrix.
Here below is my script:
x-matrix(scan(h:/data/xxx.dat,0),nrow=46,ncol=561,byrow=TRUE)
year - NULL
year - cbind(year,as.matrix(x[,1]))
or even:
head(read.table(xy.dat, header = TRUE), -2)
On Jan 28, 2008 10:52 AM, Gabor Grothendieck [EMAIL PROTECTED] wrote:
If you don't mind reading it in twice its just:
DF - read.table(xy.dat, header = TRUE, nrow =
length(readLines(xy.dat)) - 3)
tail(DF)
# or
DF - read.table(xy.dat,
Dear Thierry
On 28/01/2008, ONKELINX, Thierry [EMAIL PROTECTED] wrote:
Dear Rainer,
I think you'll have to recode the labels of the levels to get the output
that you want.
I haven't thought about using levels - that was probably the reason why it
threw the sorting out. Thanks
Rainer
Come join me on MaLaYsiaN BesT!
Mari kita berkongsi cerita mengenai apa sahaja mengenai malaysia. cth: gambar
lucu, gambar aneh, gambar hantu, artis malaysia, info bola sepak malaysia,
sejarah bola sepak, berita sukan bola sepak terkini dan arena bola sepak .
Click here to join:
See the subplot function in the TeachingDemos package.
Hope this helps,
From: [EMAIL PROTECTED] on behalf of lamack lamack
Sent: Mon 1/28/2008 8:51 AM
To: r-help@r-project.org
Subject: [R] sub-plot in a plot
From: [EMAIL PROTECTED]: [EMAIL PROTECTED]:
Thanks for all who commented on this question.
It turns out that there is a nice set of functions in library BiodiversityR
that translate a stack or list input data file to a community matrix:
import.from.Excel(file = file.choose(), sheet = community, sitenames =
sites,
column = species,
This is not a bug; it is deliberately designed this way.
There are circumstances when you want to drop levels on subsetting and
other circumstances where you don't, so the default behaviour can't make
everyone happy. However, there is an option to get the behaviour you want
Thanks very much for the enlightenment. Very interesting indeed, and
I am glad
to find Nightingale exonerated of her purported crime.
cheers,
Rolf
On 29/01/2008, at 8:25 AM, Greg Snow wrote:
I had heard the same thing about Florence Nightingale, but it seems
Thank you!
cheers,
Rolf
On 29/01/2008, at 8:38 AM, roger koenker wrote:
Howard Wainer (Graphical Discovery, PUP, 2005, p 20) gives
this dubious honor to Playfair (1759- 1823). Nightingale (1820-
1910) was far too enlightened for this sort of thing, see for example
Christophe Genolini wrote:
Hi the list.
I do not understand the philosophy behind numeric and integer.
- 1 is numeric (which I find surprising)
- 2 is numeric.
- 1:2 is integer.
Why is that ?
I hope I can answer your question at least partly:
Numeric means double, i.e. internally
If anyone is interested in seeing Nightingale's Coxcomb a.k.a. Nightingale's
Rose, it can be seen at
http://www.economist.com/world/europe/displaystory.cfm?story_id=10278643.
Best regards,
Peter.
-Oprindelig meddelelse-
Fra: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] På vegne af Rolf
Thomas and David,
Thanks for your answers. The svyvar() function does the trick. It does not
provide a vcov() method for the variance-covariance matrix of these
estimates themselves. But I guess I can bootstrap them using
library(boot)...
Thanks again, daniel
On Jan 25, 2008 6:13 PM, Thomas
Hello
I'm trying to implement a grid search for a threshold autoregressive
model, it is a model in which the regression coefficients are different
according to the regimes (under the lower threshold, between lower and
upper, over the upper threshold).
Estimation of the threshold is made with
Dr. Walter H. Schreiber wrote:
Dear list,
I wonder if somebody has succeeded in installing R on an eeePC
(Xandros desktop). Searching via Rseek (term eeePC) and in eeePC
forums (term Cran) left me without proper hits.
Try just looking for 'eee' at
http://tolstoy.newcastle.edu.au/~rking/R/
Filling bars with lines/grids/points is legacy back to the days when the only
way to get high quality plots was to use a pen plotter (on the screen you would
see bars made of '*' or similar). The pen plotter would use a mechanical arm
to draw the lines using a pen/marker. it was easy to have
On 28/01/2008, at 11:14 PM, [EMAIL PROTECTED] wrote:
Hello;
Please coud you advise me of a simple way to select points (x,y
coordinates) that fall within a polygon.
I've got a set of polygons, each one defined by an arbitrary number of
points, and several points inside each polygon.
I know
leo_aries wrote:
Hi, all
I've written some R script to calculate the linear regression of a
matrix.
Here below is my script:
x-matrix(scan(h:/data/xxx.dat,0),nrow=46,ncol=561,byrow=TRUE)
year - NULL
year - cbind(year,as.matrix(x[,1]))
lm.sol-lm(x~year)
I was asked for the following information and hope it might help those who
could answer my question...
To import the table I used:
AsMi05test=read.csv(C:/AsMi_Site05_1998.csv)
str(AsMi05test)
`data.frame': 12 obs. of 8 variables:
$ X : Factor w/ 6 levels
On 28/01/2008, at 12:07 AM, Peter Dalgaard wrote:
Jean lobry wrote:
snip
about an hour North of Paris. Her father inquired -
coincidentally during the cheese course - what work I was
doing in Paris; I replied that I was researching the
activities of a Scot, William Playfair,
On 28-Jan-08 21:23:12, Roland Rau wrote:
Christophe Genolini wrote:
Hi the list.
I do not understand the philosophy behind numeric and integer.
- 1 is numeric (which I find surprising)
- 2 is numeric.
- 1:2 is integer.
Why is that ?
I hope I can answer your question at least partly:
I had heard the same thing about Florence Nightingale, but it seems that this
is a confusion of different graphs.
Nightingale developed a graph based on a circle, but all the angles were equal
and the different values were encoded by using different radii of the slices
(and she did the
Howard Wainer (Graphical Discovery, PUP, 2005, p 20) gives
this dubious honor to Playfair (1759- 1823). Nightingale (1820-
1910) was far too enlightened for this sort of thing, see for example
her letter to Galton about endowing an Oxford professorship
in social statistics (reprinted in Karl
Roger,
I tried your suggestions and ran into same problems as Pieter did before.
On 28.01.2008 12:35 (UTC+1), Roger Bivand wrote:
pieterprovoost at gmail.com writes:
No, I get the same error message there...
Please do not needlesly delete the thread content. Your original question was:
http://cran.fhcrc.org/bin/macosx/old/powerpc/R-2.2.0.dmg
But that is sooo old.
Maybe he could compile from the source (he'd need to install XCode and
a fortran compiler).
I assume that once these requirements are installed, the compilation
would be as simple as:
./configure
make
Hello,
I am Mr. Ellis David by name and I?m the CEO/OWNER of Davidellis
Textiles company Inc. in West Africa. We supply textiles to over
hundred of customers spread all over the world.We need someone
that will be representing my company as a
Manager/Representative, who will attend to customers in
RSiteSearch(R asus)
G.
On Mon, Jan 28, 2008 at 03:31:45PM +0100, Dr. Walter H. Schreiber wrote:
Dear list,
I wonder if somebody has succeeded in installing R on an eeePC (Xandros
desktop). Searching via Rseek (term eeePC) and in eeePC forums (term
Cran) left me without proper hits.
Rolf Turner wrote:
On 26/01/2008, at 10:54 AM, Carson Farmer wrote:
Dear List,
I am attempting to perform a harmonic analysis on a time series of snow
depth, in which the annual curve is essentially asymmetric (i.e. snow
accumulates slowly over time, and the subsequent melt occurs
Hi all
I am trying to generate a normal unbalanced data to estimate the coefficients
of LM, LMM, GLM, and GLMM and their standard errors. Also, I am trying to
estimate the variance components and their standard errors. Further, I am
trying to use the likelihood ratio test to test H0: sigma^2_b
On 28-Jan-08 22:40:02, Peter Dalgaard wrote:
[...]
AFAIR, space is/was more of an issue. If you do something like
for i in 1:1e7
some.silly.simulation()
then you have 40 MB sitting there doing nothing, and 80 MB if
it had been floating point.
Hmmm ... there's something to be said
On Mon, 28 Jan 2008, Rainer Hurling wrote:
Roger,
I tried your suggestions and ran into same problems as Pieter did before.
Rainer:
I don't think that we know what platform Pieter was using, apart from the
working directory that looked like Windows, and where the problem could
not be
(Ted Harding) wrote:
Further to the above: The help
?:
says:
Value:
For numeric arguments [as opposed to factors],
a numeric vector. This will be of type 'integer'
if 'from' and 'to' are both integers and
representable in the integer type, otherwise of
type 'numeric'.
By
I am trying to install the quantreg Package on MacOSX. When I choose CRAN
mirror, it gives me the error:
Error in read.dcf(file = tmpf) : Line longer than buffer size
In addition: There were 44 warnings (use warnings() to see them)
Is there any known current problem with this package?
Sincerely
Further to the above: The help
?:
says:
Value:
For numeric arguments [as opposed to factors],
a numeric vector. This will be of type 'integer'
if 'from' and 'to' are both integers and
representable in the integer type, otherwise of
type 'numeric'
???
This is very
I know nothing of Macintosh, so please be patient.
My student has a Macintosh with OSX 10.3.9. The R for Mac says it is
for 10.4.4 or higher.
Aside from saying get a new Mac, what can be said to my student?
Can you point me at the newest version of R that did work on 10.3 ?
pj
--
Paul E.
Hello,
I'm trying to compile Shogun 0.4.4 from source on Redhat Enterprise
Linux 4 (update 5). I've compiled and installed R-2.6.1 already and
that seems to work fine. I had some initial problems finding R
development files, which were fixed by setting CFLAGS='-I/tmp/R-2.6.1/
include' and
Dear R community,
I'm curious to know how people go about estimating standard errors for
parameter estimates after model selection by ridge regression and the
lasso. Do you have any practical or theoretical advice?
Warmly,
Andrew
--
Andrew Robinson
Department of Mathematics and Statistics
hits=-2.6 tests=BAYES_00
X-USF-Spam-Flag: NO
On Mon, 2008-01-28 at 12:17 -0700, Michelle DePrenger-Levin wrote:
I was asked for the following information and hope it might help those who
could answer my question...
That looks fine to me Michelle.
You will have problems with as.matrix on this
On Jan 28, 2008 1:25 PM, Greg Snow [EMAIL PROTECTED] wrote:
I had heard the same thing about Florence Nightingale, but it seems that this
is a confusion of different graphs.
Nightingale developed a graph based on a circle, but all the angles were
equal and the different values were encoded
I have a data vector as following:
z
[1] 183.1370 201.9610 113.7250 140.7840 156.2750 42.1569 42.1569 42.1569
[9] 240.1960 308.4310 42.1569 42.1569 42.1569 42.1569 42.1569 42.1569
[17] 42.1569 42.1569 42.1569 42.1569 279.8040 42.1569 42.1569
when I sort, it gave me the right
that seems right
order() gives you the indexes idx such that x[idx] == sort(x)
set.seed(123)
x - rnorm(10)
idx - order(x)
identical(x[idx], sort(x))
[1] TRUE
best
b
On Jan 28, 2008, at 8:19 PM, Waverley wrote:
I have a data vector as following:
z
[1] 183.1370 201.9610 113.7250
On Mon, 28 Jan 2008, Waverley wrote:
I have a data vector as following:
z
[1] 183.1370 201.9610 113.7250 140.7840 156.2750 42.1569 42.1569 42.1569
[9] 240.1960 308.4310 42.1569 42.1569 42.1569 42.1569 42.1569 42.1569
[17] 42.1569 42.1569 42.1569 42.1569 279.8040 42.1569
well if you want to find the spectral density aka what frequencies
explain most of the variance then I would suggest the spectral
density. This can be implemented with spec.pgram(). This is
conducted with the fast fourier transform algorithm.
a-ts(data, frequency = 1) #make the time series
On Jan 28, 2008, at 11:30 AM, Liaw, Andy wrote:
The original paper is:
Cook, J. R. Stefanski L. A. (1994) Simulation--extrapolation
estimation in parametric measurement error models. Journal of the
American Statistical Association, 89, 1314-1328.
From: Michael Kubovy
Dear
Hi,
For example, given two 2x2 matrices m1 and m2. I would like to add/subtract
element by element
m1
[,1] [,2]
[1,] NA NA
[2,]12
m2
[,1] [,2]
[1,]1 NA
[2,] NA2
m1 + m2
[,1] [,2]
[1,] NA NA
[2,] NA4
How can I ignore the NA, and get
Hi,
does anyone know how to simulate two seasonal data series that are
cointegrated?
Thanks!
--
Tom
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
Perhaps you could explain the motivation behind this. At any rate here
are three different solutions:
ifelse(is.na(m1), ifelse(is.na(m2), NA, m2), ifelse(is.na(m2), m1, m1 + m2))
apply(array(c(m1, m2), c(2,2,2)), 1:2, function(x) sum(c(na.omit(x), NA)[1]))
na.m1 - is.na(m1)
na.m2 - is.na(m2)
Your problem is that your newdata data frame only specifies what the
value for A is in the predictions. The values for W1 and W2 are
unspecified. To predict from fitted models, you need to specify the
values you wish to use for *all* the predictor variables, not just the
one (I presume) is
According to the R documentation for NLMINB, the returned value of convergence
is 0 for successful convergence. When I got another code (1), I looked up the
PDF that linked from the documentation
(http://netlib.bell-labs.com/cm/cs/cstr/153.pdf), which said that a return code
under 3 was
Dear all;
What can be wrong with this simple example?
library(boot)
d1-c(rnorm(10,mean=10))
fm-function(d,i) mean(d[i])
bd1-boot(d1,fm,1)
Error: evaluation nested too deeply: infinite recursion /
options(expressions=)?
Thanks for any idea
__
I wonder what your sessionInfo() is...
I copied and pasted your example and could not reproduce the problem
(ie, worked just fine).
b
On Jan 29, 2008, at 12:27 AM, Pedro Mardones wrote:
Dear all;
What can be wrong with this simple example?
library(boot)
d1-c(rnorm(10,mean=10))
Hello,
I used 'arrows' function to plot the error bar but it dosen't work always.
For example:
Similar commands on two different data sets gives first result OK and second
result NOT (attached plot1 is OK and plot2 is NOT). Plots as well as data
file is attached here. Below is the script which
94 matches
Mail list logo
