Return to Rnumerical.problem CVS log

Up to [LON-CAPA] / loncom / homework / templates

Annotation of loncom/homework/templates/Rnumerical.problem, revision 1.3

1.1       www         1: <problem>
1.3     ! bisitz      2: 
1.1       www         3: <script type="loncapa/perl">
                      4: $seed=&random(1,500,1);
                      5: $n=&random(15,25,1);
                      6: $offset=&random(2,5,0.1);
                      7: $slope=&random(0.6,2.5,0.1);
                      8: # construct a data set using R
1.2       www         9: # dump is for debugging, print to screen to see data structure
                     10: ($data,$dump)=&cas_hashref('R',"set.seed($seed);x<-1:$n;w<-1+sqrt(x)/2;data.frame(x=x,y=$offset+$slope*x+rnorm(x)*w);");
                     11: @x=&cas_hashref_array($data,'x');
                     12: @y=&cas_hashref_array($data,'y');
1.1       www        13: $datax=join(',',@x);
                     14: $datay=join(',',@y);
                     15: # calculate the right answer using R
1.2       www        16: ($answerdata,$dump)=&cas_hashref('R',"x<-c($datax);y<-c($datay);dataset<-data.frame(x=x,y=y);fm<-lm(y~x,data=dataset);");
                     17: $answer=&cas_hashref_entry($answerdata,'coefficients','x').'*x+'.&cas_hashref_entry($answerdata,'coefficients','(Intercept)');
1.1       www        18: </script>
1.2       www        19: 
                     20: <startouttext />
                     21: Consider the plotted data set.<br />
                     22: <endouttext />
1.3     ! bisitz     23: 
1.2       www        24: <gnuplot width="400" solid="0" plotcolor="monochrome" gridlayer="off" bmargin="default" font="9" alttag="dynamically generated plot" bgcolor="xffffff" texfont="22" transparent="off" plottype="Cartesian" rmargin="default" gridtype="Cartesian" minor_ticscale="0.5" fontface="sans-serif" grid="on" align="left" texwidth="93" height="300" border="on" samples="100" fgcolor="x000000" major_ticscale="1" tmargin="default" lmargin="default" fillstyle="empty">
1.1       www        25:     <curve linestyle="points" linetype="solid" color="x000000" pointtype="3" limit="closed" pointsize="2" linewidth="1">
                     26:         <data>@x</data>
                     27:         <data>@y</data>
                     28:     </curve>
                     29: </gnuplot>
                     30: 
                     31: <startouttext />
1.2       www        32: <br />
1.1       www        33: Give a linear function approximating the data.<br />
                     34: <tt>y(x)=</tt>
                     35: <endouttext />
1.3     ! bisitz     36: 
1.2       www        37: <formularesponse id="11" answer="$answer" samples="x@1:$n#20">
1.1       www        38:     <responseparam name="tol" default="2%" description="Numerical Tolerance" type="tolerance" />
                     39:     <textline size="25" readonly="no" />
                     40: </formularesponse>
1.3     ! bisitz     41: 
1.1       www        42: </problem>