text_combined: carrascotorrecilla roman r romancarrasco06imperialacuk dear trying fit non linear model dispersal data seems sometimes fit model data good start getting singular gradient errors however modify slightly data wont occurr also tried changing initial parameter values algorithm fitting nls didnt help ended programming software fitting us using searching algorithm map parameters minimising sum squares get fits r models could fit data well still dont know r could fit models wonder calculate standard errors pvalues non linear regression parameters trying find information find linear regression could tell please calculated nls r also wrong nls able fit model many thanks advance roman carrasco original message rhelpbouncesrprojectorg mailtorhelpbouncesrprojectorg behalf rhelprequestrprojectorg sent 04 august 2008 1100 rhelprprojectorg subject rhelp digest vol 66 issue 4 send rhelp mailing list submissions rhelprprojectorg subscribe unsubscribe via world wide web visit httpsstatethzchmailmanlistinforhelp via email send message subject body help rhelprequestrprojectorg reach person managing list rhelpownerrprojectorg replying please edit subject line specific contents rhelp digest todays topics 1 r command history like matlabs peter dalgaard 2 output components gam gavin simpson 3 gaps time series gabor grothendieck 4 drop1 seems give unexpected results compare anova thomas p c chu 5 get inset graph ie graph within graph mark difford 6 get inset graph ie graph within graph gabor grothendieck 7 losing rownames matrix operations gabor grothendieck 8 neighbors element matrix rostam shahname 9 computer run 64 bit ubuntur tom la bone 10 installing r apache web server ajay ohri 11 access object inside function jonas 12 access object inside function duncan murdoch 13 running strucchange burgs 14 gmaps projection inset john p burkett 15 access object inside function patrick burns 16 access object inside function jonas 17 sans objet pir2jv 18 drop1 seems give unexpected results compare anova prof brian ripley 19 sans objet dirk eddelbuettel 20 bubble plots franz mueter 21 computer run 64 bit ubuntur prof brian ripley 22 computer run 64 bit ubuntur dirk eddelbuettel 23 access object inside function duncan murdoch 24 missing f statistic anovagam michael milligan 25 computer run 64 bit ubuntur prof brian ripley 26 convert date decimal days cls59 27 convert loop apply charles c berry 28 convert date decimal days prof brian ripley 29 convert date decimal days cls59 30 determining model parameters rkevinburtoncharternet 31 determining model parameters ben bolker 32 changing values andrew ramsey 33 changing values jim holtman 34 bubble plots michael bibo 35 determining model parameters ben bolker 36 changing values j dougherty 37 sweave ggplot2 sornnorngdpivicgovau 38 decomposing tests interaction terms mixedeffects models andrew robinson 39 95ci around median fernando marmolejo ramos 40 95ci around median simon blomberg 41 graph morneau franois 42 major difference multivariate analyses spss r draga r 43 sweave ggplot2 onkelinx thierry 44 decomposing tests interaction terms mixedeffects models peter dalgaard 45 major difference multivariate analyses spss r peter dalgaard 46 95ci around median gavin simpson 47 unexpected nls behaviour solved keith jewell 48 format output file way want pierre8rou 49 source script file straight subversion repository martin maechler message 1 date sun 03 aug 2008 122648 0200 peter dalgaard subject r r command history like matlabs gad abraham cc rhelprprojectorg prof brian ripley messageid 489587e870108biostatkudk contenttype textplain charsetiso88591 formatflowed quite true cant type anything happens still reverseisearch get effect r l c reverseisearchl ls reverseisearchl ls reverseisearchl ls reverseisearchld levelsftpain3 listnonenoneintermediatecmildmediumseveresevere snaps press c twice try typing ctrlc reverse search mode nothing displays console beep sounds doesnt matter many times press ctrlc happens think depends history least typing something ls previously way press either enter shows previously highlighted command ls runs use key ctrlc actually true must hit 2nd c im using ubuntu hardy fedora 9 im pretty sure generic systems use libreadline o__ peter dalgaard ster farimagsgade 5 entrb c _ dept biostatistics po box 2099 1014 cph k university copenhagen denmark ph 45 35327918 pdalgaardbiostatkudk fax 45 35327907 message 2 date sun 03 aug 2008 124550 0100 gavin simpson subject r output components gam andreasgwinteryahoocom cc rhelprprojectorg messageid 121776395026506camelgraptoleberisgeoguclacuk contenttype textplain gamgam mgcvgam latter understand want gamout contains gam model fitted using gam mgcv outterm predictgamout type terms gives contributions covariate response fitted model read predictgam details type terms actually returns note wary accessing components model fits directly instead get various bits using accessor functions coefs coefficients fitted fitted values etc often components dont hold exactly think processed correctly accessed appropriate accessor function hth g sat 20080802 1847 0700 andreas winter wrote would like request help following trying use generalized additive model gam examine density distribution fish function latitude longitude continuous variables year categorical variable model written gamout gamdensity slat slon asfactoryear fitted model prediction link function gamout linearpredictors presumably gamoutlinearpredictors must derived combination original predictor variables lat lon year corresponding coefficients intercept gamoutcoefficients smooth outputs gamoutsmooth andor gamoutsp comparison glm model glmout glmdensity lat lon asfactoryear simply glmoutlinearpredictors glmoutcoefficientsintercept glmout coefficients year glmoutcoefficientslat x lat glmout coefficientslon x lon problem figure get equivalent gam model would like know decompose gamout linearpredictors components evaluate effects different predictor variables separately would appreciate comments help thank andreas winter blacksburg va alternative html version deleted ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code message 3 date sun 3 aug 2008 090224 0400 gabor grothendieck subject r gaps time series rkevinburtoncharternet cc rhelprprojectorg messageid 971536df0808030602x6abdb431gaccd0ae3a3251b66mailgmailcom contenttype textplain charsetiso88591 using zoo libraryzoo setseed1 x zoornorm7 asnumeric17 zoornorm4 c1 2 6 7 xy mergexy xyy xyx want difference times existing x sat aug 2 2008 516 pm wrote like fact subtracting two time series objects effort align series time series begins 1 one begins 2 subtraction operation makes sure proper values subtracted unclear best way build time series holes say data day 1267 one time series 1234567 another construct time series indicate missing data 345 first series eventually would like find difference two series difference series 2 minues series one indexes 345 resultant difference would series 2 value resultant series would length 7 thank kevin ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code message 4 date sun 03 aug 2008 095359 0400 thomas p c chu subject r drop1 seems give unexpected results compare anova rhelprprojectorg messageid 8cac38783a1f557dcc1db6webmailma04sysopsaolcom contenttype textplain charsetusascii formatflowed thanks mr dalgaard advice everyone else contributed inclusion error term end simsety line cure problems drop1 step suppose inexperience carrying simulations caused gaffe thomas ________________________________________________________________________ aol email goes mobile read aol emails whilst move sign free aol email account unlimited storage today message 5 date sun 3 aug 2008 065541 0700 pdt mark difford subject r get inset graph ie graph within graph rhelprprojectorg messageid 18798809posttalknabblecom contenttype textplain charsetusascii hi arthur done quite easily using appropriate arguments listed par approaches readymade functions exist several packages tend use addscatter package ade4 short function easy customize works well straight box hth mark arthur roberts wrote hi get inset graph ie graph within graph input greatly appreciated best wishes art university washington department medicinal chemistry ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code view message context httpwwwnabblecomhowdoigetaninsetgraph28iegraphwithina graph29tp18796563p18798809html sent r help mailing list archive nabblecom message 6 date sun 3 aug 2008 100354 0400 gabor grothendieck subject r get inset graph ie graph within graph arthur roberts cc rhelpstatmathethzch messageid 971536df0808030703u4c74d306ub84cc6a741a670b0mailgmailcom contenttype textplain charsetiso88591 sun aug 3 2008 351 arthur roberts wrote get inset graph ie graph within graph input greatly appreciated see subplot teachingdemos package message 7 date sun 3 aug 2008 100724 0400 gabor grothendieck subject r losing rownames matrix operations rcoder cc rhelprprojectorg messageid 971536df0808030707y3b83a302l4f56bf653971441dmailgmailcom contenttype textplain charsetiso88591 constructing zoo object well covered zoo three vignettes vignettezoo vignettezooquickref vignettezoofaq sat aug 2 2008 1027 rcoder wrote thanks reply gabor already dates column dataframe column rownames possible preserve whilst still making data set suitable rollapply thanks rcoder gabor grothendieck wrote regular convert ts zoo irregular convert zoo reason expect rollapply work objects classes read ts zoo ts note start frequency arguments sat aug 2 2008 750 rcoder wrote hi everyone data frame following format matdate rownames id1 id2 id3 1 date1 2 date1 3 date3 etc perform rollapply statement matrix without converting matrix time series ie mattstsmatdate use start frequency arguments see ts rollapply statement work matminrollapplymatts 2by2 min narmf apply rollapply statement dataframe get following error error could find function rollapply problem convert dataframe matrix time series matrix lose dates rownames column want know anyone could suggest way get around problem ie keep rownames column place use rollapply statement thanks rcoder view message context httpwwwnabblecomlosingrownamesinmatrixoperationstp18788509p1 8788509html sent r help mailing list archive nabblecom ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code view message context httpwwwnabblecomlosingrownamesinmatrixoperationstp18788509p1 8789688html sent r help mailing list archive nabblecom ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code message 8 date sun 3 aug 2008 120358 0300 rostam shahname subject r neighbors element matrix rhelprprojectorg messageid contenttype textplain hi r users wonder function would render value neighbors given element ij matrix thanks rostam alternative html version deleted message 9 date sun 3 aug 2008 080404 0700 pdt tom la bone subject r computer run 64 bit ubuntur rhelprprojectorg messageid 18799376posttalknabblecom contenttype textplain charsetusascii reading 64bit systems software still somewhat uncertain things dell dimension xps 400 dual core intel pentium 940 32 ghz 4 gb memory currently dual boot system xp professional ubuntu 804 32 bit simply install ubuntu 804 64 bit system run 64 bit linux version r would appreciate advice might willing offer tom view message context httpwwwnabblecomwillthiscomputerrun64bitubunturtp18799376 p18799376html sent r help mailing list archive nabblecom message 10 date sun 3 aug 2008 204519 0530 ajay ohri subject r installing r apache web server rhelprprojectorg messageid 4bc14e460808030815y2b6767f0j7d1b2c383e7e6c85mailgmailcom contenttype textplain hi list run r windows terminal server works fine need install apache web server amazon ec 2 server test r gui computing cloud using framework wrote aboutintent benchmark r superior algorithms big files 12 gb mainstream software also data compression encryption technique package r intent reduce bandwidth transfer time costs regards ajay decisionstats com alternative html version deleted message 11 date sun 3 aug 2008 172948 0200 jonas subject r access object inside function rhelprprojectorg messageid 1875de320808030829u70eadd9ao6b2a3908f78c2930mailgmailcom contenttype textplain charsetiso88591 hi apologies covered numerous times lack search andor r skills stopping finding solution im trying access object defined created inside function object returned function cant seem get work keep getting object found error message thought solution found changing frameenvironment cant seem understand correctly want something like myfunction functionfrompackagexy plot objectinsidepackagefunction im using r 271 linux sincerely jonas message 12 date sun 03 aug 2008 113432 0400 duncan murdoch subject r access object inside function jonas cc rhelprprojectorg messageid 4895d00860405statsuwoca contenttype textplain charsetiso88591 formatflowed 03082008 1129 jonas wrote hi apologies covered numerous times lack search andor r skills stopping finding solution im trying access object defined created inside function object returned function cant seem get work keep getting object found error message thought solution found changing frameenvironment cant seem understand correctly want something like myfunction functionfrompackagexy plot objectinsidepackagefunction pseudocode doesnt really make sense didnt save result functionfrompackage anywhere change myfunction z functionfrompackagexy plot z perfectly standard question expecting r look find objectinsidepackagefunction duncan murdoch im using r 271 linux sincerely jonas ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code message 13 date sun 3 aug 2008 084622 0700 pdt burgs subject r running strucchange rhelprprojectorg messageid 18799748posttalknabblecom contenttype textplain charsetusascii thanks replies appreciated prof brian ripley wrote sat 2 aug 2008 burgs wrote yes im using windows im getting error use source command seems like path r file im putting although correct paththanks response anyway see rwfaq q216 want send code r console run paste provided complete line newline end run otherwise says need hit returnenter run macqueen wrote time use source command sourcepath_to_file maybe press return key pasting console assume youre windows dont use basically idea 245 pm 0700 8208 burgs wrote thank post think question simply run code contained r text file manual see highlightedselected run using ctrl character however simply seems paste code rconsole without actually anything macqueen wrote would suggest imitating example help page function youre trying use since havent shown example tried including error messages theres really way anyone help beyond generalities like mine please review posting guide whose url bottom every email rhelp 211 pm 0700 8208 burgs wrote greetings im complety new r question ive read couple manuals im problem getting something run properly id like attempt use strucchange package sample data however im trouble understanding proper syntax commands basically want take data run strucchange able observe plots suggestions greatly appreciated view message context http www nabblecomrunningstrucchangetp18791489p18791489html sent r help mailing list archive nabblecom ______________________________________________ rhelprprojectorg mailing list https statethzchmailmanlistinforhelp please read posting guide http www rprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code macqueen lawrence livermore national laboratory livermore ca usa 9254231062 macqllnlgov ______________________________________________ rhelprprojectorg mailing list https statethzchmailmanlistinforhelp please read posting guide http www rprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code view message context http www nabblecomrunningstrucchangetp18791489p18793670html sent r help mailing list archive nabblecom ______________________________________________ rhelprprojectorg mailing list https statethzchmailmanlistinforhelp please read posting guide http www rprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code macqueen lawrence livermore national laboratory livermore ca usa 9254231062 macqllnlgov ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code view message context httpwwwnabblecomrunningstrucchangetp18791489p18795602html sent r help mailing list archive nabblecom ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code brian ripley ripleystatsoxacuk professor applied statistics httpwwwstatsoxacukripley university oxford tel 44 1865 272861 self 1 south parks road 44 1865 272866 pa oxford ox1 3tg uk fax 44 1865 272595 ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code view message context httpwwwnabblecomrunningstrucchangetp18791489p18799748html sent r help mailing list archive nabblecom message 14 date sun 03 aug 2008 115921 0400 john p burkett subject r gmaps projection inset rhelprprojectorg messageid 4895d5d98030100uriedu contenttype textplain charsetiso88591 formatflowed running r version 261 gentoo linux im trying produce thematic map usa using gmaps package result thus far two problems point view first projection miller cylindrical elongates southern states flattens northern ones unattractively id prefer albers conic projection something similar second inset alaska placed south contiguous states id prefer put north would grateful anyone could suggest change projection andor b relocate inset relevant section r code pasted needed supply entire program data john code starts librarygmaps gridnewpage gridframenamemap statenames listalabama alaska arizona arkansas california colorado connecticut florida georgia hawaii idaho illinois indiana iowa kansas kentucky louisiana maine maryland massachusetts michigan minnesota mississippi missouri montana nebraska nevada new hampshire new jersey new mexico new york north carolina north dakota ohio oklahoma oregon pennsylvania rhode island south carolina south dakota tennessee texas utah vermont virginia washington west virginia wisconsin wyoming delaware delaware last missing value shaded map gridpackmapusalevelplotstatesstatenameslevelsspccolfungreens normalizetheightunit1null code ends john p burkett department environmental natural resource economics department economics university rhode island kingston ri 028810808 usa phone 401 8749195 message 15 date sun 03 aug 2008 170556 0100 patrick burns subject r access object inside function duncan murdoch cc rhelprprojectorg messageid 4895d7647010504pburnsseanetcom contenttype textplain charsetiso88591 formatflowed duncan murdoch wrote 03082008 1129 jonas wrote hi apologies covered numerous times lack search andor r skills stopping finding solution im trying access object defined created inside function object returned function cant seem get work keep getting object found error message thought solution found changing frameenvironment cant seem understand correctly want something like myfunction functionfrompackagexy plot objectinsidepackagefunction pseudocode doesnt really make sense didnt save result functionfrompackage anywhere change myfunction z functionfrompackagexy plot z perfectly standard question expecting r look find objectinsidepackagefunction suspect op doesnt care r looks long finds fact everything functionfrompackage except z sucked black hole time plot scene going make search impossible one solution create myfunctionfrompackage change return value include objectinsidepackagefunction patrick burns patrickburnsstatcom 44 020 8525 0696 httpwwwburnsstatcom home poetry guide unwilling user duncan murdoch im using r 271 linux sincerely jonas ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code message 16 date sun 3 aug 2008 181340 0200 jonas subject r access object inside function patrick burns cc rhelprprojectorg duncan murdoch messageid 1875de320808030913n372546bfq2c7e7a17800c59bcmailgmailcom contenttype textplain charsetiso88591 duncan patrick sorry clear nonworking code example im actually trying achieve creating wrapper around nelsonsiegel function fbonds package order among things create plot function however objects would like access inside myfunction returned nelsonsiegel function local inside nelsonsiegel function way inside wrapper function access objects myfunction z nelsonsiegelyield maturity doplotfalse would like access objects returned nelsonsiegel jonas sun aug 3 2008 605 pm patrick burns wrote duncan murdoch wrote 03082008 1129 jonas wrote hi apologies covered numerous times lack search andor r skills stopping finding solution im trying access object defined created inside function object returned function cant seem get work keep getting object found error message thought solution found changing frameenvironment cant seem understand correctly want something like myfunction functionfrompackagexy plot objectinsidepackagefunction pseudocode doesnt really make sense didnt save result functionfrompackage anywhere change myfunction z functionfrompackagexy plot z perfectly standard question expecting r look find objectinsidepackagefunction suspect op doesnt care r looks long finds fact everything functionfrompackage except z sucked black hole time plot scene going make search impossible one solution create myfunctionfrompackage change return value include objectinsidepackagefunction patrick burns patrickburnsstatcom 44 020 8525 0696 httpwwwburnsstatcom home poetry guide unwilling user duncan murdoch im using r 271 linux sincerely jonas ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code message 17 date sun 3 aug 2008 181954 0200 pir2jv subject r sans objet rhelprprojectorg messageid contenttype textplain charsetiso88591 delspyes formatflowed peuton poser des questions en franais trs simple je fais un programme trucr sous emacs ce programme ralise des graphiques je voudrais les enregistrer au fur et mesure sinon quand la session est finie ces graphiques sont perdus jai essay pdftrucpdf mais le fichier trucpdf cr est vide comment faiton merci et salutations jacques vernin message 18 date sun 3 aug 2008 172347 0100 bst prof brian ripley subject r drop1 seems give unexpected results compare anova thomas p c chu cc rhelprprojectorg messageid contenttype textplain charsetusascii formatflowed sun 3 aug 2008 thomas p c chu wrote thanks mr dalgaard advice everyone else contributed inclusion error term end simsety line cure problems drop1 step suppose inexperience carrying simulations caused gaffe r 280 warn asked compute f tests essentially 00 essentially relative total sum squares examples future easier notice brian ripley ripleystatsoxacuk professor applied statistics httpwwwstatsoxacukripley university oxford tel 44 1865 272861 self 1 south parks road 44 1865 272866 pa oxford ox1 3tg uk fax 44 1865 272595 message 19 date sun 3 aug 2008 162445 0000 dirk eddelbuettel subject r sans objet pir2jv cc rhelprprojectorg messageid 20080803162444ga16797masterdebianorg contenttype textplain charsetiso88591 sun aug 03 2008 061954pm 0200 pir2jv wrote peuton poser des questions en franais pas vraiment rhelp est une liste organisee en anglais trs simple je fais un programme trucr sous emacs ce programme ralise des graphiques je voudrais les enregistrer au fur et mesure sinon quand la session est finie ces graphiques sont perdus jai essay pdftrucpdf mais le fichier trucpdf cr est vide comment faiton cest une faq il te faut un devoff apres le ou les commandes plot lines donc pdftrucpdf plot110 typeb mainma graphique devoff merci et salutations de rien mais reviens en anglais silteplait dirk three two people difficulties fractions message 20 date sun 3 aug 2008 083047 0800 franz mueter subject r bubble plots cody hamilton messageid 000f01c8f586496b4a10dc41de30net contenttype textplain charsetusascii believe symbols looking course need convert categories numeric values plotting example adding actual data data frame cbindd p runif15 symbolsasnumericdtime asnumericdy circlesdp inches02 annf xaxtn yaxtn use axis get proper axis labels axis1 13 levelsdtime axis2 14 levelsdy franz original message rhelpbouncesrprojectorg mailtorhelpbouncesrprojectorg behalf cody hamilton sent friday august 01 2008 707 pm rhelprprojectorg subject r bubble plots way create bubble plot r example define following data frame containing level observed 5 patients three time points timecreptime 15reptime 25reptime 35 ycabcdabcadaaabcd ddataframecbindytime would like display percentage subjects level time point bubble whose size proportional percentage subjects given level given time point thus case data frame plot would levels abcd yaxis levels time time 1time 2 time 3 xaxis four bubbles time point eg size bubble bottom left corner plot would proportional percentage patients ya timetime 1 running r 271 windows regards cody ________________________________ message contains information may confidentdropped8 message 21 date sun 3 aug 2008 173801 0100 bst prof brian ripley subject r computer run 64 bit ubuntur tom la bone cc rhelprprojectorg messageid contenttype textplain formatflowed charsetusascii pemtium 940 supports emt64 according httpwwwpcstatscomarticleviewcfmarticleid1988 googled key emt64 intels name x86_64 box run x86_64 linux believe debianubuntu refer amd64 understandably intel favour gnu configure use linux r compiles readily x86_64 linux main development platform far 64 bit linux version r really run precisely able compile r sources also install amd64 ubuntu package r sun 3 aug 2008 tom la bone wrote reading 64bit systems software still somewhat uncertain things dell dimension xps 400 dual core intel pentium 940 32 ghz 4 gb memory currently dual boot system xp professional ubuntu 804 32 bit simply install ubuntu 804 64 bit system run 64 bit linux version r would appreciate advice might willing offer tom brian ripley ripleystatsoxacuk professor applied statistics httpwwwstatsoxacukripley university oxford tel 44 1865 272861 self 1 south parks road 44 1865 272866 pa oxford ox1 3tg uk fax 44 1865 272595 message 22 date sun 3 aug 2008 164331 0000 dirk eddelbuettel subject r computer run 64 bit ubuntur prof brian ripley cc tom la bone rhelprprojectorg messageid 20080803164331ga2156masterdebianorg contenttype textplain charsetusascii sun aug 03 2008 053801pm 0100 prof brian ripley wrote pemtium 940 supports emt64 according httpwwwpcstatscomarticleviewcfmarticleid1988 googled key emt64 intels name x86_64 box run x86_64 linux believe debianubuntu refer amd64 understandably intel favour gnu configure use linux r compiles readily x86_64 linux main development platform far 64 bit linux version r really run precisely able compile r sources also install amd64 ubuntu package r id rather grab current prebuild binaries cran mirror eg httpcranrprojectorgbinlinuxubuntu vincent michael really excellent job rebuilding debian packages ubuntu x86 amd64 typically within day goodies see readme aforementioned directory dirk three two people difficulties fractions message 23 date sun 03 aug 2008 124347 0400 duncan murdoch subject r access object inside function jonas cc rhelprprojectorg patrick burns messageid 4895e0433070909statsuwoca contenttype textplain charsetiso88591 formatflowed 03082008 1213 pm jonas wrote duncan patrick sorry clear nonworking code example im actually trying achieve creating wrapper around nelsonsiegel function fbonds package order among things create plot function however objects would like access inside myfunction returned nelsonsiegel function local inside nelsonsiegel function way inside wrapper function access objects myfunction z nelsonsiegelyield maturity doplotfalse would like access objects returned nelsonsiegel patrick said function returns locals z disappear need get source function modify duncan murdoch jonas sun aug 3 2008 605 pm patrick burns wrote duncan murdoch wrote 03082008 1129 jonas wrote hi apologies covered numerous times lack search andor r skills stopping finding solution im trying access object defined created inside function object returned function cant seem get work keep getting object found error message thought solution found changing frameenvironment cant seem understand correctly want something like myfunction functionfrompackagexy plot objectinsidepackagefunction pseudocode doesnt really make sense didnt save result functionfrompackage anywhere change myfunction z functionfrompackagexy plot z perfectly standard question expecting r look find objectinsidepackagefunction suspect op doesnt care r looks long finds fact everything functionfrompackage except z sucked black hole time plot scene going make search impossible one solution create myfunctionfrompackage change return value include objectinsidepackagefunction patrick burns patrickburnsstatcom 44 020 8525 0696 httpwwwburnsstatcom home poetry guide unwilling user duncan murdoch im using r 271 linux sincerely jonas ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code message 24 date sun 3 aug 2008 102527 0700 pdt michael milligan subject r missing f statistic anovagam rhelprprojectorg messageid 4185884927qmweb50305mailre2yahoocom contenttype textplain charsetutf8 hello encountered results sure interpret using anovagam compare 2 different models certain tests results include f associated pstatistic happens comparing certain models others discern patten explaining test works output tests x1 x2 1variable vectors x matrix several variables results compare model additively separable x1 x2 model assumed additively separable model 1 y1 sx1 sx2 x model 2 y1 x sx1 x2 resid df resid dev df deviance f prf 1 38155111 298606 2 38103577 298988 51534 382 f statistic computed though statistic computed dependent variables used results similar analysis different dependent variable model 1 y2 x1 x2 x model 2 y2 sx1 sx2 x resid df resid dev df deviance f prf 1 3822000 33970 2 3819535 33921 2465 49 22257 009578 signif codes 0 0001 001 005 01 1 anovaeqquadeqgamastestf analysis deviance table model 1 y2 x1 x1sq x2 x2sq x1x2 x model 2 y2 sx1 sx2 x resid df resid dev df deviance f prf 1 381900000 33955 2 381953502 33921 053502 34 model 1 y2 sx1 sx2 x model 2 y2 x sx1 x2 resid df resid dev df deviance f prf 1 38195350 33921 2 38218323 33967 22973 46 22391 009863 f statistic reported comparing linear model additively separable semiparametric model comparing additively separable model nonadditvely separable model comparing partially quadratic model xsq means x2 additively separable semiparametric model im happy provide information dataset estimation dont know might helpful really dont understand cause problem dataset small 3800 observations say x2 many observations value 0 appreciate light anyone shed issue thank much michael milligan phd candidate university new mexico message 25 date sun 3 aug 2008 182910 0100 bst prof brian ripley subject r computer run 64 bit ubuntur dirk eddelbuettel cc tom la bone rhelprprojectorg messageid contenttype textplain charsetusascii formatflowed sun 3 aug 2008 dirk eddelbuettel wrote sun aug 03 2008 053801pm 0100 prof brian ripley wrote pemtium 940 supports emt64 according httpwwwpcstatscomarticleviewcfmarticleid1988 googled key emt64 intels name x86_64 box run x86_64 linux believe debianubuntu refer amd64 understandably intel favour gnu configure use linux r compiles readily x86_64 linux main development platform far 64 bit linux version r really run precisely able compile r sources also install amd64 ubuntu package r id rather grab current prebuild binaries cran mirror eg httpcranrprojectorgbinlinuxubuntu vincent michael really excellent job rebuilding debian packages ubuntu x86 amd64 typically within day goodies see readme aforementioned directory meant amd64 ubuntu package r brian ripley ripleystatsoxacuk professor applied statistics httpwwwstatsoxacukripley university oxford tel 44 1865 272861 self 1 south parks road 44 1865 272866 pa oxford ox1 3tg uk fax 44 1865 272595 message 26 date sun 3 aug 2008 111001 0700 pdt cls59 subject r convert date decimal days rhelprprojectorg messageid 18800462posttalknabblecom contenttype textplain charsetusascii hello quick question formatting date strings currently debugging matlab code someone else wrote since bad go line line figured would rewrite thing r code produces plots wave spectra decimal days since epoch x axis wave period axis able convert date stamp input file posixct object following call date asposixctheader12tzgmtformatymd_k date 20080711 030000 gmt every thing looks groovy using asdoubledate gives number seconds since epoch reproduces correct date fed date u r command line however cant seem find clean way converting date decimal days since epoch using asdoubledate86400 gives 1407112 catastrophic loss significance number 14071125 tried playing around format command see could get r format date differently luck suggestions charlie view message context httpwwwnabblecomconvertdatetodecimaldaystp18800462p18800462h tml sent r help mailing list archive nabblecom message 27 date sun 3 aug 2008 112343 0700 charles c berry subject r convert loop apply billvenablescsiroau cc rhelprprojectorg anhnttranuclaedu messageid contenttype textplain charsetusascii formatflowed length levels a1cat large structure suggested toy example holds actual case probes a2 overlap order according a2st a2en might unsplitlapplylevelsa1cat functionx tmp1 subset a1 catx tmp2 subset a2 catx findinterval tmp1en tmp2st findinterval tmp1st tmp2en1 a1cat 1 1 1 2 2 0 1 hth chuck sun 3 aug 2008 billvenablescsiroau wrote way speed toy example _________________ a1 dataframeid 16 cat pastecat rep13 c231 st c1 7 30 40 59 91 en c5 25 39 55 70 120 a2 dataframeid pasteprobe 18 cat pastecat rep13 c233 st c1 9 20 38 53 70 80 95 en c6 15 36 43 58 75 85 98 c1 outera2st a1en c2 outera2en a1st c3 outera2cat a1cat a1coverage colsumsc1c2c3 __________________ wont work one step a1 30000 rows a2 200000 unless memory size approximately infinite need loop suppose handle 1000 probes time might able get away something like __________________ elided yahoo spam checkcoverage functiona1 a2 colsumsoutera2st a1en outera2en a1st outera2cat a1cat coverage numericn nrowa2 m2 0 whilem2 n m1 m2 1 m2 minm2 chunk n coveragem1m2 checkcoveragea1 a2m1m2 a1coverage coverage __________________ warning untested code failing go c code done btw think apply going useful bill venables csiro laboratories po box 120 cleveland 4163 australia office phone email preferred 61 7 3826 7251 fax absolutely necessary 61 7 3826 7304 mobile 61 4 8819 4402 home phone 61 7 3286 7700 mailtobillvenablescsiroau httpwwwcmiscsiroaubillvenables original message rhelpbouncesrprojectorg mailtorhelpbouncesrprojectorg behalf anh tran sent saturday 2 august 2008 404 pm rlist subject r convert loop apply hi alli know topic came multiple times ive never fully understand apply function anyway im asking help convert loop apply 2 data frames following information a1 fragment need covered a2 probes cover specific fragment need count number probes cover every given fragment need cat id fragment a1dataframeidc16 catccat 1cat 1cat 2cat 2cat 2cat 3 stc1730405991 enc525395570120 a2dataframeidpasteprobec18 catccat 1cat 1cat 2cat 2cat 2cat 3cat 3cat 3 stc19203853708095 enc615364358758598 a1coveragenull came loop basically probe starts fragment end end fragment start cover fragment 1lengtha1id a1coverageilengtha2a2sta1enia2ena1stia2cata1cat iid a1coverage 1 1 1 2 2 0 1 loop runs awefully slow 200000 probes 30000 fragments anyway speed apply time loop scan first 20 record dataset user system elapsed 2264 0501 2770 think room improvement idea thanks regards anh tran alternative html version deleted ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code charles c berry 858 5342098 dept familypreventive medicine e mailtocberrytajoucsdedu uc san diego httpfamprevmeducsdedufacultycberry la jolla san diego 920930901 message 28 date sun 3 aug 2008 205352 0100 bst prof brian ripley subject r convert date decimal days cls59 cc rhelprprojectorg messageid contenttype textplain charsetusascii formatflowed think find catastrophic loss significance printing default level try printx digits15 sun 3 aug 2008 cls59 wrote hello quick question formatting date strings datetime strings methinks want days use asdate asposixct currently debugging matlab code someone else wrote since bad go line line figured would rewrite thing r code produces plots wave spectra decimal days since epoch x axis wave period axis able convert date stamp input file posixct object following call date asposixctheader12tzgmtformatymd_k date 20080711 030000 gmt every thing looks groovy using asdoubledate gives number seconds since epoch reproduces correct date fed date u r command line however cant seem find clean way converting date decimal days since epoch using asdoubledate86400 gives 1407112 catastrophic loss significance number 14071125 tried playing around format command see could get r format date differently luck suggestions charlie view message context httpwwwnabblecomconvertdatetodecimaldaystp18800462p18800462h tml sent r help mailing list archive nabblecom ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code brian ripley ripleystatsoxacuk professor applied statistics httpwwwstatsoxacukripley university oxford tel 44 1865 272861 self 1 south parks road 44 1865 272866 pa oxford ox1 3tg uk fax 44 1865 272595 message 29 date sun 3 aug 2008 131906 0700 pdt cls59 subject r convert date decimal days rhelprprojectorg messageid 18802117posttalknabblecom contenttype textplain charsetusascii youre right problem print function displaying enough precision thanks view message context httpwwwnabblecomconvertdatetodecimaldaystp18800462p18802117h tml sent r help mailing list archive nabblecom message 30 date sun 3 aug 2008 135507 0700 subject r determining model parameters rhelprprojectorg messageid 20080803165507gxq3d82014rootmp17 contenttype textplain charsetutf8 may begining question please bear data based historgram plots looks like beta distribution function functions within r help determine model parameters distirbution similarily common distirbutions poissonlambda chisquaredegrees freedom chisquare statistic etc normal distribution relatively straight forward mean variance computed failrly readily far reading parameters less given kind inference using prior model seems key parameters least estimated accurately thank kevin message 31 date sun 3 aug 2008 210631 0000 utc ben bolker subject r determining model parameters rhelpstatmathethzch messageid contenttype textplain charsetusascii charternet writes data based historgram plots looks like beta distribution function functions within r help determine model parameters distirbution librarymass fitdistr message 32 date sun 3 aug 2008 165443 0400 andrew ramsey subject r changing values rhelprprojectorg messageid 58c18421686a4d5ea95229d547e07434schoolphumassedu contenttype textplain charsetusascii delspyes formatflowed hello relatively new user r find information need please help large data set values including letters numbers symbols sometimes within vector value ie x9 ive imported data using readfwp arrives list format id like change letters symbols numbers ie x9 00911 every entry would recommend try thank message 33 date sun 3 aug 2008 175517 0400 jim holtman subject r changing values andrew ramsey cc rhelprprojectorg messageid 644e1f320808031455s3b08c840v3338774f06aed855mailgmailcom contenttype textplain charsetiso88591 exactly translation x9 00911 done unique mappings character sequences numbers many different ones might leading zeros result changed numeric default printing results 911 provide idea data looks like eg stryourdata done selected columns need provide bit background understand problem trying solve sun aug 3 2008 454 pm andrew ramsey wrote hello relatively new user r find information need please help large data set values including letters numbers symbols sometimes within vector value ie x9 ive imported data using readfwp arrives list format id like change letters symbols numbers ie x9 00911 every entry would recommend try thank ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code jim holtman cincinnati oh 1 513 646 9390 problem trying solve message 34 date sun 3 aug 2008 231346 0000 utc michael bibo subject r bubble plots rhelpstatmathethzch messageid contenttype textplain charsetusascii cody hamilton edwardscom writes way create bubble plot r example define following data frame containing level observed 5 patients three time points timecreptime 15reptime 25reptime 35 ycabcdabcadaaabcd ddataframecbindytime would like display percentage subjects level time point bubble whose size proportional percentage subjects given level given time point thus case data frame plot would levels abcd yaxis levels time time 1time 2 time 3 xaxis four bubbles time point eg size bubble bottom left corner plot would proportional percentage patients ya timetime 1 sounds like function balloonplot package gplots might looking however bug labelling plot meaning contact maintainer change lines 364 376 respectively labelsynames labelsylab labelsxnames labelsxlab work specific purpose included full code amended balloonplot function end email copy paste r console required given dataframe attachd balloonplotproptabletabletimey100 gives percentages play proptable syntax get percentages row column totals play titles margin totals etc amended balloonplot function allpurpose fix works specific purpose id balloonplotr 908 20060302 214324z warnes balloonplot functionx usemethodballoonplotx balloonplottable functionx xlab ylab zlab showzeros false showmargins true obj x tmp asdataframex x tmp1 tmp2 z tmp3 tableflag true ifmissingxlab xlab namesdimnamesobj1 ifmissingylab ylab namesdimnamesobj2 ifmissingzlab zlab freq balloonplotdefaultx z xlabxlab ylabylab zlabzlab showzeros showzeros showmargins showmargins balloonplotdefault functionxyz xlab ylab zlabdeparsesubstitutez dotsize2maxstrwidth19strheight19 dotchar19 dotcolorskyblue main labeltrue labeldigits2 scalemethodcvolumediameter colsrtparsrt rowsrtparsrt colmar1 rowmar2 showzerosfalse showmarginstrue cummarginstrue sortedtrue labellinestrue funfunctionxsumxnarmt hideduplicatestrue ifisnullnamesx xnames ascharactersubstitutex iflengthxnames1 xnames xnames1 else xnames namesx ifisnullnamesy ynames ascharactersubstitutey iflengthynames1 ynames ynames1 else ynames namesy ifmissingxlab xlab paste xnames collapse ifmissingylab ylab paste ynames collapse handle arguments scalemethod matchargscalemethod anyz 0 warningz values zero detected balloons displayed cells ifmissingmain ifscalemethodvolume main pasteballoon plot pastexnames collapse pasteynames collapse narea proportional zlab sep else main pasteballoon plot pasteynames collapse pasteynames collapse ndiameter proportional zlab sep iflengthdotcolor subject r determining model parameters rkevinburtoncharternet cc rhelprprojectorg messageid 48964faa4020203ufledu contenttype textplain charsetutf8 formatflowed beta distribution applies data bounded 0 1 cases strictly bounded ie 0 wrote charternet writes data based historgram plots looks like beta distribution function functions within r help determine model parameters distirbution librarymass fitdistr message 36 date sun 3 aug 2008 200015 0700 j dougherty subject r changing values rhelprprojectorg messageid 20080803200015726jwdsurewestnet contenttype textplain charsetiso88591 sunday 03 august 2008 015443 pm andrew ramsey wrote hello relatively new user r find information need please help large data set values including letters numbers symbols sometimes within vector value ie x9 ive imported data using readfwp arrives list format id like change letters symbols numbers ie x9 00911 every entry would recommend try thank andrew information little sparse limits helpful list members purpose recoding values readfwp mean readfwf fixedwdith format looking example leading zeros example value would seem mean value still alphanumeric string rather number per se might order neatly printout common ascii really different still list rule would recommend handling tasks related data coding dbms postgresql access afterall john message 37 date mon 4 aug 2008 145445 1000 sornnorngdpivicgovau subject r sweave ggplot2 rhelprprojectorg messageid contenttype textplain hi ive trying run sweave r code embedded using ggplot2 package particular qplot command appears problem sweave picking qplot function anybody else tried use qplot sweave successful help would much appreciated kind regards sorn notice email attachments may contain information tdropped14 message 38 date mon 4 aug 2008 163214 1000 andrew robinson subject r decomposing tests interaction terms mixedeffects models rhelp discussion messageid 20080804063214gn3754msunimelbeduau contenttype textplain charsetusascii dear r colleagues friend trying develop modest workflow problem decomposing tests higherorder terms interpretable sets tests lower order terms conditioning example interaction 3 levels c 2 levels significant may interest ask whether significant level 1 c level 2 c textbook response seems subset data perform tests two subsets using rss df original fit trying answer question using new explanatory variables approach seems work fine aov lme example build example dataset setseed101 block gl6 6 36 lab pasteb 16 sep gl3 2 36 lab pastea 13 sep c gl2 1 36 lab pastec 12 sep example dataframeblock c exampley rnorm362 asnumericexampleaasnumericexamplec2 3 reprnorm6 each6 withexample interactionplota c lme requirenlme anovalmeyac random 1block data example summaryaovy errorblock ac data example significant 2way interaction would like test effect c1 effect c2 construct two new variables decompose interaction lrt possible exampleac exampleac1 exampleac2 interactionexamplea examplec exampleac1examplec c1 a1c1 constant c1 exampleac2examplec c2 a1c2 constant c2 lme fails lmer lmey ac1 ac random 1 block data example aov seems fine summaryaovy errorblock ac1 ac data example ac significant significant c1 summaryaovy errorblock ac2 ac data example ac significant significant c2 experimentation suggests lme doesnt like partial confounding approach using rather variables constructed lmey ac random 1 block data example lmey ac random 1 block data example anything obviously wrong another approach goal trying achieve many thanks andrew andrew robinson department mathematics statistics tel 61383446410 university melbourne vic 3010 australia fax 61383444599 httpwwwmsunimelbeduauandrewpr httpblogsmbsedufishinginthebay message 39 date mon 4 aug 2008 161928 0930 fernando marmolejo ramos subject r 95ci around median rhelprprojectorg messageid 12178325684896a678c5c5ewebmailadelaideeduau contenttype textplain charsetiso88591 dear people ive learnt using boxplotstats command grdevices library get 5number summaries boxplot plus important information like confidence interval around median im interested knowing actual formula used package calculate confidence interval someone help cheers fernando message 40 date mon 04 aug 2008 170009 1000 simon blomberg subject r 95ci around median fernando marmolejo ramos cc rhelp messageid 12178332091771340camelsibsblomber01dsibuqeduau contenttype textplain see fivenum stats package type statsfivenum get code crucial calculations last lines simon mon 20080804 1619 0930 fernando marmolejo ramos wrote dear people ive learnt using boxplotstats command grdevices library get 5number summaries boxplot plus important information like confidence interval around median im interested knowing actual formula used package calculate confidence interval someone help cheers fernando ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code simon blomberg bsc hons phd mappstat lecturer consultant statistician faculty biological chemical sciences university queensland st lucia queensland 4072 australia room 320 goddard building 8 61 7 3365 2506 httpwwwuqeduauuqsblomb email sblomberg1_at_uqeduau policies 1 analyse data 2 deadline problem combination data aching desire answer ensure reasonable answer extracted given body data john tukey message 41 date mon 4 aug 2008 092640 0200 morneau franois subject r graph elyakhlifi mustapha cc rhelprprojectorg messageid 5e3d22a4869bb94aa1138ab97660d8b6013cf672populusifnfr contenttype textplain charsetiso88591 hello segments looking franois message dorigine de rhelpbouncesrprojectorg mailtorhelpbouncesrprojectorg de la part de elyakhlifi mustapha envoy vendredi 1 aot 2008 1715 rhelprprojectorg objet r graph hello dont know ouput segments points like chart attached help please thanks ________________________________________________________________________ ____ yahoofr message 42 date mon 4 aug 2008 095602 0200 draga r subject r major difference multivariate analyses spss r messageid contenttype textplain dear colleagues know spss compute linear mixed models used r computing multivariate analyses never encountered major difference outcome spss r spss pearson correlation variable 1 variable 2 31 p0001 spss binary logistic regression gives us odds ratio or49 95 ci 2790 p0001 n338 lower upper gender 1120 0565 2221 age 0985 0956 1015 variable 2 4937 2698 9032 r multilevel logistic regression using statistical package lmer gives us odds ratio102 95 ci 6314 p024 n338 groups group 1 98 group 2 84 lower upper gender 2295 2840 7430 age 0003 70047 70054 variable 2 10176 6295 14056 crosstabs gives us variable var b 0 1 0 156 108 1 17 57 would somebody know possible spss get p0001 r get p024 r 95 ci odds ratio 62141 pvalue024 thanks ronald alternative html version deleted message 43 date mon 4 aug 2008 100234 0200 onkelinx thierry subject r sweave ggplot2 messageid 2e9c414912813e4eb981326983e0a1040527d64binboexchinbobe contenttype textplain charsetusascii dear sorn hard guess problem dont provide sample code guess graphics empty use printqplot qplot latter wont work need printqplot hth thierry ir thierry onkelinx instituut voor natuur en bosonderzoek research institute nature forest cel biometrie methodologie en kwaliteitszorg section biometrics methodology quality assurance gaverstraat 4 9500 geraardsbergen belgium tel 32 54436 185 thierryonkelinxinbobe wwwinbobe call statistician experiment done may asking perform postmortem examination may able say experiment died sir ronald aylmer fisher plural anecdote data roger brinner combination data aching desire answer ensure reasonable answer extracted given body data john tukey oorspronkelijk bericht van rhelpbouncesrprojectorg mailtorhelpbouncesrprojectorg namens sornnorngdpivicgovau verzonden maandag 4 augustus 2008 655 aan rhelprprojectorg onderwerp r sweave ggplot2 hi ive trying run sweave r code embedded using ggplot2 package particular qplot command appears problem sweave picking qplot function anybody else tried use qplot sweave successful help would much appreciated kind regards sorn notice email attachments may contain information dropped10 message 44 date mon 04 aug 2008 101738 0200 peter dalgaard subject r decomposing tests interaction terms mixedeffects models andrew robinson cc rhelp discussion messageid 4896bb2270302biostatkudk contenttype textplain charsetiso88591 formatflowed andrew robinson wrote dear r colleagues friend trying develop modest workflow problem decomposing tests higherorder terms interpretable sets tests lower order terms conditioning example interaction 3 levels c 2 levels significant may interest ask whether significant level 1 c level 2 c textbook response seems subset data perform tests two subsets using rss df original fit trying answer question using new explanatory variables approach seems work fine aov lme example build example dataset setseed101 block gl6 6 36 lab pasteb 16 sep gl3 2 36 lab pastea 13 sep c gl2 1 36 lab pastec 12 sep example dataframeblock c exampley rnorm362 asnumericexampleaasnumericexamplec2 3 reprnorm6 each6 withexample interactionplota c lme requirenlme anovalmeyac random 1block data example summaryaovy errorblock ac data example significant 2way interaction would like test effect c1 effect c2 construct two new variables decompose interaction lrt possible exampleac exampleac1 exampleac2 interactionexamplea examplec exampleac1examplec c1 a1c1 constant c1 exampleac2examplec c2 a1c2 constant c2 lme fails lmer lmey ac1 ac random 1 block data example aov seems fine summaryaovy errorblock ac1 ac data example ac significant significant c1 summaryaovy errorblock ac2 ac data example ac significant significant c2 experimentation suggests lme doesnt like partial confounding approach using rather variables constructed lmey ac random 1 block data example lmey ac random 1 block data example anything obviously wrong another approach goal trying achieve looks like generic issue lmelmer happy rank deficiencies design matrix heres dirty trick really supposed know hidden inside stats namespace modelmatrixac1ac dataexample m2 statsthincolm summarylmey m2 1 random 1 block data example thincol thinrow rank support functions anovamlm et al maybe document put open o__ peter dalgaard ster farimagsgade 5 entrb c _ dept biostatistics po box 2099 1014 cph k university copenhagen denmark ph 45 35327918 pdalgaardbiostatkudk fax 45 35327907 message 45 date mon 04 aug 2008 104457 0200 peter dalgaard subject r major difference multivariate analyses spss r draga r cc rhelprprojectorg messageid 4896c1897080501biostatkudk contenttype textplain charsetiso88591 formatflowed draga r wrote would somebody know possible spss get p0001 r get p024 r 95 ci odds ratio 62141 pvalue024 asked still providing sufficient information harold already pointed excluding important variance component spss could easily explain difference p values however discrepancy ci p value suggests could misreading output lmer andor misspecifying model simply telling us data like esp variables vary within levels output version lmer used also confusing us variable 1 variable 2 variable var b o__ peter dalgaard ster farimagsgade 5 entrb c _ dept biostatistics po box 2099 1014 cph k university copenhagen denmark ph 45 35327918 pdalgaardbiostatkudk fax 45 35327907 message 46 date mon 04 aug 2008 095725 0100 gavin simpson subject r 95ci around median simon blomberg cc rhelp messageid 1217840246265510camelgraptoleberisgeoguclacuk contenttype textplain mon 20080804 1700 1000 simon blomberg wrote see fivenum stats package type statsfivenum get code crucial calculations last lines give code calculate five number summary fernando wants know confidence interval calculated boxplotstats see code type boxplotstats followed return command line r relevant line conf doconf stats3 c158 158 iqrsqrtn working median stats3 details computation boxplotstats fernandos first port call two works cited support calculation full references references section help page hth g simon mon 20080804 1619 0930 fernando marmolejo ramos wrote dear people ive learnt using boxplotstats command grdevices library get 5number summaries boxplot plus important information like confidence interval around median im interested knowing actual formula used package calculate confidence interval someone help cheers fernando ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code message 47 date mon 4 aug 2008 095814 0100 keith jewell subject r unexpected nls behaviour solved rhelpstatmathethzch messageid hi everyone id omitted nonoptional parameters argument selfstart making change ssbatch gives successful result two calls nls ssbatchselfstart modelfunctionbatch coeffs coeffsbatch initialfunctionmcall data lhs estimate coefficients mean batch xy sortedxydatamcallbatch lhs data batch dataascharactermcallbatch check batch successive integers starting 1 minxyx 1 anydiffxyx1 stop batch successive integers sequence lval listxyy nameslval mcallcoeffs lval parameters ccoeffs sorry wasting anyones time keith jewell keith jewell wrote message news hi everyone thought selfstart function two exactly equivalent nlsaform df nlsaform df startgetinitialaform df example case r although splus v62 ssbatchselfstart modelfunctionbatch coeffs coeffsbatch initialfunctionmcall data lhs estimate coefficients mean batch xy sortedxydatamcallbatch lhs data batch dataascharactermcallbatch check batch successive integers starting 1 minxyx 1 anydiffxyx1 stop batch successive integers sequence lval listxyy nameslval mcallcoeffs lval df dataframeac09 11 19 20 21 29 30 batchc1122233 aform formulaassbatchbatchca nlsaform df startgetinitialaform df nlsaform df dont ask id want silly selfstart thats long story guess wherever would used nlsaform df could use nlsaform df startgetinitialaform df seems clumsy anyone point mistake limitation nls r hesitate use bg word thanks advance keith jewell dont think relevant completeness sessioninfo version 270 20080422 i386pcmingw32 locale lc_collateenglish_united kingdom1252lc_ctypeenglish_united kingdom1252lc_monetaryenglish_united kingdom1252lc_numericclc_timeenglish_united kingdom1252 attached base packages 1 stats graphics grdevices datasets tcltk utils methods base attached packages 1 xlsreadwrite_132 svsocket_095 svio_095 r2html_158 svmisc_095 svide_095 loaded via namespace attached 1 tools_270 vgam_077 message 48 date mon 4 aug 2008 024044 0700 pdt pierre8rou subject r format output file way want rhelprprojectorg messageid 18808276posttalknabblecom contenttype textplain charsetusascii think found answer something better please write pierre8r r code libraryquantmod libraryxts lines 2008070102001576115766157601576365 2008070102151576215765157571576195 2008070102301576215765157581575958 2008070102451575815758157451574691 20080701030015745157511574415744117 20080701031515742157421572715729100 2008070103301573015736157301573661 2008070103451573515740157351573955 quotes readcsvtextconnectionlines headerfalse x asxtsquotes12 asposixctpastequotes1quotes2formatymd hm colnamesx copenhighlowclosevolume xopen sprintf504f xopen xhigh sprintf504f xhigh xlow sprintf504f xlow xclose sprintf504f xclose xvolume sprintf0f xvolume writetablexfile koutputfilecsv quotefalse colnamesfalse rownamesformatindexxymdhm sep view message context httpwwwnabblecomhowtoformattheoutputfilejustthewayiwant tp18790926p18808276html sent r help mailing list archive nabblecom message 49 date mon 4 aug 2008 115053 0200 martin maechler subject r source script file straight subversion repository steven mckinney cc marc_schwartzcomcastnet rhelprprojectorg murdochstatsuwoca messageid 1858253501790713250345statmathethzch contenttype textplain charsetusascii sm steven mckinney fri 1 aug 2008 183854 0700 writes sm thanks duncan murdoch marc schwartz sm excellent help sm dont yet apache http interface svn sm running yet svn export access mechanism yes note svnrprojectorg r sources recommended packages ones runs apache secure http interface https security reasons big advantage exist user logins machine hence compromised user ssh information misused break elided yahoo spam martin maechler eth zurich _______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code end rhelp digest vol 66 issue 4 ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code tue 05 aug 2008 234404 0100 r p values non linear regression singular gradients using nls