similar to: Poisson mixed effects model

Displaying 20 results from an estimated 200 matches similar to: "Poisson mixed effects model"

2010 May 10
1
how to get p-value from ave
Hi there, I checked google for aov. usually one uses summary to see whether the p-value is small. but I want to put aov in my script. how can I get the p-value, (0.1115, 0.6665, 0.6665 in the following example)? thanks YU > datafilename="http://personality-project.org/r/datasets/R.appendix2.data" > data.example2=read.table(datafilename,header=T) > aov.ex2 =
2009 Nov 08
2
Simple 2-Way Anova issue in R
Hello, I'm new to R and have been following many guides including the two-way anova (http://www.personality-project.org/r/r.anova.html). Using that walkthrough including the supplied data I do str(data.ex2) and receive the appropriate types of data as follows: > str(data.ex2) 'data.frame': 16 obs. of 4 variables: $ Observation: int 1 2 3 4 5 6 7 8 9 10 ... $ Gender :
2008 Dec 12
0
repeated two ways anova
Hi there, I have a question regarding a data set that looks like the following: 8 repeated measurements for 3 dosage levels. However, the measurements for each dosage levels are taken on a different group of subjects. example: subject time dosage measure 11         1          1          154 11         2          1          150 :          :           :             : 11        
2011 Dec 13
0
snpStats imputed SNP probabilities
Hi, Does anybody know how to obtain the imputed SNP genotype probabilities from the snpStats package? I am interested in using an imputation method implemented in R to be further used in a simulation study context. I have found the snpStats package that seems to contain suitable functions to do so. As far as I could find out from the package vignette examples and its help, it gives the
2010 Sep 20
2
interpreting one-way anova tables
Hi, I am trying to reconcile anova table in R (summary(lm)) with individual t.test. datafilename="http://personality-project.org/R/datasets/R.appendix1.data" data.ex1=read.table(datafilename,header=T) #read the data into a table summary(lm(Alertness~Dosage,data=data.ex1)) gives: Call: lm(formula = Alertness ~ Dosage, data = data.ex1) Residuals:    Min     1Q Median     3Q    Max
2004 Sep 28
0
My testimonial about skuper viakgra
I'm 48, and I can't believe I'm having the best of my life now thanks to cliallis Soft Tabs. The best part is no longer worrying about whether I'll be hard. With this stuff, I am always stiff, without fail. You have to experiment with the dosage and time lag. I started with 20 mg (good to rev my engine, I think) but ultimately found that 10 mg about 15 minutes before is ideal
2012 Oct 23
2
Data type in a data frame
Hi all, How does R know to regard a variable as a factor and not a character? For example, consider the following table: Observation Gender Dosage Alertness 1 m a 8 2 m a 12 3
2023 Dec 18
3
Function with large nested list
Hello list, I want to make a large rulebased algorithm, to provide decision support for drug prescriptions. I have defined the algorithm in a function, with a for loop and many if statements. The structure should be as follows: 1. Iterate over a list of drug names. For each drug: 2. Get some drug related data (external dataset). Row of a dataframe. 3. Check if adaptions should be made to
2011 Feb 28
1
Logwatch reporting spamassassin messages as unmatched entries
I've recently switched to using spamassassin via a sendmail milter, rather than using procmail to invoke it. This means that I get a number of messages appearing in my maillog, and then being reported by logwatch as unmatched entries. An example of such a messages is: Feb 27 04:33:09 quail sendmail[24780]: p1R4X46P024780[2]: URIBL blacklist\n\t* [URIs: tablettoxicspillsrx.ru]\n\t* 1.5
2019 Sep 19
5
deprecations leading up to C8
Was reading the long list of things being deprecated in the RHEL 8 release notes t'other day, and saw my old friend sendmail on the list. I've stuck by her because I've learned how she likes to be stroked to get the best out of her, and quail at the thought of taking on a new MTA and having to learn it all all over again. But I guess I'm going to have to, soon or late, probably
2011 Mar 18
1
Help with setting the y-axis in a plot.
Hi all. I'm working on an assignment for a psychology class and I am not sure how to adjust the y-axis so it displays the range: 0, 5, 10, 15 The code below is almost ideal: lsd = c(3, 5, 13) mar = c(1, 2, 3) g_range <- range(0, lsd, mar) plot(lsd, type="o", col="blue", ylim=g_range, axes=FALSE, ann=FALSE) axis(1, at=1:3,
2011 Jul 13
1
How to remove last 3 letters
Hello to everybody I have vector of drug dosage in string format 1.5 mg , 50 mg , 100 mg , 0.5 mg ... The task is to remove last free symbols: ' mg' Any suggestions are highly appreciated
2019 Sep 20
0
deprecations leading up to C8
Does that mean it's gone from the CentOS repos or just that PostFix is the default choice? On 9/19/19 6:04 PM, Fred Smith wrote: > Was reading the long list of things being deprecated in the RHEL 8 > release notes t'other day, and saw my old friend sendmail on the list. > > I've stuck by her because I've learned how she likes to be stroked to > get the best out of
2019 Sep 21
2
deprecations leading up to C8
On Fri, Sep 20, 2019 at 12:45:03PM -0500, SternData wrote: > Does that mean it's gone from the CentOS repos or just that PostFix is > the default choice? I believe it means that it won't be in Centos anymore. I also believe the Postifx is already the default choice, Sendmail isn't in the default installation. > > On 9/19/19 6:04 PM, Fred Smith wrote: > > Was
2019 Sep 20
0
deprecations leading up to C8
When I migrated our mx boxes from CentOS 5 to 7 I made the leap to postfix. I built a test server with postfix etc on it and threw everything I could think of at it before going live. The key is, like in learning a new language, to start thinking in postfix terms instead of thinking in sendmail terms. It takes a while to 'get it' but I'm glad we made the link. We search, read as
2003 Mar 13
0
Repeated measures 2-way anova -- robustness question
I'm a journalist, wondering what questions to ask about a study that contrasted the impact on serum cholesterol of two drugs. This was a 40 dog study: 5 treatment blocks of 4 dogs each, randomized to: a control block, two blocks at different doses of drug1 and two at corresponding doses of drug2. Analysis was 2-factor repeated measures ANOVA on treatment group and sampling time. Linear
2011 Jun 07
0
ID50 i) comparisons ii) dose.p vs Reed-Muench
I'm analysing some ID50 data for 2 different groups and had already calculated this by hand using Reed-Muench formula, when I came across the dose.p function in R. I have 2 queries: 1) dose.p gives me a different answer to Reed-Muench, and actually I suspect wrong answer, given that the dose.p result dosage stated to infect 50% is actually stronger than the dose used in my experiments caused
2012 Jul 03
3
design matrix creation in R
Hello, I want to create a design matrix using R. Can you explain the code which creates the following please? I understand the first part. b=g1(?) does what? dd <- data.frame(a = gl(3,4), b = gl(4,1,12)) # balanced 2-way dd a b 1 1 1 2 1 2 3 1 3 4 1 4 5 2 1 6 2 2 7 2 3 8 2 4 9 3 1 10 3 2 11 3 3 12 3 4 I am using the tree dataset in R. I want to form a reparameterized design
2012 Jul 25
0
Integrate: compound distribution
Dear All, thanks in advance for your attention. I'm struggling with this for a few weeks now, and figured it's time to get some help. I want to integrate(f(x), lower=-Inf, upper=Inf) with f(x) = ((gamma(K+1)/(gamma(r+1)*gamma(K-r+1)))*(q(x)^r)*(((1-q(x))^(K-r))*phi(x), where phi(x) is the standard normal pdf, and q(x) the logistic CDF (with inverse scale parameter, so a little
2007 Feb 05
1
manova discriminant functions?
Hello, I've been playing with the manova() function to do some pretty straightforward multivariate analyses, and I can't for the life of me figure out how to get at the discriminant functions used. When predicting several variables simultaneously, it's important to be able to gauge how much each variable is contributing to the analysis...a simple p-value isn't really enough. I