search for: dosages

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 =

...andard error: 4.924 on 15 degrees of freedom Multiple R-squared: 0.5396,     Adjusted R-squared: 0.4782 F-statistic: 8.789 on 2 and 15 DF,  p-value: 0.002977 As far as I understand it the lines "Dosageb" and "DosageC" represent the difference between DosageA and the other two dosages. My question is this: are these differences and the p-values associated with them the same as a t.test or pairwise.t.test on these groups? If I do t.tests, I get different values for t and p-value from those in the anova table above. Can someone please explain what the discrepancy is? Thanks...

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 :

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

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(...

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

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        

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

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

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) >

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

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,

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

...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 coul...

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

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

> 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. > >

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 <-

...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....

Processing commands for control at bugs.debian.org: > # bugs with submitter debian at abeckmann.de > submitter 479445 ! Bug #479445 [gnuplot] gnuplot: wrong line numbers for missing files reported during replot Changed Bug submitter to 'Andreas Beckmann <anbe at debian.org>' from 'Andreas Beckmann <debian at abeckmann.de>' > submitter 681568 ! Bug #681568