Re: [R] Sum weights of independent variables across models (AIC)

2011-07-14 Thread Michael Just
This is what I was looking for. When I initially read about model.avg I didn't recognize it also provided variable scores. Thank you kindly, Mike __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

[R] Sum weights of independent variables across models (AIC)

2011-07-13 Thread Michael Just
Hello, I'd like to sum the weights of each independent variable across linear models that have been evaluated using AIC. For example: library(MuMIn) data(Cement) lm1 - lm(y ~ ., data = Cement) dd - dredge(lm1, beta = TRUE, eval = TRUE, rank = AICc) get.models(dd, subset =

Re: [R] Format of Output of Residuals

2010-10-04 Thread Michael Just
    -51.6454545 3    -43.6454545 4    -38.6454545 5    -26.6454545 ... HTH, Dennis On Thu, Sep 30, 2010 at 9:09 PM, Michael Just mgj...@gmail.com wrote: An excerpt from dataset ChickWeight:     weight Time Chick Diet 1   42    0 1    1 2   51    2 1    1 3   59    4

[R] Format of Output of Residuals

2010-09-30 Thread Michael Just
An excerpt from dataset ChickWeight:     weight Time Chick Diet 1   42    0 1    1 2   51    2 1    1 3   59    4 1    1 I am interested in the residuals of the dataset.  Specifically in saving them to another format. I have been creating text files with sink. CW.lm -

Re: [R] Adjusting Font Size: lattice / bwplot

2010-09-20 Thread Michael Just
Hello, I know I am supposed to get 'it' when I ?bwplot but your example really made things clearer for me. Thanks, Cheers, Michael On Mon, Sep 20, 2010 at 12:50 AM, David Winsemius dwinsem...@comcast.netwrote: On Sep 20, 2010, at 12:43 AM, Michael Just wrote: Hello, If you run

[R] Adjusting Font Size: lattice / bwplot

2010-09-19 Thread Michael Just
Hello, If you run the following code with lattice installed: bwplot(voice.part ~ height, data=singer, xlab=Height (inches)) There will be some text in the graph e.g. Soprano 1, Soprano 2, etc - 60, 65, etc and the title: Height (Inches). How can one make the font for this text larger and/or

[R] bargraph.ci - CI and color question.

2009-11-03 Thread Michael Just
(mean(x)-sd(x), mean(x) + sd(x)) ) Am I getting 95% CI bars? Thank you kindly, Michael Just __ 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

Re: [R] bargraph.ci - CI and color question.

2009-11-03 Thread Michael Just
appreciate an indication of where it is. Thank you kindly, Michael Just On Tue, Nov 3, 2009 at 8:53 AM, Manuel Morales manuel.a.mora...@williams.edu wrote: On Tue, 2009-11-03 at 03:51 -0600, Michael Just wrote: Hello, When using bargraph.CI in package sciplot can the bars for each group

Re: [R] bargraph.ci - CI and color question.

2009-11-03 Thread Michael Just
Thank you. On Tue, Nov 3, 2009 at 3:08 PM, Manuel Morales manuel.a.mora...@williams.edu wrote: On Tue, 2009-11-03 at 10:20 -0600, Michael Just wrote: Hello, Thanks for the responses.  Yes, I did try to use ?bargraph.CI for the colors. When I said bars,  I meant the main bars on the graph

[R] Select only cases with negative values

2008-10-18 Thread Michael Just
d for my regression, how could I make that selection? Thanks, Michael Just __ 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

Re: [R] Select only cases with negative values

2008-10-18 Thread Michael Just
, dat=dat, subset= c(d 0)) Error in eval(expr, envir, enclos) : invalid 'envir' argument Any other suggestions? Thanks, Michael Just On Sat, Oct 18, 2008 at 5:15 PM, Henrique Dallazuanna [EMAIL PROTECTED] wrote: Try this: dat.lm -lm(c~d, subset=d 0) On Sat, Oct 18, 2008 at 7:03 PM, Michael

Re: [R] Select only cases with negative values

2008-10-18 Thread Michael Just
PM, Johannes Hüsing [EMAIL PROTECTED] wrote: Am 18.10.2008 um 23:03 schrieb Michael Just: Hello, I was wondering if there was a way to only select cases my from data frame that contained a negative value? c-c(1,2,3,4,5,6,7,8,9,10) d- c(-1,2,-3,-4,5,6,-7,8,-9,10) f - cbind(c,d) dat

[R] Using source()

2008-10-16 Thread Michael Just
Hello, I have never used source and I am a R beginner. If I have text file that contains 1,000's of lines of code. Can I use source to bring this code into R and execute the code? Does it run the code one line at a time? Is there a best way to setup source() for maximum effieciency? After reading

Re: [R] Using source()

2008-10-16 Thread Michael Just
On Thu, Oct 16, 2008 at 1:42 AM, Prof Brian Ripley [EMAIL PROTECTED] wrote: On Thu, 16 Oct 2008, Michael Just wrote: Hello, I have never used source and I am a R beginner. Not a 'beginner' R-help poster, though. (32 posts this month so far under this name.) I hope this isn't bad form

[R] AIC score

2008-10-14 Thread Michael Just
Hello, I ran AIC for some competing models I created. I get df and an AIC score from the AIC procedure. Can I use the models with the lowest AIC scores from this procedure to choose my 'best' models? If not, what else do I need to do (and know) and how can I do it in R to chose the 'best' models?

[R] Maximum number of pasted 'code' lines?

2008-10-14 Thread Michael Just
Hello, I write most of my R code in excel and then paste it into R. I am wondering if there is a limit to how much I can paste? I want to paste about 19,000 lines of code should this work? I am doing this because when I did it chunks it took about an hour and half. I thought if I could insert it

Re: [R] Maximum number of pasted 'code' lines?

2008-10-14 Thread Michael Just
:35 PM, Erik Iverson [EMAIL PROTECTED]wrote: Michael Just wrote: Hello, I write most of my R code in excel and then paste it into R. Do you actually use Excel as a text editor? Is this common? What benefits do you get by writing code in Excel? [[alternative HTML version deleted

Re: [R] AIC score

2008-10-14 Thread Michael Just
that there is no simple answer to your question. John Frain 2008/10/14 Michael Just [EMAIL PROTECTED]: Hello, I ran AIC for some competing models I created. I get df and an AIC score from the AIC procedure. Can I use the models with the lowest AIC scores from this procedure to choose my 'best

[R] Add notes to sink output

2008-10-13 Thread Michael Just
Hello, How can I add notes (i.e. text) to a sink output? sink(test.txt) #This text will describe the test summary(x) sink() How can I add that text above to the sink output? Thanks, Michael [[alternative HTML version deleted]] __

Re: [R] Add notes to sink output

2008-10-13 Thread Michael Just
for more information. HTH, Jorge On Mon, Oct 13, 2008 at 4:02 PM, Michael Just [EMAIL PROTECTED] wrote: Hello, How can I add notes (i.e. text) to a sink output? sink(test.txt) #This text will describe the test summary(x) sink() How can I add that text above to the sink output? Thanks

[R] LM intercept

2008-10-13 Thread Michael Just
What is the difference when including or not including the intercept when using lm()? x.noint - lm(weight ~ group - 1))# omitting intercept x - lm(weight ~ group)) This has nothing to do with forcing the intercept to 0, correct? Thank you kindly, Michael [[alternative HTML version

Re: [R] LM intercept

2008-10-13 Thread Michael Just
Great, Thanks, Michael On Mon, Oct 13, 2008 at 3:56 PM, Rolf Turner [EMAIL PROTECTED]wrote: On 14/10/2008, at 9:42 AM, Michael Just wrote: What is the difference when including or not including the intercept when using lm()? x.noint - lm(weight ~ group - 1))# omitting intercept x - lm

Re: [R] LM intercept

2008-10-13 Thread Michael Just
change of what question is being answered. Hope this helps, From: [EMAIL PROTECTED] [EMAIL PROTECTED] On Behalf Of Michael Just [EMAIL PROTECTED] Sent: Monday, October 13, 2008 2:59 PM To: Rolf Turner Cc: r-help Subject: Re: [R] LM intercept Great

[R] Correlation among correlation matrices cor() - Interpretation

2008-10-10 Thread Michael Just
Hello, If I have two correlation matrices (e.g. one for each of two treatments) and then perform cor() on those two correlation matrices is this third correlation matrix interpreted as the correlation between the two treatments? In my sample below I would interpret that the treatments are 0.28

[R] Interpretation in cor()

2008-10-09 Thread Michael Just
])) Does this operation above tell me how correlated each of the three treatments is? Because this how I am interpreting it. Thanks, Michael Just [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

[R] Exporting symnum() result from cor()

2008-10-09 Thread Michael Just
them with non-R users? Is it possible to export the way things look in View() ? Thanks, Michael Just [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] Exporting symnum() result from cor()

2008-10-09 Thread Michael Just
? --- On Thu, 10/9/08, Michael Just [EMAIL PROTECTED] wrote: From: Michael Just [EMAIL PROTECTED] Subject: [R] Exporting symnum() result from cor() To: r-help r-help@r-project.org Received: Thursday, October 9, 2008, 5:04 AM Hello, I am trying to export the results from symnum() while

Re: [R] Interpretation in cor()

2008-10-09 Thread Michael Just
. Thanks, Michael On Thu, Oct 9, 2008 at 8:35 AM, John Kane [EMAIL PROTECTED] wrote: Please provide a small working example. It is difficult to see what you are doing from the description below. --- On Thu, 10/9/08, Michael Just [EMAIL PROTECTED] wrote: From: Michael Just [EMAIL PROTECTED

Re: [R] Exporting symnum() result from cor()

2008-10-09 Thread Michael Just
Richard, Thanks this works well. -Michael Just On Thu, Oct 9, 2008 at 2:43 PM, Michael Just [EMAIL PROTECTED] wrote: Hello John, Is it possible to write more than one data frame to a single write.csv or write.table? Thanks, Michael On Thu, Oct 9, 2008 at 8:38 AM, John Kane [EMAIL

[R] Write multiple objects / data frames into a single .csv -- write.csv

2008-10-09 Thread Michael Just
.LN.sl.vs.mean) sink() Could I put these all into a single .csv, in a similar format to sink()? Thanks, Michael Just [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] Write multiple objects / data frames into a single .csv -- write.csv

2008-10-09 Thread Michael Just
: +86-(0)10-82509086 Mobile: +86-15810805877 Homepage: http://www.yihui.name School of Statistics, Room 1037, Mingde Main Building, Renmin University of China, Beijing, 100872, China On Fri, Oct 10, 2008 at 5:54 AM, Michael Just [EMAIL PROTECTED] wrote: Hello, Is there a way to put multiple

Re: [R] Plot means with error bars - A novice needs help

2008-10-08 Thread Michael Just
assistance you are able to provide me. My limited knowledge in R grows quite rapdily and things become so much more successful for me when I am granted tips and snippets of code form R users like you. Sincerely, Michael Just On Wed, Oct 8, 2008 at 11:25 AM, hadley wickham [EMAIL PROTECTED] wrote

Re: [R] Plot means with error bars - A novice needs help

2008-10-08 Thread Michael Just
= c(0,20), legend = TRUE) On Wed, Oct 8, 2008 at 2:28 PM, Manuel Morales [EMAIL PROTECTED] wrote: On Wed, 2008-10-08 at 12:01 -0500, Michael Just wrote: Thank you all for you suggestions. They are all helpful. However, I have come to a more fundamental problem. Preparing my data to even make

[R] Choose subset for plot use (bwplot)

2008-10-08 Thread Michael Just
How can I choose only two groups from 'Pri_No'? Thank you kindly, Michael Just [[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

Re: [R] Choose subset for plot use (bwplot)

2008-10-08 Thread Michael Just
X6 X7 slope [46] std.slopeCI.SLOPe X.CI.SLOPE 0 rows (or 0-length row.names) Somehow my Pri_No got reclassified to TRUE and FALSE. I did it over and it works. Sorry about the waste of time. But it helped me out. Thanks, Michael Just On Wed, Oct 8, 2008 at 7:37 PM

Re: [R] Multiple hist(ograms) - One plot

2008-10-07 Thread Michael Just
Deepayan, Thanks for the advice. -Micahel On Fri, Oct 3, 2008 at 1:24 AM, Deepayan Sarkar [EMAIL PROTECTED]wrote: On 10/2/08, Michael Just [EMAIL PROTECTED] wrote: Dieter and Thierry: Per you suggestions I have tried: ggplot2 from Thierry: p - ggplot(dat, aes(x=bbContag, y

Re: [R] Multiple hist(ograms) - One plot

2008-10-07 Thread Michael Just
and an aching desire for an answer does not ensure that a reasonable answer can be extracted from a given body of data. ~ John Tukey -- *Van:* Michael Just [mailto:[EMAIL PROTECTED] *Verzonden:* vrijdag 3 oktober 2008 6:52 *Aan:* ONKELINX, Thierry *CC:* r-help@r

Re: [R] Multiple hist(ograms) - One plot

2008-10-02 Thread Michael Just
for an answer does not ensure that a reasonable answer can be extracted from a given body of data. ~ John Tukey -Oorspronkelijk bericht- Van: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Namens Michael Just Verzonden: donderdag 2 oktober 2008 7:29 Aan: r-help@r-project.org Onderwerp: [R

Re: [R] Multiple hist(ograms) - One plot

2008-10-02 Thread Michael Just
? I know I can do: sc_recov.21 - dat[dat$sc_recov=21,] but how could I select all cases where sc_recov = 21, 31 or 41? On Thu, Oct 2, 2008 at 3:33 AM, Michael Just [EMAIL PROTECTED] wrote: Hello, If I use: p - ggplot(dat, aes(x=bbContag, y=..density..)) + geom_histogram() and then: p

[R] Multiple histograms one graphic device groupings

2008-10-01 Thread Michael Just
Hello, I am trying to produce some graphics to visualize my data. I think I want histograms. I have a grouping variable that has 48 different groups. I would like to produce a graphic that contains three of of these groups in the same device (they are time steps). I would like create graphics for

[R] Aggregate issues with subset

2008-10-01 Thread Michael Just
Hello, I used the aggregate function with success on my data frame all. aggregate(all, list(sc_recov), mean) I then made a subset of this data frame: s.all -subset(all, select= c(3,4,5,6,7,8,9,10,11,23)) but when I tried to to run the aggregate function on this subset: aggregate(s.all,

Re: [R] Add title in sink output and possibility of plot-like output for text

2008-09-26 Thread Michael Just
] On Behalf Of Michael Just Sent: Wednesday, September 24, 2008 11:46 AM To: r-help@r-project.org Subject: [R] Add title in sink output and possibility of plot-like output for text Hello, I have been using sink to create text file outputs. sink(summ_model1and2.txt) summary(model1

Re: [R] Add title in sink output and possibility of plot-like output for text

2008-09-26 Thread Michael Just
, something like this should work: sink(summ_model1and2.txt) cat('Summary model1','\n','\n') summary(model1) cat('\n','\n') cat('Summary model2','\n','\n') summary(model2) sink() For the second one, ?savePlot could work. HTH, Jorge On Wed, Sep 24, 2008 at 1:45 PM, Michael

Re: [R] Options for viewing / using results from lm

2008-09-24 Thread Michael Just
. -- *From:* Michael Just [EMAIL PROTECTED] *Sent:* Tuesday, September 23, 2008 2:49 PM *To:* Greg Snow *Cc:* r-help@r-project.org; [EMAIL PROTECTED] *Subject:* Re: [R] Options for viewing / using results from lm Greg, Thanks for the pseudocode and cbind suggestions I think

[R] Add title in sink output and possibility of plot-like output for text

2008-09-24 Thread Michael Just
Hello, I have been using sink to create text file outputs. sink(summ_model1and2.txt) summary(model1) summary(model2) sink() Q1: Is there a way I could add a line of the text above the summary to act like a title? Also, I have been using the following to save plots from the lm function:

[R] Options for viewing / using results from lm

2008-09-23 Thread Michael Just
Hello, I would like to state what I am (trying) to do. I have data set. It has 5749 rows (including the header) and 23 columns. The data contains values related to spatial aspects of the 412 landscapes (over various years). I will be making 2 groups from the data based on spatial extent. I will

Re: [R] Options for viewing / using results from lm

2008-09-23 Thread Michael Just
Hello, Why do you think it is better done in another program? Keeping it in R saves you from the exporting, which you say you are having trouble with. I think it might be better for viewing capabilities. If I had columns of residuals right next to each other I could spot differences amongst

Re: [R] Options for viewing / using results from lm

2008-09-23 Thread Michael Just
Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] project.org] On Behalf Of Michael Just Sent: Tuesday, September 23, 2008 2:06 PM To: r-help@r-project.org Cc: [EMAIL PROTECTED] Subject: Re: [R] Options for viewing / using results from lm Hello, Why do you think

[R] Create groups from data to compute lm?

2008-09-22 Thread Michael Just
Hello, Below are the first two rows from my dataset and the header. This dataset has 5749 rows and I want to select only certain rows to be used based on existing grouping values. I am trying to group the data based on the values under 'ex_bin'. (e.g a group for 250, 251, 252, 500, 501, 502) I

Re: [R] Create groups from data to compute lm?

2008-09-22 Thread Michael Just
Hello again, After getting some help from a reader of this board my issue has been resolved. To group I will now use something like this: sub.all - subset(all, ex_bin==250, select= c(pfor, ED, PD)) before I perform any further analysis, such as lm. Thanks, M Just Hello,