TOR TO HOLD ZK VALUES
>for(x in 1:nrow(xk)) {
> k <- intersect(i1[x,], i2[x,])
> zk[x] <- mean(unlist(k), na.rm = TRUE)
>}
>xk$zk <- zk
>data <- na.omit(xk)
>--
>View this message in context:
>http://www.nabble.com/memory.
search that part of the archive for "relative risk" you will find
a number of suggestions. Unfortunately the responses are not all
threaded so you need to search the whole thing.
>Tony.
> [[alternative HTML version deleted]]
Michael Dewey
h
clues.
>Thanks in advance..
>
>Michael
>
>
Michael Dewey
http://www.aghmed.fsnet.co.uk
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org
our space
bar. Computer storage is cheap.
>Thanks, Elva.
>
>_
>Got a favourite clothes shop, bar or restaurant? Share your local knowledge
>
>
Michael Dewey
http://www.aghmed.fsnet.co.uk
_
optim(start, fmin, gmin, method = "BFGS", hessian = Hess, ...) :
>initial value in 'vmmin' is not finite
>
>Something to do with the optimitation procedure... but, how can I fix
>it? Any help would be greatly appreciated.
>
>Thanks.
Michael Dewey
sufficient to use one of the results of
RSiteSearch("random multivariate normal", restrict = "functions")
or have I completely misunderstood what you want? (I appreciate this
is not exactly what you say you want.)
>Thanks for the help.
>
>Greg
>
>
>
ge,digits=digit)
> windowrange <- seq(from=lowppm,to=highppm,by=10^(-digit))
> if (sum(round(windowrange,digits=digit) %in%
>round(presence,digits=digit)) < 1) {
> newframe <- rbind(newframe,frame[int,])
> presence <- c(presence,windowrange)
> }
>
t;In addition: Warning messages:
>1: NaNs produced in: log(x)
>2: fitted values close to 0 or 1 in: tfun(mu = mu, y = y, w = w, res =
>FALSE, eta = eta, extra)
>3: 19 elements replaced by 1.819e-12 in: checkwz(wz, M = M, trace = trace,
>wzeps = control$wzepsilon)
>
>I will a
com/Difference-in-linear-regression-results-for-Stata-and-R-tf4069072.html#a11563283
>Sent from the R help mailing list archive at Nabble.com.
Michael Dewey
http://www.aghmed.fsnet.co.uk
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mail
thz.ch 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.
>
>
>--
>"The scientists of today think deeply
anybody knows a elegant way to solve the problem?
If I understand you correctly
?merge
should help you here
>Thanks in advance
>
>Thomas H.
>
>
Michael Dewey
http://www.aghmed.fsnet.co.uk
__
R-help@stat.math.ethz.ch mailing list
https://
receive
the emails in MIME format and the responses do get properly threaded
as far as I can see from the mailing list archives.
>Thank you!
>
>Moshe Olshansky
>[EMAIL PROTECTED]
Michael Dewey
http://www.aghmed.fsnet.co.uk
__
R-help@stat.m
91611504 3.455925e-10 2.089893e-01 5.531904e-02
> 4
>36 10171895 0.0089116669 2.956742e-10 2.153844e-01 8.614259e-02
> 4
>42 10198295 0.0078515166 3.147140e-10 2.437943e-01 7.314111e-02
> 5
>
>
>Thanks
>
>Rina
>
>
>
>
>
but is MatrixExp in the msm
package (available from CRAN) not relevant here?
[snip other suggestion]
Michael Dewey
http://www.aghmed.fsnet.co.uk
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the po
r VM can easily deal with))
>which will most probably just appear as one long e-mail in
>Thunderbird, with "table of contents" of all the subject lines,
>but nothing "clickable"
>
>Regards,
>Martin
>
>Martin Maechler, ETH Zurich
Michael Dewey
[EMAIL PROTECT
chive/42884.html
a help
>Cheers
>Robert Biddle
>
>
Michael Dewey
http://www.aghmed.fsnet.co.uk
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-gu
ber 2004)
>the function can be called by:
>
> temp(Sat ~ Infl + Type + Cont, housing)
>
>where all data is available from MASS, as it is
>an example in R Help on 'polr'.
>
>Results are:
>
> Re-fitting to get Hessian
>
>Error in eval(e
gt;Someone knows how I can solve the problem?
By giving us a short reproducible example?
Specifically we do not know:
1 - what formula1 is
2 - what the structure of data1 is
3 - what the enclosing function looks like
4 - what dat is
>Thanks for any help.
>Paolo
Michael Dewey
http://www.
> > __
> > R-help@stat.math.ethz.ch 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, min
rmance is not,
>for me in this case, important. Rather, I
>seek an easy R solution to find the time windows
>(starting & ending indicies) where 5 or more
>small numbers in my search_range were measured
>all in a row.
>
>Advice welcome and many thanks in advance.
>
>Ed
interval:
> 7.481272 11.734136
>sample estimates:
>odds ratio
>9.42561
>
>References
>
>Eronen, M. et al. (1996 - 1) Mental disorders and homicidal behavior
>in Finland. Archives of General Psychiatry, 53, 497-501
>
>Eronen, M et al (1996 - 2). Schizophrenia &am
tion (R's is based on the noncentral
>hypergeometric distribution), but the classical asymptotic formula gives
>
> > exp(log(41*2936210/920/20068)+qnorm(c(.025,.975))*sqrt(sum(1/M)))
>[1] 4.767384 8.918216
>
>which is closer, but still a bit narrower.
>
Michael Dewey
ht
to 0.5
>95 percent confidence interval:
> 0.03115856 0.05795744
>sample estimates:
> p
>0.04266389
>
> >
> > ci.p <- prop.test(920, 961)$conf
> > ci.odds <- ci.p/(1-ci.p)
> > ci.odds
>[1] 16.25404 31.09390
>attr(,"conf.level"
ent:
>>This email message, including any attachments, is for the so...{{dropped}}
>>
>>__
>>R-help@stat.math.ethz.ch mailing list
>>https://stat.ethz.ch/mailman/listinfo/r-help
>>PLEASE do read the posting guide http://ww
set parameter(s)
>basicstyle=\small, % print whole in small
>stringstyle=\ttfamily, % typewriter type for strings
>numbers=left, % numbers on the left
>numberstyle=\tiny, % Tiny numbers
>stepnumber=2, % number every second line of code
>numbersep=5pt, % 5pt seperation between number
translate your example (traddutore, traditore and all that)
>Thank you
>
>Massimo Cressoni
>
>run;
Michael Dewey
http://www.aghmed.fsnet.co.uk
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEAS
e manual pages which usually describe what is rather
than tell you explicitly what to do to achieve salvation. Perhaps I
could suggest that the if page tells you what error message you get
when you wanted ifelse and that the ifelse page has a health warning
along the lines of
stent
>
>
>
>Büro:
>
>IREBS Immobilienakademie
>
>c/o ebs Immobilienakademie GmbH
>
>Berliner Str. 26a
>
>13507 Berlin
>
>
>
>Privat:
>
>Dunckerstraße 60
>
>10439 Berlin
>
>
>
>Mobil: 0176 7002 6687
>
>Email: <mailto:
an Anselmo, CA 94960-2674
>tel 415.407.1357
>fax 415.485.1065
>email: [EMAIL PROTECTED]
>web:www.pinest.org
> www.pinestreetfoundation.org
> www.medepi.net/meta
>
>
Michael Dewey
http://www.aghmed.fsnet.co.uk
__
of your analysis
>thx
>
>ej
>
>On 12/3/06, Michael Dewey <[EMAIL PROTECTED]> wrote:
>>At 05:55 03/12/2006, Ethan Johnsons wrote:
>> >I have the clinical study data.
>> >
>> >Year 0 Year 3
>> &g
e estimates:
> prop 1prop 2
>0.8687259 1.000
>
>Warning message:
>Chi-squared approximation may be incorrect in: prop.test(c(2.25,
>1.53), c(2.59, 1.53))
>
>Any ideas, please?
>
>thx
>
>ej
>
>
Michael Dewey
http://www.aghmed.fsnet.co.uk
___
re file and then
input your data. It also enables you to implement checks during data
entry which in my view is an important feature. It enables export in
a number of formats but its native output format, the .rec file, can
be read into R using the foreign package.
>Duncan Murdoch
.19098 -1.225 0.221718
>
>SST NA NA NA NA
>
>
>
>Can someone explain the output "Coefficients: (1 not defined because of
>singularities)"? What does this mean? I'm guessing that something to do
>with my SST da
,data=p5,family=binomial,control=list(usePQL=FALSE,msV=1))
>>> > p5.random.p1 <- lmer(Y
>>>~1+(1|p),data=p5,family=binomial,control=list(usePQL=FALSE,msV=1))
>>
>>Well, try them and see. Then if you cannot understand the output tell us
>>a) what you found
model. Is this a cell means model?
>m00 <- glm(sc ~aa-1,data = p5)
>m000 <- glm(sc ~1+aa-1,data = p5)
See above
>thanks,
>
>Aimin Yan
>
>
Michael Dewey
http://www.aghmed.fsnet.co.uk
__
R-help@stat.math.ethz.ch mailing li
UE)
do what you wanted.
>Thanks,
>
>Aimin
>
>
Michael Dewey
http://www.aghmed.fsnet.co.uk
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/
se the studlab parameter what happens, and how does
that differ from what you wanted to happen?
Michael Dewey
http://www.aghmed.fsnet.co.uk
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guid
At 14:43 10/11/2006, Michael Dewey wrote:
After considerable help from list members and some digging of my own
I have prepared a summary of the findings which I have posted (see
link below). Broadly there were four suggestions
1 - Wald-type intervals,
2 - transforming the odds ratio confidence
t parameters
>
>Of course if someone have better idea to import data - not trough
>csv file no problem any format is OK - data frame or something else
>and to make forest plot ... no problem ... I need them on the forest plot.
>Greetings...
>Peter Boolton
>MCA adcc..
Michael
At 12:35 12/11/2006, Peter Dalgaard wrote:
>Michael Dewey <[EMAIL PROTECTED]> writes:
>
> > At 15:54 10/11/2006, Viechtbauer Wolfgang (STAT) wrote:
> >
> > Thanks for the suggestion Wolfgang, but whatever the original authors
> > did that is not it.
>
http://www.wvbauer.com/
>
>
> > -Original Message-
> > From: [EMAIL PROTECTED] [mailto:r-help-
> > [EMAIL PROTECTED] On Behalf Of Michael Dewey
> > Sent: Friday, November 10, 2006 15:43
> > To: r-help@stat.math.ethz.ch
> > Subject: [R] Confidence in
authors of the
paper I am refereeing that
I think they are wrong unless I can help them with what they
should have done.
Is there another package I should have tried?
Is there some other way of doing this?
Am I doing something fundamentally wrong-headed?
Michael Dewey
http://www.aghmed.fsnet.
cond issue, I would also like to offset the second set of points
>along the x-axis so that the error bars will be visible.
>
>Any help is appreciated.
>
>Regards,
>Sumit
>
>
Michael Dewey
http://www.aghmed.fsnet.co.uk
__
R-help@stat.m
t
>Thank you for your help.
>
>Karen
>
> [[alternative HTML version deleted]]
Michael Dewey
http://www.aghmed.fsnet.co.uk
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the p
2226.0825 8.7
>2.6 9 10 0 2226.0825 8.7
>2.7 10 11 0 2226.0825 8.7
>2.8 11 12 0 2226.0825 8.7
>2.9 12 13 0 2226.0825 8.7
>2.1 13 14 0 2226.0825
> > table <- read.table('latestWithNumber.txt', header=T)
>Error in "row.names<-.data.frame"(`*tmp*`, value = c("1.1", "1.2", "1.3", :
> duplicate 'row.names' are not allowed
>
>Yongchuan
Michael Dewey
http://
to. Biologia Geral/Lab. Ecologia Evolutiva
>| : :' : Campus Universitário Prof. Darcy Ribeiro, Vila Mauricéia
>| `. `'` CP: 126, CEP: 39401-089, Montes Claros - MG - Brasil
>| `- Fone: (38) 3229-8190 | [EMAIL PROTECTED]
>| ICQ#: 5692561 | LinuxUser#: 205366
Michael Dewey
h
an(link=log))
>Error in eval(expr, envir, enclos) : cannot find valid starting
>values: please
>specify some
Have you tried specifying some starting values as it asks?
Try start=c(-3,1) for instance
Michael Dewey
http://www.aghmed.fsnet.co.uk
_
orm?
What version of mvtnorm do you have?
What version of R?
>Thanks a lot.
>
>Yonghai Li
Michael Dewey
http://www.aghmed.fsnet.co.uk
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read th
ing to answer, but I find that when
I am tempted to do this it is almost always better to make a list (in
your case of data frames) and then access the elements of the list
>Example:
>week <- 28
>test(week) <- data.frame(a,b,s,c);
>test28
Michael Dewey
http:
t;
>CPData Optimum
And please do not send HTML.
> [[alternative HTML version deleted]]
Michael Dewey
http://www.aghmed.fsnet.co.uk
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
v1v2v3v4
> v1v2v3v4
> etc...
>
>Any advice will be appreciated.
>
>Regards
>Michael.
>
>
>
>
>
>________
>
>Coming soon: Celebrity Survivor - 11 celebrities, 25 days, unlimit
>D-20146 Hamburg
>Germany
>
>phone: +49-(0)40-42838.7498 (office)
>+49-(0)40-42838.4591 (Billon)
>fax: +49-(0)40-42838.2344
>email: [EMAIL PROTECTED]
>www:
>http://www2.jura.uni-hamburg.de/instkrim/kriminologie/Mitarbeiter/Enzmann/Enzmann.html
>
>
Michael
about other packages with the same aim and why I should use this one
tell me where to go for more information
>but great vignettes, websites or book are also of interest.
>
>What is your favourite bit of R documentation, and why?
>
>Thanks,
>
>Hadl
; > 100 8th Avenue S.E.
> > St. Petersburg, FL 33701
> > (727)896-8626 Ext. 4305
> > [EMAIL PROTECTED]
> >
> > __
> > R-help@stat.math.ethz.ch mailing list
> > https://stat.ethz.ch/mailman/listinfo/r-help
> >
that I'm doing something wrong. Must weights be integers, or at
>least greater than one?
>
>Regards,
>
>--
>Ramón Casero Cañas
>
>http://www.robots.ox.ac.uk/~rcasero/wiki
>http://www.robots.ox.ac.uk/~rcasero/blog
Michael Dewey
http://www.aghmed.fsnet.co.uk
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
on R and R-help) that if I type cor at the prompt what I see looks
like source code but is not _the_ source code.
Michael Dewey
http://www.aghmed.fsnet.co.uk
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLE
g and reading at least one of the many books
on R. When I started I found the free documents useful but I made most
progress when I bought MASS. I do realise that liking books is a bit last
millennium.
Michael Dewey
http://www.aghmed.fsnet.co.uk
"&", "\&", temp)
[1] "R & D"
> sub("&", "\\&", temp)
[1] "R & D"
> sub("&", "\\\&", temp)
[1] "R & D"
> sub("&", "&", temp)
[1] &
ossible, considering NAs as a level).
If they are ordered have you considered downloading the polycor package?
Michael Dewey
[EMAIL PROTECTED]
http://www.aghmed.fsnet.co.uk/home.html
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mai
sperando una respuesta desde la universidad del valle, quien
> le escribe. Andres Felipe
>
>
>
>-
> 1GB gratis, Antivirus y Antispam
> Correo Yahoo!, el mejor correo web del mundo
> Abrí tu cuenta aquí
> [[alter
fill in the details in the same way as I have done for Eudora.
The draft is available at
http://www.aghmed.fsnet.co.uk/r/digest.htm
Any other comments welcome of course.
Michael Dewey
http://www.aghmed.fsnet.co.uk
__
R-help@stat.math.ethz.ch mailing
ect of the message.
Hope that helps.
>Regards,
>Martin
Michael Dewey
http://www.aghmed.fsnet.co.uk
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
ntes:
>
>ImportationPays Produit Tonnage
> Côte d'IvoireCoconut10123
>
>... et je n'y arrive pas. Ce doit pourtant être basique.
>
>Merci pour votre aide.
>
>Jacques Vernin
>PiR2
>10, Boulevard de Brazza
>13008 M
tion does I lack
the skill to modify it myself. I will have to wait for someone more gifted.
If there is to be a successor edition to Pinheiro and Bates perhaps I could
suggest that this topic merits a bit more discussion?
Michael Dewey
[EMAIL PROTECTED]
http://www.aghmed.fsnet.co.uk/home.html
2.1.0 under XP (SP2) with the versions of boot and nlme
which shipped with the binary.
Michael Dewey
[EMAIL PROTECTED]
http://www.aghmed.fsnet.co.uk/home.html
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLE
the combination of n things taking k at time.
How do I do that in R ?
I wrote my own function but in this form isn't easy for my
students.
Can you help me ?
Michael Dewey
[EMAIL PROTECTED]
http://www.aghmed.fsnet.co.uk/home.html
__
R-help@stat.math.et
code?
(d) I appreciate this last item is not strictly an R question, but my
intention is to use these as input into the sem package for structural
equation models. If anyone thinks that is misguided I would be intersted to
here.
Michael Dewey
[EMAIL PROTECTED
At 17:58 01/10/04, Michael Dewey wrote:
I think it is easiest to describe
what I want in terms of the concrete
problem I have.
I have data from a number of countries
in each of which a sample of people was
interviewed. In presenting the results
in a forthcoming collaborative publication
much
crude,
as long as Italy looks roughly like
a boot that is fine. I am an epidemiologist
not a geographer.
Michael Dewey
[EMAIL PROTECTED]
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
3 6 3
> B
a id
1 1 4
2 2 5
3 3 6
> AandB
a id b
1 1 1 4
2 1 4 NA
3 2 2 5
4 2 5 NA
5 3 3 6
6 3 6 NA
>
Give you what you wanted? You can always strip out the "id" column if you
wish (after sorting on it if you would like the original order back.
Thank you for your h
ion 1.
What exactly do you mean by question 2? Do you mean it is not reading the
hours column (Hora) or do you mean it takes a long time? It looks to me as
though column one, despite its label of day (Día) actually contains some
much more complex Excel format.
Michael Dewey
[EMAIL PROTECTED
e in R
Message: 56
Hello, it want to know if there is a library that transform a file in
excell (*. xls) to text file(*.txt, *,prn, *.csv). Thanks
library(RODBC) will read your Excel file into R. You do not need a copy of
Excel to do this.
Michael Dewey
[EMAIL PRO
ttice available. I
find myself making so many more graphical displays because they are (a)
easy (b) aesthetically pleasing (c) revealing about my dataset.
Deepayan
Michael Dewey
[EMAIL PROTECTED]
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.et
{Springer-Verlag}
}
Michael Dewey
[EMAIL PROTECTED]
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
(few thousand cases, few
hundred variables). In order to facilitate working with colleagues I need
to stick with windows even if linux would be more efficient
Michael Dewey
[EMAIL PROTECTED]
http://www.aghmed.fsnet.co.uk/home.html
__
[EMAIL PROTECTED
76 matches
Mail list logo