[R] nlme package, fixing variance.covariance matrix of residuals

in nlme? Are there any changes planned to incorporate this feature? Best wishes, I am looking forward to your feedback Anna-Lena W�lwer [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https

[R] Modifying graphs in 'survrec' package

"col" or "lty" has not worked: plot(fit,ylim=c(0,1),xlim=c(0,2000), col=c("red", "blue", "orange"), lty=3) I would be very grateful for any help! Thank you very much in advance! Kind regards, Anna McLean [[alternative HTML version del

Re: [R] seqinr ?: Splitting a factor name into several columns. Dealing with metabarcoding data.

and only used R for the statistics part. I intend to make a permanent swith to R as I can see the benefits. I am just severely frustrated by my own inabilities. Therefore, once again many thank´s for your time! kind regards Anna º`•. . • `•. .• `•. . º`•. . • `•. .• `•. .º

Re: [R] seqinr ?: Splitting a factor name into several columns. Dealing with metabarcoding data.

Thank you! That was a easy and fast solution! May I post a follow-up question? (I am not sure if this would rather should be posted as a new question, but I post it here and then I can re-post it if this is the wrong place to ask this). I am ever so grateful for your help! /Anna

[R] seqinr ?: Splitting a factor name into several columns. Dealing with metabarcoding data.

)), B1298712 = factor(rep(LETTERS[1:2], each = 5))) # Many thank's and with kind regards Anna Zakrisson �`�. . � `�. .� `�. . �`�. . � `�. .� `�. .�`�. . � `�. .� `�. .� Anna Zakrisson Braeunlich PhD student Department of Ecology, Environment and Plant Sciences Stockholm University Svante

[R] working with matrices

column of the matrix and apply the same function, get the single number and place it to (2,3) entry of the matrix and so on how can i do it? thanks anna [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https

[R] lattice, latticeExtra: Adding moving averages to double y plot

in insering this into lattice (I am not an experienced lattice user). I want to keep the data points in the plot. Find below dummy data and the script as well as annotations further describing my question. thank you in advance! Anna Zakrisson mydata- data.frame( Year = 1980:2009, Type

[R] Groupwise sum

OUTPUT? Kindly guide.  Anna [[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/posting-guide.html and provide commented

[R] gRain on R 3.1.0

I am using the gRain package but I can't get it to work under R 3.1.0. It is no longer available in the CRAN. Does anyone have suggestions for how to get a successful installation? A [[alternative HTML version deleted]] __

[R] Please help! Matrices and parameter of time

help and idea!Thank you!Anna S. [[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

[R] ggplot2: removing empty factor levels in facet_grid

()) p Anna Zakrisson Braeunlich PhD student Department of Ecology, Environment and Plant Sciences Stockholm University Svante Arrheniusv. 21A SE-106 91 Stockholm Sweden/Sverige Lives in Berlin. For paper mail: Katzbachstr. 21 D-10965, Berlin - Kreuzberg Germany/Deutschland E-mail: anna.zakris

[R] gam(mgcv) predicting across random effect and understanding narrow CI

Hi everyone, I am new to additive modelling and am surprised by the results of a model I'm working on. I wanted to check with more experienced users to make sure I'm not misunderstanding something basic. *Data:* I have 10 replicated runs from an evolutionary simulation model, measuring the

[R] ggplot2: how to add new names to legend

., Aphanizomenon sp., N. spumigena)) with kind regards Anna Zakrisson Braeunlich PhD student Department of Ecology, Environment and Plant Sciences Stockholm University Svante Arrheniusv. 21A SE-106 91 Stockholm Sweden/Sverige Lives in Berlin. For paper mail: Katzbachstr. 21 D-10965, Berlin

[R] Smooth limits for a surface plot

are, so that the plot does look nice afterwards. Thank you very much for your help! A great new years eve to everybody and a wonderful year 2014. Greetings, Anna [[alternative HTML version deleted]] __ R

[R] Problem with loop, matrix and data frame

! Simeonidou Anna [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] ggplot2: geom_boxplot. Mapping aes factor but with different color scale and hatching

), labels=c(cyanobacteria, seston)) #if I wanted to have seston in italics - how would I do that? with kind regards Anna Zakrisson PhD student Department of Ecology, Environment and Plant Sciences Stockholm University Svante Arrheniusv. 21A SE-106 91 Stockholm Sweden/Sverige Lives in Berlin

[R] Loop over factor returns NA

(not the single parts) What do I do wrong? And how can i fix it? Thank you so much in advance. Kind regards, Anna [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https

[R] Problem with makePSOCKcluster R3.0.1

: In socketConnection(master, port = port, blocking = TRUE, open = a+b, : ip-10-158-31-12:11883 cannot be opened Execution halted Can you help me? Best regards Anna Longari [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https

[R] saving as TIFF - problem with compression

, plotP, add.ylab2 = TRUE) dev.off() this works: tiff(file=myplot8.tiff, bg=white, res=800, width=13.3, height=9.45, units=cm, pointsize=6, compression = lzw) plot(1:100) dev.off() Anna Zakrisson Braeunlich PhD student Department of Ecology, Environment and Plant Sciences Stockholm

[R] ggplot2: two y-axis - any change?

lattice for every second plot. What are your thoughts on this for the future and are there already some possibilities in place? I prefer continusing using ggplot2 because it is brilliant. with kind regards Anna Zakrisson Braeunlich PhD student Department of Ecology, Environment and Plant Sciences

[R] lattice: double y - problem changing axis color after doubleYScale

? such as: # --O-- Nc on the left y-axis (but with lty=1: I could not do a non-dotted line on the keybord). with kind regards Anna Zakrisson Braeunlich PhD student Department of Ecology, Environment and Plant Sciences Stockholm University Svante Arrheniusv. 21A SE-106 91 Stockholm Sweden/Sverige Lives

[R] ggplot2: changing shapes facet_grid - special case.

, angle=00)) + theme(axis.title.x = element_text(size=20, colour=black, family=serif, angle=00))+ theme(axis.title.y = element_text(size=20, colour=black, family=serif, angle=90))+ facet_wrap( ~ year) p Anna Zakrisson Braeunlich PhD student Department of Ecology, Environment and Plant

[R] ggplot2: changing strip text in facet_grid and a legend text problem

(axis.title.y = element_text(size=16, colour=black, family=serif, angle=90)) + scale_fill_discrete(guide = guide_legend(), labels=c(cyanobacteria, zooplankton)) p with kind regards Anna Zakrisson Braeunlich PhD student Department of Ecology, Environment and Plant Sciences Stockholm University

[R] ggplot2: connecting medians of boxes using facet_wrap. Changing median marker.

connected and # and ses medians connected and zpl connected. Not cy-ses-zpl, but that is perhaps obvious. ### Thank you for your time! # with kind regards # A. Zakrisson Anna Zakrisson Braeunlich PhD student Department of Ecology, Environment and Plant Sciences Stockholm University Svante

Re: [R] ggplot2: connecting medians of boxes using facet_wrap. Changing median marker.

)) + theme(axis.text.y = element_text(size = 12, colour=black, family=serif, angle=00)) + geom_hline(yintercept=0, linetype=3) #draws dotted line at 0 p # Anna Zakrisson Braeunlich PhD student Department

Re: [R] lattice: bwplot: getting two plots into one.

, col=black), x = .35, y = 1.06, corner = c(2, 2)), #legend position ylab = expression(paste(,delta^{15}, N)), scales=list(x=list(at=(3*(1:12)-1), labels=c(A, B, C Anna Zakrisson Braeunlich PhD student Department of Ecology, Environment

[R] ANOVA nested model with random and fixed effects

sample - one sediment type. What is the correct way make the model: My guess is something like: aov(Count ~ Sediment+Error(Lake/Zone/Sample)) Sincerely, Anna [[alternative HTML version deleted]] __ R

[R] lattice: bwplot: getting two plots into one.

, ..., group.number) { panel.bwplot(x + (group.number-1.8)/3, y, ...) }) trellis.focus(panel, 1, 1) panel.abline(h=0, lty=2) #lwd trellis.unfocus() trellis.focus(panel, 2, 1) panel.abline(h=0, lty=2) #lwd trellis.unfocus() Anna Zakrisson Braeunlich PhD student Department of Ecology

Re: [R] lattice: bwplot: getting two plots into one.

other - again, this will look terribel on this data, but will be OK for mine. Is this possible in lattice? Maybe using different data frames and make different plots and superpose them on top of each other? Ideas? ggplot2? with kind regards Anna Zakrisson Braeunlich PhD student Department

Re: [R] Lattice: bwplot - changing box colors in legend and plot when using panel.groups = function... and panel = panel.superpose

)), superpose.polygon=list(col=c(gray,white)), # legend par.sub=font.settings) as adviced by you. With kind regards Anna Zakrisson Braeunlich PhD student Department of Ecology, Environment and Plant Sciences Stockholm University Svante Arrheniusv. 21A SE-106 91 Stockholm Sweden/Sverige Lives

[R] Lattice: bwplot - changing box colors in legend and plot when using panel.groups = function... and panel = panel.superpose

, ...) }) Anna Zakrisson Braeunlich PhD student Department of Ecology, Environment and Plant Sciences Stockholm University Svante Arrheniusv. 21A SE-106 91 Stockholm Sweden/Sverige Lives in Berlin. For paper mail: Katzbachstr. 21 D-10965, Berlin - Kreuzberg Germany/Deutschland E-mail: anna.zakris...@su.se

Re: [R] Ggplot2: Moving legend, change fill and removal of space between plots when using grid.arrange() possible use of facet_grid?

. once again, thank you for your time. Anna Summ - ddply(mydata, .(factor3,factor1), summarize, mean = mean(var1, na.rm = FALSE), sdv = sd(var1, na.rm = FALSE), se = 1.96*(sd(var1, na.rm=FALSE)/sqrt(length(var1 Summ$Grouping - c(AB, AB, CD

Re: [R] Ggplot2: Moving legend, change fill and removal of space between plots when using grid.arrange() possible use of facet_grid?

Hi, have managed to get rid of the facet labels (so do not spend your time explaining that to me). Tthere were some old code out there which did not work. My only remaining issue is how to add the axis labels to the plot without labels. Anna Summ - ddply(mydata, .(factor3,factor1

[R] ggpliot2: reordering of factors in facets facet.grid(). Reordering of factor on x-axis no problem.

automatically from my dataframe using facet.grid() previously. I am simply afraid of displaying the wrong things without noticing it. Here is my code and comments: I want to have my factor 1 in the order: F, E, C, D, A, B instead of alphabetical. My normal methods do not work. with kind regards Anna

Re: [R] Ggplot2: Moving legend, change fill and removal of space between plots when using grid.arrange() possible use of facet_grid?

GREAT! Thank you! Will try this! Anna Anna Zakrisson Braeunlich PhD student Department of Ecology Environment and Plant Sciences Stockholm University Svante Arrheniusv. 21A SE-106 91 Stockholm Sweden Lives in Berlin. For paper mail: Katzbachstr. 21 D-10965, Berlin - Kreuzberg Germany

Re: [R] ggpliot2: reordering of factors in facets facet.grid(). Reordering of factor on x-axis no problem.

, thank you so much! Another day closer to disputation... :-S Anna Anna Zakrisson Braeunlich PhD student Department of Ecology Environment and Plant Sciences Stockholm University Svante Arrheniusv. 21A SE-106 91 Stockholm Sweden Lives in Berlin. For paper mail: Katzbachstr. 21 D-10965, Berlin

[R] Ggplot2: Moving legend, change fill and removal of space between plots when using grid.arrange() possible use of facet_grid?

(p1, p2, p3, ncol=2) Anna Zakrisson Braeunlich PhD student Department of Ecology Environment and Plant Sciences Stockholm University Svante Arrheniusv. 21A SE-106 91 Stockholm Sweden Lives in Berlin. For paper mail: Katzbachstr. 21 D-10965, Berlin - Kreuzberg Germany/Deutschland E-mail

Re: [R] Ggplot2: Moving legend, change fill and removal of space between plots when using grid.arrange() possible use of facet_grid?

this problem. I also found: theme(legend.justification=c(1,0), legend.position=c(1,0)) for the legend justification. Now this starts to look like how I want it to! Thank you so much for your time! Best regards Anna Zakrisson Braeunlich PhD student Department of Ecology Environment and Plant Sciences

Re: [R] Ggplot2: Moving legend, change fill and removal of space between plots when using grid.arrange() possible use of facet_grid?

Dear Stephen, thank you for the tip regarding Incskape. I had never heard of it before. It looks extremely useful! With kind regards Anna Zakrisson Braeunlich PhD student Department of Ecology Environment and Plant Sciences Stockholm University Svante Arrheniusv. 21A SE-106 91 Stockholm Sweden

Re: [R] Ggplot2: Moving legend, change fill and removal of space between plots when using grid.arrange() possible use of facet_grid?

tweaking grid.arrange. This will come in handy! With kind regards Anna Zakrisson Braeunlich PhD student Department of Ecology Environment and Plant Sciences Stockholm University Svante Arrheniusv. 21A SE-106 91 Stockholm Sweden Lives in Berlin. For paper mail: Katzbachstr. 21 D-10965, Berlin

[R] How to plot 2 continous variables on double y-axis with 2 factors: ggplot2, gplot, lattice, sciplot?

this in ggplot2? I have also tried plotmeans() in the sciplot package, but was unsuccessful. Sincerely Anna Zakrisson library(plotrix) ?brkdn.plot par(family=serif,font=1) brkdn.plot(y1, factor1,factor2, data=mydata, mct=mean,md=sd, main=, cex=0.8, stagger

[R] introducing jitter in overlapping graphs using ggplots (plotmeans). Also sciplot.

this and the frustration slowly approaches anger (in myself) :-) Yours sincerely Anna Zakrisson Braeunlich PhD Student Department of Systems Ecology Stockholm University Svante Arrheniusv. 21A SE-106 91 Stockholm Lives in Berlin. For fast paper-mail responses, please use the below address: Katzbachstr. 21 D-10965

[R] cptable too long for source -- gRain problem

I have a code using the gRain package and there cptable is used to get the conditional probability tables. My problem is that when I want to use this table I get an error message that the file is too long to keep source. I have tried to change the option keep.source but without luck. Anyone

[R] worldmap_region/country problem

to not display borderlines between some countries but to show the borders between others? Thank you very much for your help! Anna Fechner PricewaterhouseCoopers Aktiengesellschaft Wirtschaftsprüfungsgesellschaft Vorsitzender des Aufsichtsrates WP StB Dr. Norbert Vogelpoth Vorstandsmitglieder WP

[R] R

Hi! How Can I create a directory in R to save my files? [[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] ensembleBMA pit function warnings

Hello Chris,  I 've found two other issues  with MAE and CRPS, giving warning  when running examples. I've the same issue on my data.  Hope that you could find some time to take a look here. Thank you Anna library(ensembleBMA) Loading required package: chron example(MAE) MAE   data

Re: [R] plotting principal components on geographic map

Hi Marco, a small example will be helpful to get better to the point. But I suggest you either to address this question to  R geo mailing list  https://stat.ethz.ch/mailman/listinfo/r-sig-geo and eventually take a look to GeoXp library. Cheers Anna   Anna Freni Sterrantino Department

Re: [R] GAM without intercept

it helps Anna Anna Freni Sterrantino Department of Statistics University of Bologna, Italy via Belle Arti 41, 40124 BO. Da: SAEC sergio.es...@uach.cl A: r-help@r-project.org Inviato: Giovedì 11 Ottobre 2012 0:22 Oggetto: [R] GAM without intercept Hi everybody, I

[R] interaction testing in GAM

Hello! I'd like to know if it is correct to test with anova two models specified like this: m1=y~x1+s(x2,by=x3),family=poisson m0=y~x1+s(x2),family=poisson anova(m1,m0) Cheers Anna   Anna Freni Sterrantino Department of Statistics University of Bologna, Italy via Belle Arti 41, 40124 BO

[R] I: ensembleBMA pit function warnings

? Anna Freni Sterrantino Ph.D Student Department of Statistics University of Bologna, Italy via Belle Arti 41, 40124 BO. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] inverse binomial in R

(function(x) pbinom(k, n, x) - p, c(0, 1)) :    f() values at end points not of opposite sign while stata gen p3=invbinomial(50,50, 0.4) . display p3 0 . gen p4=invbinomial(50,50, 0.6) . display p4 0 Thanks Cheers  Anna Da: Duncan Murdoch murdoch.dun

Re: [R] inverse binomial in R[solved]

Hello Peter, Duncan, Dan and RList! How I suspected... I managed to do the  homework and forced invbinomal to return zero ( as in Stata) once that k==n. Thank you for all you replies! Cheers Anna Da: peter dalgaard pda...@gmail.com Cc: Duncan Murdoch

[R] inverse binomial in R

=guides:tutorials:regression:table I tried to replicate using qbinom the results obtained in  invbinomial(10,5, 0.5) .54830584 but with no success. Thank you Cheers Anna Anna Freni Sterrantino Department of Statistics University of Bologna, Italy via Belle Arti 41, 40124 BO

[R] 3D Black-Scholes Graph Help!

Hello all! I would like to create a 3d plot, with the option price explained by the underlying price and time. Unfortunately, I can't quite get it to work. I would very much appreciate your help! Thanks, Anna # Black-Scholes Option Graph library(lattice) blackscholes - function(s, k

[R] help

) to run this with source(name_of_file), but it doesn't work. Shouldn't a plot come up automatically when I run it? What am I doing wrong? It knows what x and y is, but I don't get an error of what might be wrong. source(name_of_file) x [1] 10.4 5.6 3.1 6.4 21.7 Best, Anna

Re: [R] help

of '1G'. Am I not allowed to use numerical values? Best, Anna On Tue, 10 Jan 2012 23:02:04 +0100, Anna Olofsson anol2...@student.su.se wrote: Hi, I'm pretty new at programming and with the R language. I'm just trying to get familiar with R and wrote a script in gedit (should I use emacs instead

[R] Autocorrelation values? How to extract?

is to multiply the dependent variable autocorrelation with the independent variable autocorrelation and then multiply by (N-j)/N where N is the sample size and j is the lag...calculate z-value...adjust my p-value...Sincerely Anna Zakrisson Braeunlich PhD Student Department of Systems Ecology Stockholm

[R] corCompSymm in gamm()?

! Anna Zakrisson Braeunlich PhD Student Department of Systems Ecology Stockholm University Svante Arrheniusv. 21A SE-106 91 Stockholm E-mail: a...@ecology.su.se Tel work: +46 (0)8 161103 Mobile: +46-(0)700-525015 Web site: http://www.ecology.su.se/staff/personal.asp?id=163 º`•. . • `â

[R] Labels in xy-plots

should be carried out on the data defined by the groups option (here: year) When I call this, I cannot separate the lines from my different years. I would like them to have different colours/other markings (irrelevant for me) and identified on the side of the graph. Thank's! Anna Zakrisson

Re: [R] Contour on top of 2d histogram

Hi, thanks for the suggestion! I had tried it before, but it did not work - this was probably because I was using the image function to plot the 2d histogram. When I use directly hist2d and then contour with add=T it works. Thanks again Anna Von: R

[R] Contour on top of 2d histogram

? Thanks, Anna [[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/posting-guide.html and provide commented, minimal, self

[R] converting commas to points

, Anna -- Der Inhalt dieser E-Mail ist vertraulich. Sollte Ihnen die E-Mail irrtümlich zugesandt worden sein, bitte ich Sie, mich unverzüglich zu benachrichtigen und die E-Mail zu löschen. This e-mail is confidential. If you have received it in error, please notify me immediately and delete

Re: [R] converting commas to points

Sorry guys, I allready found the solution. Excell showed some of the numbers in the format: 1,90053-E05 and R interpreted it as 1.9... I changed that in Excel Cheers, Anna Am 6. Oktober 2011 17:48 schrieb Uwe Ligges lig...@statistik.tu-dortmund.de: On 06.10.2011 17:39, Anna Lee wrote: Hello

Re: [R] Deleting Rows based on Factor and Time Period

, that isoWeekYear is great for week and year but not so great for month. If you have any ideas, please let me know! Thanks! Anna On Tue, Sep 13, 2011 at 6:55 PM, Mikkel Grum mi2kelg...@yahoo.com wrote: The following will get you the first stock in each week. Is that useful? install.packages(surveillance

Re: [R] Deleting Rows based on Factor and Time Period

Alright - for the first stock in the month, I pretty much copied Mikkel's code and got what I wanted. Thanks! Anna install.packages(surveillance) library(surveillance) alldat$year - isoWeekYear(alldat$mydate)$ISOYear alldat$months-months(alldat$mydate) alldat - alldat[order(alldat$year, alldat

[R] Deleting Rows based on Factor and Time Period

the index that would show me which rows in divall/alldat need to be deleted. I hope this has been clear. Please let me know if you need any more information! Thank you very much! Anna [[alternative HTML version deleted]] __ R-help@r-project.org

[R] If NA Problem!

project. The thing I most want to understand is how, after specifying a certain condition, one may save certain data that occurs when that condition is met. I hope I have been clear enough! Thank you very much for your help! Anna biglist-list(a=1:4,b=2:6) lilwin-list(x=NA,y=2) lilloss-list(m=1,n

[R] defining id argument in geeglm

(month,site)* or *id=interaction(month,type*)) Thanks for your help... Anna [[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] defining id argument in geeglm

.) On Mon, 29 Aug 2011, Anna Mill wrote: Hi all, I am trying to do a generalized estimating equation (GEE) with the geepack package and I am not 100% sure what exactly the id argument means. It seems to be an important argument because results differ considerably defining different clusters. I

[R] Testing Specific Hypothesis

to do. I have researched for a couple of hours regarding this problem but to no avail! Thanks! Anna otestme values ind 10.001008012 AAPL.UW.Equity 20.015518087 AAPL.UW.Equity 30.013221459 AAPL.UW.Equity 40.012195734 AAPL.UW.Equity 5 -0.026750298

[R] Fwd: Interpreting parameters of sigmoid fct

: modPoplar- nls(Diameter ~ d*(1-exp(-b *Age))^a ,start=list(a=20,b=0.9,d=33)) I attached the graph, too. Hoping for your answers! Best, Anna -- Telefon: 01577-7805329 E-Mail: gretschel.a...@googlemail.com

[R] what characteristics of model curve do parameters denote

: modPoplar- nls(Diameter ~ d*(1-exp(-b *Age))^a ,start=list(a=20,b=0.9,d=33)) I attached the graph, too. Hoping for your answers! Best, Anna -- Der Inhalt dieser E-Mail ist vertraulich. Sollte Ihnen die E-Mail irrtümlich zugesandt worden sein, bitte ich Sie, mich unverzüglich zu benachrichtigen

[R] what characteristics of model curve do parameters denote

: modPoplar- nls(Diameter ~ d*(1-exp(-b *Age))^a ,start=list(a=20,b=0.9,d=33)) I attached the graph, too. Hoping for your answers! Best, Anna -- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Telefon: 01577-7805329 E-Mail: gretschel.a

[R] Loop trouble with Excel Serial Numbers!

)){ mynewdate[i]-as.Date(thedate[i,1],origin=1899-12-30) } print(mynewdate) } exdate() Thank you very much! Anna [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

[R] Interpreting parameters of sigmoid fct

: modPoplar- nls(Diameter ~ d*(1-exp(-b *Age))^a ,start=list(a=20,b=0.9,d=33)) I attached the graph, too. Hoping for your answers! Best, Anna -- Der Inhalt dieser E-Mail ist vertraulich. Sollte Ihnen die E-Mail irrtümlich zugesandt worden sein, bitte ich Sie, mich unverzüglich zu benachrichtigen

Re: [R] instal tar.gz package on windows

Hi Duncan, I have tried to install a tar.gz package following your instructions (https://stat.ethz.ch/pipermail/r-help/2008-August/169599.html) but without success. Here are the steps I followed: I installed the last version of Rtools and ran Rcmd INSTALL rJava_0.8-8.tar.gz and got the error

[R] vglm: warnings and errors

literature that would help me understand my problems. Thanks Anna __ 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

[R] Add permanently environment variable

Hello everyone, I know how to add a folder path to my EV path but it only works for the current R session. Is there a way to add it permanently? Here is my code: Sys.setenv(PATH=paste(C:\\Program Files\\Java\\jre1.6.0_13\\bin;, Sys.getenv(x=PATH), sep=)) Thanks a lot! -- View this message in

[R] Errorbars

-as.vector(tapply(Sporangia,list(Host,Isolate),sd)) c-sqrt(b/36) errbar(c(1.5,2.5,3.5,5.5,6.5,7.5),(tapply(Sporangia,list(Host,Isolate),mean)),a+c,a-c,add=T) Cheers Anna [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https

[R] Error bars

,grey89),beside=T) col=c(grey39,grey64,grey89) legend(topright,inset=.05,title=Host,c(European Larch,Hybrid Larch,Japanese Larch),fill=col,horiz=FALSE) Any help would be greatly appreciated, Anna [[alternative HTML version deleted]] __ R

Re: [R] glm with binomial errors - problem with overdispersion

errors. On Mon, 13 Jun 2011, Anna Mill wrote: Dear all, I am new to R and my question may be trivial to you... I am doing a GLM with binomial errors to compare proportions of species in different categories of seed sizes (4 categories) between 2 sites. You have types and fragments

Re: [R] glm with binomial errors - problem with overdispersion

), and not independent binomial data. sorry for the dumb question: so do you think, that my data is independent and the model appropriate? Thanks, Anna -- Peter Dalgaard Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd

[R] glm with binomial errors - problem with overdispersion

models attached)? Thanks a lot for your help! Anna first model with binomial error structure: success-c(14,43,44,1,13,28,56,8) failure-c(88,59,58,101,92,77,49,97) fragment-c(1,1,1,1,2,2,2,2) type-c(1,2,3,4,1,2,3,4) y-cbind(success,failure) model-glm(y~fragment*type,binomial) summary(model

[R] biplots for PCA

the label names and if not what the best alternative is. Thanks Anna pca-biodata[,3:10] model-prcomp(pca,scale=TRUE) summary(model) biplot(model) -- Dr Anna R. Renwick Research Ecologist British Trust for Ornithology, The Nunnery, Thetford, Norfolk, IP24 2PU, UK Tel: +44 (0)1842 750050; Fax: +44 (0

[R] Comparison of two penalized spline fits in mixed model framework

of Observations: 1549 Number of Groups: 2 Warning message: In pt(q, df, lower.tail, log.p) : NaNs produced Kindly Regards, Anna-Leena Orsama __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

[R] Snow/Snowfall hangs on windows 7

, type=SOCK). When running process explorer I see the R-program creating a TCP connection in the state LISTENER but no ESTABLISHED happens. Regards Anna [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https

[R] statistical question

Dear List! I want to compare medians of non normal distributed data. Is it possible and usefull to calculate 95% confidence intervals for medians? And if so - how can this be achieved in R? Thanks a lot! Anna -- Der Inhalt dieser E-Mail ist vertraulich. Sollte Ihnen die E-Mail irrtümlich

[R] another statistical question

data and therefore comparison of means, or non normal distributed data and therefore comparison of medians? Cheers, Anna -- Der Inhalt dieser E-Mail ist vertraulich. Sollte Ihnen die E-Mail irrtümlich zugesandt worden sein, bitte ich Sie, mich unverzüglich zu benachrichtigen und die E-Mail zu

[R] generate random numbers

Hey List, does anyone know how I can generate a vector of random numbers from a given distribution? Something like rnorm just for non normal distributions??? Thanks a lot! Anna -- Der Inhalt dieser E-Mail ist vertraulich. Sollte Ihnen die E-Mail irrtümlich zugesandt worden sein, bitte ich

Re: [R] predicting values from multiple regression

Dennis: Thank you so much! I got it now - it just works perfectly. Thanks a lot! Anna 2011/3/21 Dennis Murphy djmu...@gmail.com: Hi: To amplify Ista's and David's comments: (1) You should not be inputting separate vectors into lm(), especially if you intend to do prediction. They should

Re: [R] predicting values from multiple regression

Dennis: thank you so much! I got it now and it works just perfectly. Thanks a lot to the others too! Anna 2011/3/21 Dennis Murphy djmu...@gmail.com: Hi: To amplify Ista's and David's comments: (1) You should not be inputting separate vectors into lm(), especially if you intend to do

[R] predicting values from multiple regression

)) The predicted values are however strange. Now I do not know weather just the model does not fit the data (actually all coefficiets are significant and the plot(model) shows a good shape) or wether I did something wrong with my prediction command. Does anyone have an idea??? -- Thanks a lot, Anna

Re: [R] predicting values from multiple regression

thought with the predict function the programm would calculate a value from the model function for every value of calP... ? 2011/3/20 Ista Zahn iz...@psych.rochester.edu: Hi Anna, On Sun, Mar 20, 2011 at 2:54 PM, Anna Lee ana-...@web.de wrote: Hey List, I did a multiple regression and my final

[R] R² for non-linear model

Dear List, how can I obtain the value of r suqared for a non-linear model? For linear models it can be found in the summary() of the model but for non-linear models I just don't know. Please help! Anna __ R-help@r-project.org mailing list https

Re: [R] R² for non-linear model

Am 16.03.2011 18:15, schrieb David Winsemius: On Mar 16, 2011, at 12:53 PM, Anna Gretschel wrote: Dear List, how can I obtain the value of r suqared for a non-linear model? For linear models it can be found in the summary() of the model but for non-linear models I just don't know. Please

Re: [R] R² for non-linear model

Am 16.03.2011 18:19, schrieb Joshua Wiley: Dear Anna, What is your goal in obtaining a value for R^2 ? I believe it is not provided for a non-linear model, because it does not make much sense. It certainly will not have the same interpretation as in a linear model, and all the ways it *could

Re: [R] R² for non-linear model + comparing linear + non-linear models

Dear Bert, so what can I do to obtain a goodness of fit for a non-linear model if r² does not work? And here comes my next question: is it apropriate to comopare a linear and a non-linear model with anova()? Thank you so much for answering, Anna Am 16.03.2011 18:54, schrieb Bert Gunter

Re: [R] R² for non-linear model + comparing linear + non-linear models

Am 16.03.2011 19:29, schrieb Heiman, Thomas J.: Hi Anna, AIC and BIC are good criteria for determining degree of model fit.. Sincerely, tom Thomas Heiman, PhD Info Systems Eng, Sr The MITRE Corporation | Center for Enterprise Modernization Office: 703-983-2951 | thei...@mitre.org

[R] predicting values from model function

Dear List, I have fitted a spherical function to my variogram using variofit(...) from GeoR. Now I would like to predict some data with the function predict(object,...) from package stats. Does anyone know wether this works and if it does how to do it? Thanks a lot! Anna

[R] testing power of correlation

Dear List, does anyone know how I can test the strength of a correlation? Cheers, Anna __ 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

  1   2   3   >