There are several mistakes:
Why are you using only p[1], p[2]?
Afterwards you try (for i in 1:10) p[i] so you should go to p[10]?
Then you should make a pnew wich holds 10 values, so pnew[i]
This should get you back on track
Bart
Kathy Maher wrote:
Hi,
I am trying to use the Fisher
Hi,
I am trying to use the Fisher scoring method with a geometric distribution,
with p = .07, 100 observations from the geom distrib, and 10 iterations.
I cannot quite get the code to work.
Can anyone see the mistake?
n - 100
p - 0.07
x - rgeom(n, p)
s - sum(x)
f - function(x, p)
Hi there,
Am try to find tree diameter increment in a hierarchical data with tree
diameter measurements, date of measurements, in plots that are located in
different region. How do I do it? am very new to R.
Thanks for you help.
Maina Kariuki
Level 3, 8 Nicholson Street
East Melbourne VIC
PROTECTED]
Sent: Friday, 14 March 2008 9:03 AM
To: r-help@r-project.org
Subject: [R] (no subject)
Hi there,
Am try to find tree diameter increment in a hierarchical data with tree
diameter measurements, date of measurements, in plots that are located
in
different region. How do I do it? am very new
# PROBLEM WITH ABLINE
# I have a question about a seemingly imprecise abline figure which I have
come accross now a few times.
# After reading variables a and b...
imprecise are our eyes :)
check where 4336.543 actually is and note how the range of the data
confused you.
plot(a~b)
abline(model)
abline(h=4336.543, v=10)
b
On Mar 11, 2008, at 11:27 AM, Aad Termorshuizen wrote:
# PROBLEM WITH ABLINE
# I have a question about a seemingly imprecise
On Tuesday 11 March 2008 04:27:12 pm Aad Termorshuizen wrote:
AT # So my QUESTION is: is abline imprecise or am I making errors?
Either you are not precise or I did not got your problem.
plot(a~b)
abline(model)
points(10,4336,col=red)
as you see the predicted point (red) is exactly at the
Hi all,
I am trying to read a text file in R! I have a warning message:
In read.table(file = data1.txt, header = T) : uncompleted last line of
readTableHeader in 'data1.txt'
Could you please tell me why? and how can I deal with missing value when I read
this file?
all the best
Alfons
message?
Thank you in advance!
- Original Message
From: Duncan Murdoch [EMAIL PROTECTED]
To: Alfons Sutter [EMAIL PROTECTED]
Cc: [EMAIL PROTECTED]
Sent: Monday, March 10, 2008 12:18:11 PM
Subject: Re: [R] (no subject)
On 10/03/2008 6:50 AM, Alfons Sutter wrote:
Hi all,
I am trying
On 10/03/2008 6:50 AM, Alfons Sutter wrote:
Hi all,
I am trying to read a text file in R! I have a warning message:
In read.table(file = data1.txt, header = T) : uncompleted last line of
readTableHeader in 'data1.txt'
R is warning you that the file may be incomplete. In Unix, text files
and could you tell me place
how can I avoid this warning message?
Thank you in advance!
- Original Message
From: Duncan Murdoch [EMAIL PROTECTED]
To: Alfons Sutter [EMAIL PROTECTED]
Cc: [EMAIL PROTECTED]
Sent: Monday, March 10, 2008 12:18:11 PM
Subject: Re: [R] (no subject
Dear all,
I try to import a SPSS.por dataset with about 6000 cases and 650
variables with R commander
and got error messages:
'Error: /temp/t.por use.value.labels=true,
max.value.label=inf,
'Error: data is not data frame and cannot be attached'
Any comment?
Thanks in
Trying to set up multiple spawns of R
---
Hi. I'm new to the list but I have a small problem with R. I'm trying to
set up multiple spawns of R active at the same time on an AMD with a dual
core in Vista so that I can cut down my permutation processing
[EMAIL PROTECTED] napsal dne 02.02.2008 07:10:37:
bgchen wrote:
r-help
c(2.43, 3.22, 6.9, 3.03, 5.36, 6.9, 2.29, 6.13, 6.11, 4.25, 3.85,
5.09, 7.44, 2.86, 2.82, 3.64, 3.22, 7, 2.65, 4.5, 3.73, 5.09,
5.8, 7.87, 2.87, 2.9, 6.63, 6.8, 2.45, 7.68, 2.56, 2.54, 7.35,
4.61, 2.58, 3.27,
r-help
c(2.43, 3.22, 6.9, 3.03, 5.36, 6.9, 2.29, 6.13, 6.11, 4.25, 3.85,
5.09, 7.44, 2.86, 2.82, 3.64, 3.22, 7, 2.65, 4.5, 3.73, 5.09,
5.8, 7.87, 2.87, 2.9, 6.63, 6.8, 2.45, 7.68, 2.56, 2.54, 7.35,
4.61, 2.58, 3.27, 5.8, 3.13, 3.29, 2.56, 5.79, 2.67, 7.55, 7.13,
8.4, 6.16, 3.34, 6.49, 1.85,
bgchen wrote:
r-help
c(2.43, 3.22, 6.9, 3.03, 5.36, 6.9, 2.29, 6.13, 6.11, 4.25, 3.85,
5.09, 7.44, 2.86, 2.82, 3.64, 3.22, 7, 2.65, 4.5, 3.73, 5.09,
5.8, 7.87, 2.87, 2.9, 6.63, 6.8, 2.45, 7.68, 2.56, 2.54, 7.35,
4.61, 2.58, 3.27, 5.8, 3.13, 3.29, 2.56, 5.79, 2.67, 7.55, 7.13,
8.4,
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
Hi, The following code, from Angelo Canty article on line
Resampling Methods in R: the boot Package, 2002, works fine for Angelo Canty
using R 2.6.0 on Windows XP.
It also works for me using R 1.2.1 and S-PLUS 2000 on Windows XP after
installing the S-PLUS
I've searched around and havent had much luck trying to explain why my
HTML search no longer works. If I use
?cor.test
it is unable to display the webpage. Is this a java problem? I am
running windows xp.
Thanks
--
David Kowalske Jr
Graduate
On 1/9/2008 4:28 PM, [EMAIL PROTECTED] wrote:
I have a couple of simple questions.
(a) Suppose that I have an Excel dataset stored in a file called
dataset.xls located in some folder in my hard drive. Suppose that it
has
two columns called X and Y of equal length. After I click on
Jiaming Zuo [EMAIL PROTECTED] wrote in
news:[EMAIL PROTECTED]:
I am working for the United Nations to construct a complete life
table
from an abridged table.
I want to use the code of Hydman Filter by Rob J Hydman but an
error
sentence always appears and it simply
Hi,
I don't think you have given enough information for
anyone to offer advice.
A small working (or almost working ) exmple that
reproduces the problem would help.
--- Jiaming Zuo [EMAIL PROTECTED] wrote:
Dear R Users,
I am working for the United Nations to construct
a complete
Dear R Users,
I am working for the United Nations to construct a complete life table
from an abridged table.
I want to use the code of Hydman Filter by Rob J Hydman but an error
sentence always appears and it simply doesn't run--
source(C:/R/Jamie/HymanFilter.R)
Error in
Dear R Users,
I am working for the United Nations to construct a complete life table
from an abridged table.
I want to use the code of Hydman Filter by Rob J Hydman but an error
sentence always appears and it simply doesn't run--
source(C:/R/Jamie/HymanFilter.R)
Error in
Hi,
I have a data frame that looks something like this:
id dayk
656 -1 566
663 -1 680
673 -1 773
675 -1 761
704 -1 685
714 -1 636
717 -1 421
719 -1 645
727 -1 761
731 -1 663
11 735 -1 603
12 738 -1 865
13 742 -1 594
14 744
hello,
i am a bit of a statistical neophyte and currently trying to make some sense of
confidence intervals for correlation coefficients. i am using the cor.test()
function. the documentation is quite terse and i am having trouble tieing up
the output from this function with stuff that i have
Hello,
This is my first query to the message board. I have a question on how to do a
long run multiplier in R. Our teacher has given us insight in completing this
in SAS, but not in R. Any advice would be appreciated.
Thanks,
Joshua Johns
__
You have this sequence of code in your script:
for (a in seq(1,(ncol(X)-3), by=2)){
for (b in seq(3,(ncol(X)-1), by=2)){
for (i in a){
j - a+1
for (k in b){
l - b+1
The 'for (i in a){' appears to meaningless since 'a' is just a single
value. What did you think this is supposed to do?
# Hello
# I have a question regarding pairwise calculations of a matrix using a
for-loop.
# Below I have a matrix X with 8 columns. These are Genotypic data so
Column1 Column2 is
# a unit, Column3 Column4 is a unit, Column5 Column6 is a unit, and
Coulmn7 8 is a unit.
# I have a loop designed
Hi
I have a collection of about 16 time series with occasional missing
data. A few of these time-series start later than the rest. There is a
relatively high correlation between them (they are hourly temps at
various locations around the UK). The longest series contains about
4 points
I
Hello,
I problem is in the format of the date, my time series is like this:
2006070100 1244 6162
2006070101 1221 6060
2006070102 1214 6060
2006070103 1194 5959
2006070104 1182 5858
2006070105 1178
Joao,
I can't reproduce your results. Try:
str(A)
When I read in the data frame A, I see
str(A)
'data.frame':10 obs. of 4 variables:
$ date: int 2006070100 2006070101 2006070102 2006070103 2006070104
2006070105 2006070106 2006070107 2006070108 2006070109
$ my1 : int 1244 1221 1214
Hi, All,
I try to accessing the x11() device over cgi - namely, when trying to
create a graphic using the jpeg() function, everything is fine from the
command line but it won't work over cgi, producing the error:
Unable to open connection to X11 display
Do you you know how to solve this problem
Dear all
I try to print 9 plots on a page, arranged as the code shows below.
nf - layout(matrix(c(1,0,2,0,0,3,0,4,0,5,0,6,0,0,0,0,7,0,8,9), 10,2))
layout.show(nf)
but when I try to plot, an error message
Fehler in plot.new() : Grafikränder zu groß
appears
to verify p.e. with
When you say you can not import 4.8GB, is this the size of the text
file that you are reading in? If so, what is the structure of the
file? How are you reading in the file ('read.table', 'scan', etc).
Do you really need all the data or can you work with a portion at a
time? If so, then
901 - 935 of 935 matches
Mail list logo