search for: dosage

Displaying 20 results from an estimated 22 matches for "dosage".

2010 May 10
1
how to get p-value from ave
...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 = aov(Alertness~Gender*Dosage,data=data.example2) > summary(aov.ex2)               Df  Sum Sq Mean Sq F value Pr(>F) Gender         1  76.562  76.562  2.9518 0.1115 Dosage         1   5.062   5.062  0.1952 0.6665 Gender:Dosage  1   0.063   0.063  0.0024  0.6665 Residuals     12 311.250  25.938               ...
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 -8.500 -2.437  0.250  2.687  8.500 Coefficients:             Estimate Std. Error t value Pr(>|t|)    (Intercept)   32.500      2.010  16.166 6.72e-11 *** Dosageb   ...
2009 Nov 08
2
Simple 2-Way Anova issue in R
...ing 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 : Factor w/ 2 levels "f","m": 2 2 2 2 2 2 2 2 1 1 ... $ Dosage : Factor w/ 2 levels "a","b": 1 1 1 1 2 2 2 2 1 1 ... $ Alertness : int 8 12 13 12 6 7 23 14 15 12 ... aov.ex2 = aov(Alertness~Gender*Dosage,data=data.ex2) summary(aov.ex2) Outputs: Df Sum Sq Mean Sq F value Pr(>F) Gender 1 76.562 76.562...
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 m a 13 4...
2012 Feb 05
4
nested if else statements
I have a vector of 2,1,0 I want to change to 0,1,2 respectively (the data is allele dosages) I have tried multiple nested if/else statements and looked at the ?if help and cannot work out what is wrong, other people have posted code which is identical and they state works. Any help would be greatly appreciated. > A[1:20] [1] 1 1 0 0 1 0 1 0 1 0 0 0 1 1 0 1 1 1 0 0 > B <- rep...
2023 Dec 18
3
Function with large nested list
...tions. 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 standard dosage and safety information in case of contraindications. If patient has an indication, update current dosage and safety information with the value from the dataframe row. 4. Save dosage and safety information in some lists and continue to the next drug. 5. When the iteration over all drugs is done, r...
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         8          1          158 1...
2011 Dec 13
0
snpStats imputed SNP probabilities
...e? 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 'dosage' as the output which is the expected number of risk alleles as a single value. However I need the 3 genotype probabilities (the probability of having 0, 1 and 2 risk allele copies) and not the 'dosage'. I couldn't find out how to obtain these probabilities. To illustrate this, he...
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
2006 Dec 21
0
Poisson mixed effects model
Hello, I am fitting a Poisson mixed effects model. I have the number of eggs (Eggs) laid by a quail and looking at the effect of dosage of a chemical (Dose) in the study. I have counts of eggs laid by week of the study, so I am including the week number (Week) as a random effect. I'm using the lme4 package. I have, > mod1 <- lmer(Eggs~Dose + (1|Week),family=poisson) > summary(mod1) Generalized linear mixed model fi...
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
2011 Mar 18
1
Help with setting the y-axis in a plot.
...t(lsd, type="o", col="blue", ylim=g_range, axes=FALSE, ann=FALSE) axis(1, at=1:3, lab=c("Low","Med","High")) axis(2, las=1, at=5*0:max(mar)) box() lines(mar, type="o", pch=22, lty=2, col="red") title(main="The Effect of Drug Dosage and Type on the Number of Hallucinations", col.main="red", font.main=4) title(xlab="Dose", col.lab=rgb(0,0.5,0)) title(ylab="Number of Hallucinations", col.lab=rgb(0,0.5,0)) Thanks, ~Caitlin [[alternative HTML version deleted]]
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 for me. Lower dosage means less nasal congestion and more feeling during . But 10 mg still leaves me rock hard and able to hold off climaxing just about as long a...
2011 Apr 14
1
integer and floating-point storage
...edigree format: http://pngu.mgh.harvard.edu/~purcell/plink/data.shtml#ped That might be too much to hope for. ;-) I think that the R system uses double-precision floating point numbers by default. When I impute minor-allele counts, I get posterior expected values ranging from 0 to 2 (called dosages). The imputation isn't very precise, so it would be fine to store such data using one or two bytes. (The values are used as regressors and small changes would have minimal impact on results.) I could use unsigned 8-bit integers (0 to 255), probably using only 0 to 254 so that 1 and 2 cou...
2003 Mar 13
0
Repeated measures 2-way anova -- robustness question
...ks 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 contrast tests looked at differences between the drugs, and between doses, and found p<.001 for cholesterol differences between one dosage block of drug1, the control block and the similarly-dosed block of drug2. Data on one of the drug blocks was actually just 3 dogs, because one dropped out. One critic has also claimed that the population variance of cholesterol levels in these dogs is so wide as to cast suspicion on estimates bas...
2011 Jun 07
0
ID50 i) comparisons ii) dose.p vs Reed-Muench
...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 above 50% infection! >SetA<-data.frame(c("Hi","Med","Lo"),c(8.44,7.46,6.22),c(31,27,35),c(21,14,6)) >names(SetA)<-c("Group","dose","total",&...
2004 Jul 21
1
RE: [Rcom-l] rcom 0.97 available for download
> From: Thomas Baier > > > I am confused as to what I need and must do. > > > > Your message said > > ftp://cran.r-project.org/pub/R/src/contrib/Devel/rcom_0.97.tar.gz > > > > The DESCRIPTION file inside that download says > > Package: rcom > > Version: 0.95 > > Date: 2004-07-14 > > That's too bad. Sorry for that. > >
2010 May 02
2
Calculation error
Dear Rxperts, Running the following code: ======================================================= twlo=10; twhi=20; wt=154; vd=0.5; cl=0.046; tau=6; t=3; F=1; wtkg <- wt/2.2 # convert lbs to kg vd.pt <- wtkg * vd # compute weight-based vd (L) cl.pt <- wtkg * cl # compute CL (L/hr) k <- cl.pt/vd.pt # compute k (hr^-1) cave <-
2005 Oct 24
1
Basic: setting resolution and size of an R graphic
...ic from the script below that would fit into one column of a published article and have a high resolution (as eps, or failing that tiff or png)? Thanks in advance for any advice, Peter ## Western.R ## A script to display the results of quantitative Western blotting with 6 repeats each at three dosages. ## This particular script has data from stimulation of fibroblasts with M-wt. # ---------- CONVENIENCE FUNCTIONS --------------------- # ## Define a simple function to draw the error bars. makeBars <- function(x,mean,se){ segments(x,mean - se/2,x,mean+se/2) segments(x-0.1,mean - se/2,x+0...
2013 Jan 26
0
Processed: switching email address
...xbackup-client: fails to purge - command ucf in postrm not found Changed Bug submitter to 'Andreas Beckmann <anbe at debian.org>' from 'Andreas Beckmann <debian at abeckmann.de>' > submitter 672153 ! Bug #672153 {Done: Andreas Beckmann <debian at abeckmann.de>} [dosage] dosage: unowned files after purge (policy 6.8, 10.8): /usr/lib/python2.7/dist-packages/dosage/plugins/dropin.cache Changed Bug submitter to 'Andreas Beckmann <anbe at debian.org>' from 'Andreas Beckmann <debian at abeckmann.de>' > submitter 689108 ! Bug #689108 {Done...