, I'm not
keen on using a GLM with the types predictor variables I have...
mobForest is a take-off on mob, but my initial testing indicated that it
took a VERY long time to fit the models.
Thanks!
Craig
On Thu, Aug 7, 2014 at 8:10 PM, Mitchell Maltenfort
wrote:
> I don't kn
Any ideas or suggestions greatly appreciated.
Cheers!
Craig
[[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/posti
should
look for that may potentially cause some issues?
>
> -Oorspronkelijk bericht-
> Van: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
> Namens Craig O'Connell
> Verzonden: maandag 28 april 2014 3:20
> Aan: r-help@r-project.org
> Onde
I am using a mixed model to assess the effects of various variables (i.e.
treatment, density, visibility) on bee behavior (e.g., avoidance frequency -
total avoidances per total visits; feeding frequency, and mating frequency).
Bee individuals is my random factor (n=63 different bees), whereas
I have been trying to process a data.frame (tr) on a grouping factor (Rmbr)
using diff
to find differences between srv_dt (a date) by Rmbr and keep getting the
following error message
> tmp <- by(tr,tr$Rmbr, diff, tr$srv_dt, lag=1, difference=1)
Error in r[i1] - r[-length(r):-(length(r) - lag +
t;- length(mbr.list)
> code.mat <- matrix(0,n.mbr,n.code)
> dimnames(code.mat) <- list(mbr.list,code.list)
> for(i in 1:n.mbr){
+ mbr.i <- claims.df[claims.df$member==mbr.list[i],] #EXTRACT BLOCK OF
MEMBER CLAIMS
+ code.mat[i,unique(match(mbr.i$code,code.list))] <
0 1 0 0 1
3 1 1 0 1 0
4 1 0 0 0 0
Craig McKinstry
100 Market, 6th floor
Office: 503-225-6878 | Cell: 509-778-2438
IMPORTANT NOTICE: This communication, including any attachment, contains
information that may be confid
I just received and email post with an R question. My problem is as described
above. Thanks for your help
Craig McKinstry
Cambia HealthSparq
100 Market Building, 6th Floor
100 SW Market Street
Portland, OR 97201
Ph: 503 225-6878
Email:
craig.mckins...@cambiahealth.com<mailto:craig.mck
hanks for any help given,
Peter Craig
--
##
| E-mail: p.s.cr...@durham.ac.uk Telephone: +44-191-3343076 (Work) |
| Fax: +44-191-3343051 +44-191-3737415
I am having a strange problem. I have been using python + rpy2 on my Mac
for quite sometime now. I recently upgraded to OS X 10.8 from 10.6 and
that's when the problem began.
Now, when I try to load some R packages (like zoo) using robjects.r, python
crashes with seg fault 11. Here's an example be
. I have to do
this for five scales with each scale having 6 items half of which may have
the 1 response pattern and half the 2 pattern.
Any suggestions on how best to tackle this in R would be helpful.
Craig
**
Old SPSS code sample
**
Compute Variable = 0.
IF (item1 = 1) Variable
Hello,
My name is Craig and I need help plotting a "line" for a multiple linear
regression in R.
Here is my sample data (filename: convis.txt)
Output of convis.txt is (vis and density being predictors of either
avoidance or entrance):
vis den avoid entrance
1 10 1 0.
:Â treatmentcon, feedingchum,
treatmentcon:feedingnofeed, treatmentcon:feedingsat, treatmentmag:feedingchum,
treatmentproc:feedingchum, treatmentcon:feedingchum ?  I have a strange
feeling that the issue is arising from my "family=binomial" term. Any
feedback would be great
Dear all,
  I am having a recurring problem when I attempt to conduct a GLM. Here is
what I am attempting (with fake data):
First, I created a txt file, changed the directory in R (to the proper folder
containing the file) and loaded the file:
#avoid<-read.table("avoid.txt",header=TRU
> I think this has already been done in r59493ff of the svn repos by the
> indefatigable Prof. Ripley
>
Fantastic. I guess it *is* true: the first step in writing any patch
for R is to first ask yourself if Brian Ripley has already made the
change since the last release.
Thanks Brian!
--
-cc
Hello R Community,
I don't post to these things often so excuse me if I stumble on my forum
etiquette. This is a complex problem for me, which may require two forum
entries, but I will try my best to be concise. Also, I am a self taught
coder, so if my code is not to convention, your constructi
Hello.
I am beginning to analyze my work and have realized that a simple chi-square
analysis will not suffice for my research, with one notable reason is that data
are not discrete. Since my data fit the assumptions of a logistic regression,
I am moving forward with this analysis. With t
g GLM by specifying my
own quasi family and describing the variance function using varfun, validmu,
dev.resids and initiate. However I cannot find any detailed discussion of
this method. Any help would be appreciated.
ps I've done this in SAS so I know it should be pos
Hello all,
I am trying to produce the closed-form Ait-Sahalia approximation to the
log-likelihood function of a diffusion process but the function arguments I
am passing to the function result in NaN being returned. I've checked the
'Transform','S' and 'M' functions and the problem doesn't seem to
Hello All,
I am trying to estimate the parameters of a stochastic differential equation
(SDE) using quasi-maximum likelihood methods but I am having trouble with
the 'optim' function that I am using to optimise the log-likelihood
function.
After simulating the SDE I generated samples of the simul
I am applying a GLS model and would like to look at diagnostic plots of
influence.
The function (plot(model)) that works for linear models does not seem to
function for GLS models. Is there a reason for this? Or is different code
required?
Sorry if it's a very basic question, but thanks for you
Dear R-users,
I am currently modifying a previously developed predator prey model and
was curious if there was a way to add in a disturbance to the model (let's
say at time t=100). The disturbance can be the introduction of 40 prey
(N=40) and 10 predators (Pred = 10). I would like to see my m
he TeachingDemos package
>
>
> --
> Gregory (Greg) L. Snow Ph.D.
> Statistical Data Center
> Intermountain Healthcare
> greg.s...@imail.org
> 801.408.8111
>
>
> > -Original Message-
> > From: r-help-boun...@r-project.org [mailto:r-help-boun...@
Hello R-list,
I am plotting a weighted linear regression in R. The points on my chart are
also scaled to sample size, so some points are large, some are small. I have
figured out everything I need using the plot() function: how to plot the
points, scale them by sample size, weight the linear regre
for a simple
> reproducible example? At the very least, we need your
> merge and sort code, since it sounds like that's where
> the problem truly lies.
>
> Sarah
>
> On Wed, Nov 10, 2010 at 2:02 PM, Ian Craig wrote:
> > So here is a sample of the datasets. I have
0
0 0 2009-07-23 11:15:00 141
15 38-C 217276 12/21/2008 4:31 12/22/2008 16:31 12/21/2008 4:31 0
0 0 2008-12-22 16:31:00 151
On Wed, Nov 10, 2010 at 1:12 PM, Sarah Goslee wrote:
> On Wed, Nov 10, 2010 at 12:57 PM, Ian Craig wrote:
> > Greeting
Greetings Supreme Council of R Masters,
Like toddler, I have gotten my head stuck in the banisters of R ... again.
Let it be know I am still a neophyte in the R-community forum world, so
please don't flame me too bad.
I have two sets of data, each with a set of timestamps. I would like to
someh
, as well as the p-value is given in SIR.
I was looking into running QDA and was wondering if there is a command
that will output a prediction equation and the constants for this
equation as well as a p-value of the separation?
Thanks for anyone who can help.
Craig
Daniel,
That's it! Thanks. Your help is very much appreciated. I'm hoping to
nail down the code conversion from MATLAB to R, but it seems to be a bit more
difficult that I had anticipated.
Craig
> From: djnordl...@frontier.com
> To: djnordl...@frontier.com; r-he
f the previous syntax in addition to changing the result
statement?
Thank you so much and I think my brain is fried! Happy Holiday.
Craig
> Date: Mon, 11 Oct 2010 09:59:17 +0200
> From: alain.guil...@uclouvain.be
> To: craigpoconn...@hotmail.com
> CC: pda...@gmail.com; r-help@r-proj
Thank you Peter. That is very much helpful. If you don't mind, I continued
running the code to attempt to get my answer and I continue to get inf inf
inf... (printed around 100 times).
Any assistance with this issue. Here is my code (including your corrections):
myquadrature<-functio
I need to find the area under a trapezoid for a research-related project. I
was able to find the area under the trapezoid in MATLAB using the code:
function [int] = myquadrature(f,a,b)
% user-defined quadrature function
% integrate data f from x=a to x=b assuming f is equally spaced over the
i
Please unsubscribe me from this list. Thank you.
<><
Marine Biologist
Elasmobranch Bycatch Reduction Scientist
SharkDefense Technologies, LLC
(845) 702-7087
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do re
Dear R Users,
I've never used R before and my professor has asked us to do some
pretty intense programming (or it's intense to me at least). Here is
the question: Modify the function myquadrature inside the script so
that it returns the quadrature of descrete data using the trapezoidal
rule
pen the cubic spline that is being used, or if there is an
alternative to interp() that I should be looking at. Like I said before, I'm
very new to R, so I may have missed something entirely obvious.
Thanks in advance,
Craig
__
R-help@r-pr
t; varespec
Error: object 'varespec' not found
I suspect this could be a way I have my files organised or that I have
not set a working directory? I am new to R so I suspect I have probably
made a simple mistake.
Tha
Be sure to check $estim.prec and use the parameter tol as needed.
Best regards,
Craig
Madhavi Bhave wrote:
Dear Sir,
That was GREAT!!!. Thanks a lot for the solution. Once again it showed how powerful 'R' is otherwise I was breaking my head on Newton-Raphson method.
Thanks
for solution in solve(eqn, ytm)][2]
0.0960070882662794
The third value is the only real solution.
There may be an R package for doing that, but I don't know one.
Best regards,
Craig
My query is (1) if there is any R package which will calcualte ytm or (2) is there any method in 'R&
Original Message
Subject:Re: [R] Suppress output from getGEO
Date: Fri, 29 Jan 2010 10:29:54 -0500
From: Duncan Murdoch
To: Craig P. Pyrame
References: <4b62f906.90...@gmail.com> <4b62fb32.4040...@stats.uwo.ca>
<4b62fe12.7040...@gmail.com>
#x27;s fault in this case.
Besides, I think you might want to make efforts to be less aggressive. I
don't think my use of words in my previous message was inappropriate,
but yours seems to be. I apologize if my post hurt your feelings, and
now expect you do the same,
Best regards,
Cr
noying output be modified to put an end to this?
Best regards,
Craig
__
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 commente
Dear Duncan and Rolf,
That's funny! Thanks a lot.
Best regards,
Craig
Duncan Murdoch wrote:
On 30/06/2009 5:11 PM, Craig P. Pyrame wrote:
Dear Rolf,
What do you mean?
He was talking about the fortunes package. Install it, type
fortune(), and you'll get a fortune cookie messa
Dear Rolf,
What do you mean?
Best regards,
Craig
Rolf Turner wrote:
On 1/07/2009, at 12:34 AM, Craig P. Pyrame wrote:
... it's probably not a good idea to submit bug
reports just because I misunderstand what R does.
Gotta be a fortune!!!
cheers,
Dear Barry,
Thank you for the suggestion, it does work. I think the documentation
might be improved, but it's probably not a good idea to submit bug
reports just because I misunderstand what R does.
Best regards,
Craig
Barry Rowlingson wrote:
On Tue, Jun 30, 2009 at 1:00 PM, Cr
Dear Henrique,
Thanks! This does work, and I find the following solution to my
original problem elegant enough:
> rep(list(character, integer, numeric, ...), c(3, 2, 2, ...))
Best regards,
Craig
Henrique Dallazuanna wrote:
From the help page of rep:
"Value:
An object of
one) is to use class names instead, like so:
> types = lapply(c(rep('character', 2), 'integer', 'numeric', ...),
function(type) vector(type, 0))
but I am still curious why the above doesn't work as I would expect it to.
Best regards,
Craig
_
Dear Peter,
Thank you for the explanations. Some further questions are inlined below.
Best regards,
Craig
Peter Dalgaard wrote:
Craig P. Pyrame wrote:
Stavros Macrakis wrote:
It gives me a headache, too! I think you'll have to wait for a more
expert user than me to supply explanatio
;ll try to learn more from the manuals.
Best regards,
Craig
-s
On 6/26/09, Craig P. Pyrame wrote:
Stavros Macrakis wrote:
On Thu, Jun 25, 2009 at 12:47 PM, Craig P. Pyrame
wrote:
The man page Stavros quotes states that the class attribute of the result
is
taken
Stavros Macrakis wrote:
On Thu, Jun 25, 2009 at 12:47 PM, Craig P. Pyrame wrote:
The man page Stavros quotes states that the class attribute of the result is
taken from 'test', which clearly is not the case:
Actually, the behavior is documented pretty clearly:
The m
Dear Rolf,
Rolf Turner wrote:
On 25/06/2009, at 12:27 PM, Craig P. Pyrame wrote:
Dear Stavros,
What you discuss below is somewhat scary to me as an R newbie. Is this
just an incident, a bug perhaps, or rather the way things typically go
in R, as your "Welcome to R!" seems to s
o anything with it!"
sort is gradually turning into an "I can't get why it doesn't work" and
"I can't get how to make this work" depression. I would be happy to
blame this on my incompetence and incapability, but would also like to
hear if it is not
on for the model would be the logit.
Thanks again, Professor Ripley, for your comments and suggestions.
Craig
Prof Brian Ripley wrote:
> glmmPQL can fit the same GLM families as glm() can -- it does not list
> _any_ .
>
> Howver, the beta distribution does not give a GLM family and h
inclusion of random effects.
Does anyone know of a package or code for a generalized linear mixed
model that allows a beta distribution? Transforming my data might allow
me to use another family, but I would rather not transform the data if
possible. Thanks for your help!
Sincer
be appreciated.
-Craig
__
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.
54 matches
Mail list logo
